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1.  INTRODUCTION 


The  function  of  the  round  discrimination  system  (RDS)  is  to  improve  the  scoring  accuracy  of  the 
remote  equipment  target  system  (RETS)  currently  being  used  for  live-fire  training  at  training  centers 
worldwide.  The  RDS,  which  interfaces  directly  with  the  RETS,  discriminates  among  different  main  gun 
kinetic  energy  (KE),  high-explosive  antitank  (HEAT),  and  high-explosive  incendiary  (HEI)  and  machine 
gun  (ball  and  tracer)  training  practice  (TP)  ammunition  from  the  following  armored  vehicles:  all  Ml 
Abrams  series  main  battle  tanks,  the  M60A3  main  battle  tank,  and  the  M2/3  Bradley  Rghting  Vehicle 
(BFV).  This  system  can  determine  whether  the  incoming  projectile  could  produce  a  kill  during  live  fire 
exercises. 

The  RDS  program  was  conducted  in  three  phases.  The  first  phase  was  a  feasibility  test.  Acoustic 
signatures  of  medium-  to  large-caliber  projectiles  from  the  three  types  of  armored  vehicles  were  collected; 
these  signatures  were  found  to  have  distinct  acoustic  features,  allowing  the  program  to  continue.  In  the 
second  phase,  acoustic  signatures  of  the  35-mm  tank  precision  gunnery  in-bore  device  (TPGID)  rounds 
and  small-caliber  machine  gun  rounds  were  collected  so  that  the  database  could  be  completed.  An 
algorithm  and  a  brassboard  were  developed  on  the  acoustic  signature  discriminants  found.  The  major 
discriminants  used  include  the  duration  of  the  supersonic  shock  wave,  the  location  of  the  round  relative 
to  the  acoustic  sensors,  and  the  terminal  velocity  of  the  round.  In  the  third  phase,  verificationA'alidation, 
a  field  test  was  conducted  at  Aberdeen  Proving  Ground  (APG),  MD.  A  working  brassboard  correctly 
discriminated  90%  of  the  rounds  fired.  Finally,  the  RDS  brassboard  was  demonstrated  as  a  proof-of- 
concept  to  the  sponsor,  the  Seventh  Anny  Training  Command  (7th  ATC)  in  Grafenwoehr,  Germany. 

This  report  summarizes  the  background  material,  the  round  discrimination  algorithm,  and  the  current 
RDS  brassboard  design.  The  test  data,  software  development,  and  system  components  are  mentioned  in 
the  discussioa 

2.  BACKGROUND 

2.1  RDS  Program.  The  RDS  program  was  initiated  by  the  Army  Materiel  Command  Field  Assistance 
in  Science  and  Technology  Junior  (AMC-FAST  Jr.)  engineers  at  the  U.S.  Araiy  Research  Laboratory 
(ARL),  in  response  to  a  request  by  the  AMC-FAST  Science  Advisor  at  the  7th  ATC  in  1992.  The  request 
from  the  7th  ATC  reads  as  follows:  "The  7th  ATC  requires  a  target  RDS,  which  identifies  the  round 
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location  within  an  area  which  extends  1  m  beyond  the  target  edge  and  classifies  the  round  type  into  four 
band  size  categories;  <20  mm,  20-30  mm,  30-90  mm,  and  >90  mm."  At  the  very  least,  the  Commander 
of  the  7th  ATC  wanted  a  system  able  to  discriminate  between  projectiles  <30  mm  (small  caliber  and  BFV 
main  gun  projectiles)  and  >30  mm  (Ml  A1  main  gun  projectiles)  with  a  per  umt  (system)  cost  of  around 
$5k.  Table  1  identifies  the  projectile’s  group  by  band  and  caliber.  Several  sensor  technologies  were 
explored  as  potential  candidates,  including  acoustics,  infiared,  radar,  and  video.  The  acoustic  approach 
was  chosen  for  two  reasons:  cost  and  the  ability  to  exploit  the  uniqueness  of  acoustic  signatures  of 
armored  vehicle  rotmds. 


Table  1.  Projectiles  Tested 


Projectile 

Cdiber 

(In-bore) 

(mm) 

Projectile 

C^ber 

(In-flight) 

(mm) 

Projectile 

Nomenclature 

Projectile 

Description 

Projectile 
Group  by 
Band 
(mm) 

Projectile 

Group  by 

Caliber 

(mm) 

7.62 

7.82 

M59 

Ball 

<20 

7.62 

7.62 

7.82 

M62 

Tracer 

<20 

7.62 

12.7 

12.95 

M33 

Ball 

<20 

12.7 

12.7 

12.95 

M17 

Tracer 

<20 

12.7 

25 

16.2 

M910 

TPDS-T 

<20 

25  KE 

25 

25 

M793 

TP-T 

20-30 

25HEI 

35 

35 

XM? 

TPL-T 

30-90 

35  TPGID 

35 

35 

DM18 

HEAT-TP-T 

30-90 

35TPGID 

35 

35 

DM68 

HEAT-TP-T 

30-90 

35  TPGID 

105 

63.5 

M724 

TPDS-T 

30-90 

105/120  KE 

120 

38 

M865 

TPCSDS-T 

30-90 

105/120  KE 

105 

105 

M490 

TP-T 

>90 

105/120  HE 

120 

120 

M831 

HEAT-TP-T 

>90 

105/120  HE 

2.2  RETS.  The  RETS,  as  seen  in  Figure  1,  uses  vibration  sensors  (also  called  hit  sensors)  that  are 
attached  to  the  target  and  connect  to  the  RETS  electronic  control  unit  (ECU).  When  the  sensors  detect 
a  hit,  the  ECU  conveys  a  hit  signal  to  a  range  control  tower  via  a  hardwired  or  radio  link.  This  activates 
an  electromechanical  mechanism  that  lowers  the  plywood  tank  target,  thus  indicating  a  kiU.  A  deficiency 
of  the  current  system  is  the  absence  of  a  mechanism  to  discriminate  between  rounds,  locate  the  rounds, 
and  activate  the  target  only  in  response  to  those  rounds  that  can  produce  a  kill.  A  12.7-mm  (.50  cal.) 
machine  gun  round,  for  example,  can  currently  be  credited  with  killing  a  T-72  tank  target.  RETS  can  be 
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Figure  1.  Remote  equipment  target  system. 


fooled  by  debris  from  ricocheting  rounds  or  by  small-caliber  machine  gun  fire  onto  the  target,  simulating 
a  main  gun  round  vibration  pattern.  As  a  result,  armor  kills  by  non-armor-defeating  projectiles 
erroneously  inflate  the  training  unit’s  score.  In  addition,  this  deficiency  prohibits  an  accurate  after-action 
evaluation  of  the  trainee’s  ability  to  choose  the  proper  ammunition  for  the  mission  at  hand.  The  RDS 
eliminates  this  deficiency  and  increases  the  accuracy  of  the  training  process. 

3.  RDS  CONCEPT 

An  extensive  literature  search  of  existing  scoring  systems  and  potential  signature-producing 
technologies  was  undertaken.  The  techniques  investigated  included  using  video  cameras,  colored  tracers, 
lasers,  radar,  and  acoustics.  Many  were  dismissed  due  to  cost,  practicality,  or  lack  of  technology.  The 
two  most  promising  concepts  recommended  by  the  FAST  Jr.  team  were  analysis  of  acoustical  signatures 
and  radar  signatures  of  the  passing  projectile.  Analysis  of  the  radar  signature  approach  was  subsequently 
deemed  not  feasible  due  to  high  cost.  It  was  determined  that  the  simplest  and  most  cost-effective  method 
of  discrimination  would  be  an  acoustic  approach. 

In  order  to  explore  the  feasibility  of  an  acoustic  RDS  system,  a  three-test  program  was  proposed  that 
included  participants  from  three  ARL  directorates — ^the  Weapons  Technology  Directorate  (WTD),  the 
Survivability/Lethality  Analysis  Directorate  (SLAD),  and  the  Sensors,  Signatures,  Signal  and  Information 
(SSI)  Processing  Directorate.  A  test  matrix  was  developed  to  include  all  possible  rounds  and  engagement 
distances  used  on  a  live-fire  gunnery  range  during  a  Combined  Arms  Live-Fire  Exercise  (CALFEX),  as 
dictated  by  FM  17-12-1  and  FM  23-1.  The  first  part  was  a  feasibility  test,  where  acoustical  signature  data 


were  collected  and  analyzed.  The  test  data  were  used  to  determine  if  acoustic  signal  discrimination 
between  projectile  data  at  the  extreme  firing  ranges  of  the  7th  ATC  interest  array  were  possible.  Tlie 
second  part  of  the  program  extended  the  data  collection  effort  for  projectile  discrimination  by  collecting 
data  on  projectiles  within  the  extremes  investigated  in  part  1.  During  the  third  part,  a  brassboard  RDS 
was  constmcted  and  tested  at  APG.  This  test  included  all  the  projectiles  used  on  the  live-fire  gunnery 
range.  Subsequent  to  verification  of  the  brassboard  system’s  ability  to  discriminate  effectively  and 
efficiently,  the  system  was  to  be  demonstrated  at  the  7th  ATC. 

The  projectiles  fired  during  parts  1  and  2  of  the  test  program  are  those  actually  fired  on  live-fire 
gunnery  ranges  and  are  listed  in  Table  1.  The  projectile’s  in-flight  caliber  determined  its  band  group.  No 
large-caliber  service  rounds  were  considered  in  the  test  program  due  to  cost  and  because  they  are  not  used 
in  live-fire  training  exercises.  However,  large-caliber  service  rounds  could  be  added  to  the  acoustic 
database  at  a  later  date. 

3.1  Acoustic  Round  Discrimination  Technique. 

3.1.1  Background.  As  a  supersonic  projectile  travels,  several  shocks  from  the  projectile  are 
propagated  to  the  surrounding  air.  These  shocks  are  created  from  the  projectile’s  bow,  surface  features 
on  the  projectile’s  body,  and  the  tail.  A  sensitive  air  pressure  gauge  placed  in  ttie  path  of  the  shock  waves 
will  register  what  is  typically  referred  to  as  the  N-Wave.  The  characteristics  of  the  N-Wave  are  directly 
dependent  on  the  geometry  of  the  projectile,  the  projectile’s  velocity,  the  location  of  the  shock  origination 
relative  to  the  pressure  gauge,  and  the  condition  of  the  ambient  air.  Other  discontinuous  features  of  the 
projectile,  such  as  fins  or  threads,  will  contribute  to  additional  characteristic  features  on  the  N-Wave. 
Figure  2  graphically  depicts  the  shock  cone  formed  from  supersonic  projectiles,  while  Figure  3  shows  the 
characteristic  N-Wave  pressure  profile. 

3.1.2  Theory  and  Formulation.  A  body  traveling  through  the  air  faster  than  the  ambient  speed  of 
sound  will  generate  a  conic  shock  profile.  The  shock  pressure  profile  shape  is  directly  dependent  on  the 
object’s  form.  Air  is  compressed  against  the  leading  part  of  the  object.  Since  the  object  is  traveling  faster 
than  the  time  required  for  the  air  to  re-expand  ahead  of  the  object,  a  conic  shock  front  is  generated.  A 
sensitive  pressure  gauge  some  distance  away  from  the  object’s  trajectory  will  register  the  lead  shock  as 
a  time-pressure  discontinuity.  As  the  body  passes,  the  compressed  air  expands,  and  the  local  pressure 
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decreases.  As  surface  features  on  the  object  create  shocks,  the  profile  sensed  by  the  pressure  gauge 
register  these  as  small  positive  discontinuities.  As  the  local  pressure  decreases  towards  the  tail  of  the 
object,  the  local  pressure  reaches  a  neutral  pressure  area  that  is  equivalent  to  the  ambient  pressure. 
Pressure  aft  of  the  neutral  area  continues  to  decrease  below  ambient.  If  the  object  has  surface  features 
in  the  aft  region,  the  shock  also  appears  as  positive  discontinuities  in  the  pressure  profUe. 

On  the  tail,  the  object  wake  is  made  to  fall  back  into  the  area  previously  occupied  by  the  object. 
Turbulent  eddies  are  created  by  the  air  viscosity  interaction  with  the  object  surface.  These  eddies  trail  into 
the  wake  core  and  support  a  region  of  low  pressure.  As  the  energy  in  the  eddies  dissipates  in  the  wake, 
the  region  collapses,  showing  up  as  the  final  shock  in  the  pressure  profile.  This  is  called  the  tail  shock, 
where  the  low  pressure  changes  to  ambient  Other  eddies  trailing  off  the  object  may  show  up  after  the 
tail  shock  as  strong  variations  in  pressure.  In  some  cases,  the  objects  shape  could  generate  second  and 
third  generations  of  low  pressure  wake  cores,  resulting  in  what  appears  as  several  tail  shocks. 


The  magnitude,  duration,  and  shape  of  the  N-Wave  for  an  ogive  body  of  revolution,  small-caliber 
shaped  projectile,  traveling  at  or  faster  than  the  speed  of  sound,  were  calculated  (Whitham  1952). 
Expressed  for  large  distances,  in  general  terms,  the  shock  overpressure  is. 


Po 


(1) 


where 

p  =  overpressure, 

Pjj  =  ambient  pressure, 

y  =  ratio  of  specific  heats  of  air  (1.4), 

M  =  Mach  number, 

y  =  a  characteristic  curve  relationship  between  x  and  O, 

X  =  distance  along  the  axis  from  the  bow, 

O  =  offset  distance, 

F(y)  =  a  function  used  to  solve  the  differential  equation  determining  the  flow  past  the  body  of 
revolution. 
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This  equation  formed  the  basis  of  the  theory  for  acoustic  round  discrimination.  By  simplifying  the 
generalized  results  into  simple  equations  (Skochko  1966),  one  can  calculate  the  pressure  and  total  duration 
of  the  N-Wave  by: 


_P_ 

Po 

T 


(.53)(M2-1)^/* 

D 

03/4 

(1.82)0 


D 


a(M2-i)3/8 


(2) 


where 

a  =  ambient  air  speed  of  sound  (meters/second), 
D  =  projectile  diameter  (meters), 

L  =  projectile  length  (meters), 

T  =  total  duration  (seconds). 


Theoretically,  a  properly  emplaced  sensitive  pressure  gauge  can  measure  the  N-Wave  of  passing 
projectiles.  Additionally,  if  the  ambient  air  conditions,  the  projectile  location,  and  velocity  relative  to  the 
pressure  gauge,  were  known,  the  projectile  could  be  identified.  A  system,  constructed  to  gather  data  and 
perform  analysis  for  projectile  identification. 


3.2  Acoustic  Location  and  Velocity  Technique. 

3.2.1  Background.  Determination  of  the  projectile’s  location  and  final  velocity  was  performed  with 
a  proven  acoustical  technique  constructed  of  simple  pressure  sensitive  devices  that  were  linked  to  a  clock. 
This  technology  was  developed  by  the  ARL  Human  Research  and  Engineering  Directorate  (HRED)  for 
smaU-caliber  projectiles  (Kalb,  unpublished).  Three  microphones  were  arranged  in  a  line  in  front  of  a 
target  and  placed  perpendicular  to  the  projectiles  flight  path  (see  Figure  4).  It  allowed  the  projectile’s 
location  above  the  microphone  array  to  be  determined.  On  HRED  ranges,  the  projectile  and  range  are 
known  prior  to  firing,  enabling  the  projectile’s  velocity  to  be  derived  from  ballistic  firing  tables.  Three 
additional  microphones  were  added  to  the  HRED  setup  to  create  the  RDS  microphone  setup  (see  Figure  5). 
The  additional  microphones  were  placed  directly  in  front  of  the  original  three  microphones,  yielding  a 
second  location  for  redundancy  in  case  a  microphone  was  damaged  or  failed.  The  two  center  microphones 
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Figure  4.  HRED  microphone  array. 

pix)vided  the  projectile’s  velocity.  Knowing  the  spatial  position  of  each  sensor  and  measuring  the  time 
each  was  triggered  by  the  passing  projectile’s  supersonic  shock  wave,  one  can  closely  approximate  the 
location  and  velocity  at  the  target 

3.2.2  Theory  and  Fomiulation.  The  projectile’s  shock  wave  can  be  visualized  as  an  infinite  cone 
traveling  with  the  projectile.  The  ambient  speed  of  sound  and  the  projectile’s  velocity  result  in  the  classic 
shock  angle  expression: 


a 


sin 


\.^J 


=  sin 


-1 


(±^ 


.  M  =  _, 
a 


(3) 
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Coordinate  System 


Figure  5.  Microphone  array  and  target  setup. 


where 

a  =  shock  angle  (°) 

M  =  Mach  number, 

V  =  projectile  velocity  (meters/second), 
a  =  ambient  speed  of  sound  (meters/second). 

This  relationship  assumes  that  the  cone  form  is  maintained,  subsequent  to  its  time  of  arrival,  throughout 
the  distance  of  the  shock  wave  propagation. 

The  projectile  is  assumed  to  be  traveling  parallel  to  the  Y-axis  striking  the  target  at  coordinates 
(X,Y,Z).  Microphones  1,  2,  and  3  are  placed  in  the  front  right,  front  center,  and  front  left,  respectively. 
Microphones  4,  5,  and  6  are  placed  in  the  back  right,  back  center,  and  back  left.  These  location 
microphones  (LM)  are  referred  to  as  LM1-LM6.  The  microphone  centers  for  microphones  1-6  are  located 
at  (X+1,  Y,Z),  (X,Y,Z),  (X-1,Y,Z),  (X+1,Y+1,Z),  (X,  Y+1,Z),  and  (X-1,Y+1,Z),  respectively.  Figure  5 
shows  the  microphone  placement  relative  to  the  target  and  the  coordinate  system  used. 
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On  the  arrival  of  the  shock  cone  at  each  microphone,  a  total  of  five  time  differences  (Tj,T2,T3,T4,  and 
T5)  are  recorded.  Two  time  differences  are  needed  for  each  location  array  and  one  time  difference  for  Ae 
velocity  determination.  T^-Tg  are  defined  as  follows: 

Tj  =  Time  difference  between  LM2  and  LMl, 

T2  =  Time  difference  between  LM3  and  LM2, 

T3  =  Time  difference  between  LM5  and  LM2, 

T4  =  Time  difference  between  LM5  and  LM4, 

T5  =  Time  difference  between  LM6  and  LM5. 

The  time  differences  T^  and  T2  are  measured  from  a  moving  sound  source  and  can  be  related  to  a 
hypothetical  source  that  originates  from  a  point  in  a  vertical  plane  defined  by  LMl,  LM2,  and  LM3.  The 
transformed  time  difference  from  each  microphone  pair  defines  a  hyperbola  of  ambiguity,  which  is  the 
set  of  all  point  sources  which  produce  this  same  arrival  time  difference.  The  intersection  of  these 
hyperbolas  locates  the  point  source  through  which  the  projectile  traveled  (X,Z).  Figure  6  is  an  example 
of  the  hyperbolas  of  ambiguity. 


The  altered  difference  in  arrival  time  between  the  center  and  right  microphone  pair,  Dj,  can  be 
determined  from  the  following  equation: 

(Rc  -  '^1  (4) 

Di  =  - - -  -  -3^’ 
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where 


=  in-plane  distance  from  center  microphone  to  projectile  location  (meters), 

Rj^  =  in-plane  distance  from  right  microphone  to  projectile  location  (meters), 

and  the  difference  in  arrival  time  between  the  left  and  center  microphone  pair,  D2,  is: 

^  ^  (Rl  -  Rc)  ^  ^2 

^  a  cos(a)  ’ 


where 


(5) 


Rl  =  in-plane  distance  from  left  microphone  to  projectile  location  (meters). 

It  is  useful  to  note  that  a  positive  time  indicates  that  the  source  is  to  the  right  of  the  midpoint  of  a 
microphone  pair  and  that  a  negative  time  indicates  that  the  source  is  to  the  left. 

The  sound  source  location  is  related  to  the  microphone  positions  by  using  the  law  of  cosines.  On  the  right 
side  triangle: 


Rr2  =  Rc^  +  -  2RcAcose  =  (Rc  -aDi)2  =  Rc^  -  2RcaDi  +  a^D^^ ,  (6) 


where 

A  =  microphone  separation  distance  (1  m), 

0  =  angle  from  microphone  array  axis  to  Rc- 

Similarly  for  the  left  side  triangle: 

Rj^^  =  Rc^  A^  “  2RcAcos(7r —0)  —  (Rq  **"  2Rc^i^2 

where  cos(7t  -  0)  =  -cos  (0). 
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Subtracting  equation  6  from  7  yields: 


4RcAcose  =  2Rca(D2+Di)  +  a  .  (8) 

From  geometry,  the  offset  distance  of  the  projectile  from  the  center  microphone  is: 

Rc  =  +  .  (9) 


Solving  equations  8  and  9  for  Z  yields: 


Z  =  R^cosG  =  a(D2+Di) 


a(D2-Di)^ 

2 


V. 


2A 


J 


(10) 


Rearranging  equation  9  for  X  yields: 


X  =  ^Rc^-Z^  .  (11) 

The  projectile  location  for  microphones  4,  5,  and  6  can  be  solved  using  the  same  equations  with  the 
substitution  of  and  D5. 

Knowing  that  the  position  of  the  microphones  are  fixed  and  the  projectile  is  assumed  to  be 
perpendicular  to  the  array,  one  can  approximate  the  velocity.  The  velocity  is  determined  by  taking  the 
difference  in  position  of  the  center  microphones  and  dividing  by  the  difference  in  arrival  times  of  the 
center  microphones. 


(12) 
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where 


Vp  =  predicted  projectile  velocity  (meters/second), 

Ylm2  =  y  position  of  LM2  (meters), 

^LM5  =  y  position  of  LM5  (meters). 

3.2.3  Alternate  Acoustic  Location  and  Velocity  Technique.  Another  location  and  velocity  prediction 
algorithm  based  on  microphone  detection  of  supersonic  shock  waves  was  developed  and  used  on  Yuma 
Proving  Ground  (YPG)  test  ranges  (Stallings  1992).  Stallings’  technique  does  not  require  any  special 
microphone  placement.  The  only  requirement  to  locate  projectiles  and  determine  velocity  is  to  have  at 
least  4  noncoplanar  microphones.  Three  of  the  microphones  form  a  triangle  in  the  target  plane,  and  the 
other  microphone  is  placed  in  front  of  the  others.  Generally,  the  more  microphones  used,  the  better  the 
location  accuracy  that  can  be  provided.  Better  location  accuracy  is  also  accomplished  when  the 
microphones  lie  close  to  the  target  plane  and  the  hits  lie  in  or  near  the  center  of  the  microphone  array. 
The  iiq>ut  variables  for  accurate  round  location  are  meteorological  data  (to  include  temperature,  humidity, 
wind  speed,  and  wind  direction),  a  velocity  estimate,  fall  angle  estimate,  target  position,  gun  position, 
microphone  positions,  and  microphone  trigger  times.  Location  errors  magnify  when  input  values  are  not 
precisely  known.  Stallings  showed,  through  YPG  tests,  that  errors  as  small  as  2-3  cm  are  common. 

Stallings’  technique  first  assumes  that  the  projectile  hits  the  target  at  the  center  of  mass.  The 
theoretical  times  that  the  shock  wave  was  detected  by  each  microphone  are  then  calculated.  The  actual 
microphone  times  are  compared  to  the  calculated  times.  The  output  variables,  projectile  location,  and 
velocity  are  incrementally  adjusted  around  the  initial  guess  until  the  difference  between  calculated  and 
actual  times  is  within  a  predetermined  tolerance.  This  is  a  time-intensive  algorithm  that  convolves  to  the 
solution  in  approximately  30  s  on  a  fast  personal  computer.  The  algorithm  could  be  optimized  to  reduce 
inefficiencies  in  the  program.  There  are  two  iterative  nested  processes  within  the  algorithm.  Stallings 
believes  that  one  of  these  can  be  eliminated  to  accelerate  the  program. 

The  Stallings  algorithm  cannot  predict  location  in  real  time.  Real  time  results  are  essential  for  proper 
integration  with  RETS.  Since  this  technique  was  developed  for  test  ranges,  many  of  the  input  variables 
are  known  to  a  large  accuracy  at  the  testing  range.  There  is  no  knowledge  of  the  gun  firing  position  on 
the  training  range.  For  the  purpose  of  RDS,  the  gun  position  represents  the  greatest  source  of  error  for 
this  technique  and  would  significantly  impact  scoring  accuracy. 
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4.  LIVE-FIRE  TESTING 


A  series  of  live-fire  tests  was  designed  to  investigate  the  hypothesis  that  acoustic  signatures  could  be 
used  to  discriminate  round  types.  The  smallest  and  largest  projectiles  within  the  test  matrix  were  selected 
for  the  initial  test  of  the  feasibility  of  the  hypothesis.  The  logic  was  that  if  the  data  obtained  during  the 
initial  testing  showed  that  the  characteristics  of  the  respective  N- Waves  permitted  discrimination  between 
the  projectile  extremes,  then,  and  only  then,  would  a  second  live-fire  test  be  conducted.  The  second  live- 
fire  test  would  be  a  data  collection  effort  to  compile  the  N-Wave  data  of  the  remaining  projectiles  at 
different  intervals  of  velocity.  Once  the  data  were  collected  and  analyzed,  an  algorithm  would  be 
developed  to  produce  round  discrimination.  Subsequenfiy,  the  algorithm  would  be  integrated  with  a 
computerA)readboard.  A  third  live-fire  test  would  be  conducted  to  validate  the  algorithm  and  the  RDS 
breadboard’s  performance.  Needless  to  say,  the  success  of  this  last  test  series  was  highly  dependent  on 
the  results  obtained  within  the  feasibility  and  data  collection  tests. 

These  tests  were  formally  described  as  the  feasibility  test,  data  collection  test,  and 
verification/validation  test  The  breakdown  of  rounds  and  ranges  used  within  each  test  series  is  given  in 
Table  2. 

The  data  collected  during  the  first  two  parts  of  the  live-fire  program  were  used  to  standardize  the 
method  of  collecting  the  N-Wave  data  and  the  specific  parameters  used  in  the  N-Wave  calculation.  To 
increase  the  amount  of  data  acquired,  several  sensors  were  used  for  each  shot.  The  parameters  considered 
for  the  N-Wave  calculations  were:  ambient  air  temperature,  projectile  velocity,  caliber,  length,  and  offset 
distance.  Each  live-fire  test  setup  consisted  of  a  fixed  reference  target  and  five  high-fidelity  microphones 
arranged  such  that  the  projected  distance  form  the  target  center  was  1,  3,  6,  9,  and  12  m,  respectively,  to 
the  microphones.  Each  projectile  hit  on  the  target  was  scored  using  a  video  scoring  system,  and  the  offset 
distance  from  the  projectile  to  the  microphones  was  calculated,  generating  one  parameter  for  consideration. 
The  projectile’s  velocity  was  measured  with  a  Weibel  radar  system,  designed  specifically  for  projectile 
velocity  measurement.  The  dimensions  of  each  round  were  known,  and  meteorological  data  were 
measured. 

To  ensure  that  there  was  no  unintentional  biasing  of  the  data,  the  order  of  fire  was  randomized  to  the 
extent  practical.  Theoretically,  if  an  experiment  involving  different  round  types  is  conducted,  the  test 
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Table  2.  Proposed  Round  and  Firing  Range  Breakdown  for  Each  Test 


Test  Type 

Round  Type 

Range 

(m) 

120-mm  M865 

400 

2,500 

120-mm  M831 

400 

2,500 

Feasibility 

105-mm  M724 

400 

2,500 

105-mm  M490 

400 

2,500 

25-mm  M910 

400 

2,000 

25-mm  M793 

400 

2,000 

120-mm  M865 

1,250 

120-mm  M831 

1,250 

105-mm  M724 

1,250 

105-mm  M490 

1,250 

Data  Collection 

35-mm  DM18 
35-mm  DM68 

400 

400 

1,250 

1,250 

25-mm  M910 

1,250 

2,500 

25-mm  M793 

1,250 

2,500 

12.7-mm  M33 

400 

900 

7.62-mm  M59 

400 

900 

120-mm  M865 
120-mm  M831 
105-mm  M724 
105-mm  M490 

VerificationA^  alidation 

35-mm  DM18 
35-mm  DM68 
25-mm  M910 
25-mm  M793 
12.7-mm  M33 
7.62-mm  M59 

TBD® 

®  To  be  determined. 


conditions  should  be  identical  for  all  round  types  so  that  the  round  differences  can  be  detected.  Since 
reality  dictates  that  the  process  of  extracting  data  is  never  the  same  each  time,  randomization  reduces  the 
average  effect  of  any  irrelevant  factors,  and  the  errors  in  the  response  are  independently  and  randomly 
distributed  over  the  population.  Conversely,  if  rounds  are  systematically  tested,  i.e.,  all  of  a  particular 
roimd  are  tested  at  one  time,  then  the  true  response  of  that  round  cannot  be  separated  from  the  extraneous 
factors  such  as  weather  conditions  that  contributed  to  the  response. 
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4.1  Feasibility  Test.  The  first  part  of  the  live-fire  program  was  intended  to  be  a  smaU  effort. 
Originally,  the  120-mm  HEAT  and  KE  training  projectiles  were  scheduled  to  be  fired  at  toe  closest  and 
farthest  engagement  distances,  and  toe  data  were  compared  to  toe  25-mm  armor-piercing  (AP)  and  HE 
training  projectiles  at  their  extreme  engagement  distances.  The  extremes  of  toe  test  matrix  were  designed 
to  envelop  toe  velocity  expected  in  real  use.  According  to  toe  simplified  N-Wave  calculations,  toe  N- 
Wave  pressure  and  time  duration  produced  fiom  toe  largest  projectile  at  its  slowest  velocity,  far  from  toe 
microphone,  was  expected  to  approach  toe  N-Wave  pressure  and  time  duration  of  toe  smallest  projectile 
at  its  fastest  velocity,  close  to  toe  microphone.  If  toe  round  type  of  these  two  projectiles  could  be 
discriminated,  round  discrimination  using  N-Wave  data  would  be  considered  feasible. 


4.1.1  Test  Purpose.  The  purpose  of  this  test  was  to  determine  if  toe  N-Wave  data  collected  with  a 
particular  configuration  of  high-fidelity  microphones  could  be  used  to  discriminate  toe  round  types  fired. 
The  largest  projectile,  toe  120-mm  M831,  at  its  slowest  velocity  was  to  be  compared  to  toe  25-mm  M910 
at  its  fastest  velocity.  If  toe  N-Wave  data  from  toe  same  offset  distance  were  separable,  then  these 
N- Waves  of  supersonic  projectiles  could  be  used  to  perform  round  discrimination  at  toe  target 

Due  to  the  extremely  high  operational  cost  of  Ml/MlAl  tanks,  and  the  prospect  of  incurring 
prohibitively  high  maintenance  costs  at  toe  turn  of  toe  fiscal  year,  it  was  decided  to  advance  a  portion  of 
toe  data  collection  effort  into  toe  feasibility  study.  Consequently,  all  120-mm,  105-mm,  and  25-mm 
projectiles  were  fired  during  toe  feasibility  test.  The  120-mm  rounds  were  fired  from  an  MlAl  Abrams 
tank,  toe  105-mm  projectiles  were  fired  from  a  M60A3  tank,  and  aU  toe  25-mm  rounds  were  fired  from 
an  M2  BFV.  Table  3  describes  toe  actual  rounds  fired  during  toe  feasibility  test.  With  toe  modification 
of  toe  test  schedule,  toe  sample  set  increased  by  50%  and  would  provide  better  parameter  estimates  in  toe 
statistical  analysis  for  this  phase  of  toe  test  program. 


Table  3.  Actual  Round  and  Firing  Range  Breakdown 
for  toe  Feasibility  Test 


Round  Type 

Range 

(m) 

120-mm  M865 

400 

1,000 

120-mm  M831 

400 

1,000 

HPH 

25-mm  M910 

400 

1,000 

HpH 

25-mm  M793 

400 

1,000 

105-mm  M724 

400 

1,000 

HPH 

105-mm  M490 

400 

1,000 

■EH 
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4.1.2  N-Wave  Characteristics.  The  acoustic  signature  for  each  projectile  was  studied  specifically  for 
seven  different  characteristics.  These  were  positive  peak  pressure,  negative  peak  pressure,  positive  phase 
duration,  negative  phase  duration,  total  duration,  peak  net  pressure  impulse,  and  peak  absolute  pressure 
impulse.  The  variables  are  illustrated  in  Figure  7.  The  data  were  compared  to  N-Wave  formulas  and 
analyzed  using  advanced  statistical  methods.  Shadowgraphs  of  the  projectiles  in  midflight  from  previous 
tests  were  compared  to  the  N-Wave  data  to  relate  the  signature  features  to  the  surface  features  of  the 
projectiles.  This  comparison  helped  determine  if  certain  N-Wave  features  were  truly  dependent  on  the 
projectile  geometry  or  were  the  results  of  an  environmental  influence,  such  as  a  ground  reflection. 

The  seven  characteristics  of  the  acoustic  signature  were  essential  to  the  N-Wave  data  comparisons. 
The  N-Wave  formulae  were  used  as  a  theoretical  guide  to  help  test  the  hypothesis  of  round  discrimination. 
The  real  N-Wave  data  were  used  to  empirically  refine  the  N-Wave  formulae.  These  formulae  provided 
an  engineering  basis  for  the  round  discrimination  algorithm  used  on  the  breadboard. 

4.1.3  Data  Acquisition  and  Reduction.  The  data  acquisition  system  consisted  of  two  systems  for 
recording  and  one  system  for  sensing.  The  sensors  were  high-fidelity  microphones  able  to  operate  in  the 
range  of  overpressures  encountered.  These  sensors  were  high-value  items  and  were  not  intended  as  an 
all-weather  solution  to  a  fielded  RDS.  The  priority  lay  in  obtaining  high-fidelity  signals  to  accurately 
record  the  N-Wave  phenomenon.  The  N-Wave  data  from  several  projectiles,  at  several  offset  distances, 
traveling  at  various  velocities,  had  to  be  studied  and  analyzed.  The  two  recording  systems  served  a  dual 
purpose.  Each  system  had  a  specific  deficiency — ^poor  resolution  in  either  time  or  voltage  levels.  Each 
was  able  to  complement  the  other  in  that  respect.  The  secondary  purpose  was  to  provide  a  backup 
recording  system  in  the  event  of  a  malfunction  during  the  live-fire  test. 

Reduction  of  the  data  was  performed  on  a  personal  computer  (PC).  The  data  from  each  chaimel  were 
converted  into  ASCII  files.  These  files  were  read  into  a  data  processing  program  called  DADisp  (DSP 
Development  Corporation  1991).  This  software  had  the  capability  to  read  the  data,  convert  them  into 
engineering  units,  and  perform  a  series  of  calculations.  Among  those  were  the  calculation  of  net  pressure 
impulse  and  absolute  pressure  impulse.  Along  with  these  calculations  and  plots,  DADisp  helped  to 
determine  the  pressure  peaks  and  durations.  At  times,  the  signals  were  clipped.  In  an  effort  to  maximize 
the  resolution  in  the  ordinate,  the  recording  window  was  set  as  small  as  possible.  At  times,  the  recording 
devices  were  preset  on  low.  As  a  result,  hard  data  on  the  maximum  levels  were  lost.  DADisp  software 
helped  curve-fit  the  peaks,  restoring  the  data  with  at  least  ±10%  accuracy. 
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4.1.4  Instrumentation.  B&K  microphones  are  high-grade  laboratory  instruments  with  die  following 
specifications:  sensitivity  of  -56-dB  sound  pressure  level  (SPL),  flat  frequency  response  ftom  less  than 
5  Hz  to  50  kHz,  and  dynamic  range  up  to  178  dB  SPL  with  3%  distortion.  Each  B&K  microphone  unit 
is  comprised  of  a  type  4136  1/4-in  microphone,  a  type  2639  preamplifier,  and  a  type  2801  AC  or  a  type 
2804  DC  power  supply.  The  cost  per  unit  is  $3.5k-$4k  in  small  quantities. 

The  outputs  of  the  microphones  were  connected  by  150-ft  BNC  cables  to  three  TEAC  digital  audio 
tape  (DAT)  recorders  and  a  Hewlett-Packard  (HP)  digital  oscilloscope,  model  no.  HP  541 12D,  controlled 
by  a  PC  located  behind  a  berai.  The  three  DAT  recorders  had  a  total  of  six  wide-band  channels  (two 
channels  per  DAT)  with  48-kHz  sampling  rate  and  14-bit  resolution  per  channel.  The  HP  digital 
oscilloscope  had  four  input  channels  with  1-MHz  sampling  rate  per  channel  and  8-bit  resolution  per 
channel.  The  PC  communicated  with  the  oscilloscope  through  the  IEEE  port.  A  data  acquisition  program 
written  in  C  language  was  developed  to  automatically  dump  64  kB  of  data  per  channel  onto  the  PC’s  hard 
disk  when  the  scope  was  triggered  and  to  reset  the  scope  to  wait  for  the  next  trigger.  Both  types  of 
recording  instruments  were  used  to  provide  both  time  resolution  and  amplitude  resolution. 

4.1.5  Acoustic  Signatures.  Live-fire  tests  were  conducted  at  the  H-field  facility  at  APG,  Edgewood 
Area  (EA).  Data  collected  during  the  feasibility  test  are  summarized  in  Appendix  A.  The  target  was 
positioned  up  on  top  of  a  berm  at  the  end  of  the  field.  The  firing  platform  was  moved  to  each  firing 
statioa  After  some  calibration  roimds  were  fired,  the  testing  for  each  projectile  type  was  performed.  A 
cursory  look  at  the  data,  while  still  on  the  field,  indicated  that  the  signals  for  each  round  type  was  distinct. 
Qualitatively,  it  was  surmised  that  round  discrimination  was  feasible.  A  subsequent  closer  look,  within 
the  laboratory  confirmed  this  fact,  quantitatively,  and  is  discussed  in  section  5. 

The  M831  and  M490  rounds  possessed  distinctly  different  N-Wave  shapes.  The  M490  round 
produced  the  classic  N-Wave,  as  seen  in  Figure  8,  while  the  M831  round  had  several  interval  shocks  on 
the  negative  phase  of  the  N-Wave.  As  the  offset  distance  varied,  the  level  and  occurrence  of  the  shocks 
changed.  This  indicated  that  the  shock  was  produced  from  a  surface  feature  and  was  merging  with  the 
tail  shock.  The  time  relation  of  shock  travel  and  the  distance  to  the  ground  disputed  the  idea  that  it  was 
a  ground  reflection  (Ritzel  and  Gottlieb  1982).  These  observations  were  distinct  to  the  M831  round  and 
not  the  M865. 
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It  was  speculated  that  these  pressure  differences  are  due  to  shock  interaction  with  the  fins,  and  the 
rionogival  geometry  of  the  M831  round.  The  M831  round  has  the  stand-off  fuse  and  the  concave  taper 
to  the  body,  as  opposed  to  the  convex  taper  of  the  ogive.  The  body  slims  down  aft  of  the  body  for 
aerodynamic  flight.  The  bottle  nose  feature  may  create  a  separation,  causing  a  shock,  unless  it  is 
isentropically  ramped.  If  the  surface  aft  of  the  body  turning  away  from  the  flow  is  of  sufficient  gradient, 
another  separation  would  occur,  creating  more  shocks.  A  detailed  summary  of  the  resultant  data  is 
included  in  Appendix  B.  figure  9  shows  the  same  comparison  between  the  M793  and  M910  fired  fi’om 
a  BFV.  The  slimmer  AP  rounds  shows  the  classic  N-Wave  while  the  HE  roimd  has  the  secondary  tail 
shock. 


Figure  9.  M793  and  M910  signature  comparisons. 

4.1.6  Round  Location  and  Velocity  Determination  Data.  The  RDS  acoustic  location  system  consisted 
of  six  electret  microphones.  Each  microphone  was  connected  so  that  it  acted  as  a  variable  resistor  when 
exposed  to  sound.  In  response  to  the  acoustic  shock  wave,  the  current  in  the  microphone  wires  increased. 
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The  current  was  converted  to  a  voltage  signal  by  a  signal  conditioning  card.  The  six  pulses  were  recorded 
on  two  DATs.  LMl,  LM2,  and  LM3  were  recorded  on  DAT  1,  and  LM4,  LM5,  LM6,  and  LM2  were 
recorded  on  DAT  2.  LM2  was  recorded  on  both  DATs  so  that  all  six  microphone  trigger  times  could  be 
related  to  one  time  base.  The  DAT  clock  was  synchronized  to  the  Weibel  radar  clock.  Each  channel  was 
played  back  through  a  Nicolet  digital  oscilloscope  and  manually  reduced  for  the  time  passed  until  shock 
detection  (see  Appendix  C).  The  accuracy  of  these  times  was  within  ±21  ps.  The  time  difference  pulses 
Tj-Tj  were  then  calculated. 

4.2  Data  Collection  Test.  The  result  of  the  feasibility  test  suggested  that  round  discrimination  using 
acoustic  signatures  was  possible  and  the  data  collection  test  should  proceed.  The  ammumtion  used  in  the 
data  collection  test  is  shown  in  Table  4.  This  portion  of  the  project  was  originally  intended  to  be  the 
major  data  collection  effort.  However,  since  all  of  the  main  gun  data  needed  for  the  M60A3,  Ml/Al,  and 
M2  were  collected  during  the  feasibility  test,  all  that  remained  were  the  TPGID  and  small  arms. 


Table  4.  Actual  Breakdown  of  Rounds  and  Ranges 
for  Data  Collection  Test 


Round  Type 

Range 

(m) 

35-mm  DM18 

400 

1,000 

1,700 

35-mm  DM68 

400 

1,000 

1,700 

35-mm  TPL-T 

400 

1,000 

1,400 

12.7-mm  M33 

400 

1,000 

12.7-mm  M17 

400 

1,000 

7.62-mm  M59 

400 

600 

800 

7.62-mm  M62 

400 

The  small  arms  ammunition  used  on  a  training  range  is  either  ball  or  tracer.  The  tracer  ammumtion 
was  studied  as  an  independent  round  type  since  the  effect  of  the  tracer  was  postulated  to  have  an  effect 
on  the  N-Wave.  It  was  theorized  that  the  pyrotechnic  material  would  increase  the  duration  since  it  would 
add  energy  to  the  wake.  The  degree  of  effect  was  pure  conjecture.  Only  one  type  of  35-mm  TPGID  was 
supposed  to  be  tested.  Before  the  data  collection  test  occurred,  two  new  types  of  TPGID  rounds,  DM68 
and  the  TPL-T,  were  obtained  for  testing.  These  were  surplus  obtained  from  a  lot  which  had  just 
undergone  proof-testing  with  the  U.S.  Army’s  Combat  Systems  Test  Activity  (CSTA).  These  rounds  were 
tested  independently  of  the  DM18.  Their  diameter  and  length  were  unchanged,  but  the  surface  features 
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were  dilferent.  The  DM68  and  the  TPL-T  have  a  sharp  pointed  cone  attached  to  a  cylinder  body,  as 
opposed  to  the  DM18,  which  has  an  ogive  shape.  It  was  speculated  that  there  would  be  a  slight  change 
in  the  N-Wave  but  that  the  change  would  not  be  significant 

The  purpose  of  the  data  collection  test  was  to  complete  the  database  of  projectile  signatures  expected 
on  a  U.S.  Army  training  range.  These  signatures  were  needed  to  develop  the  RDS  algorithm  and  form 
the  basis  of  aU  assumptions  and  work.  Once  the  data  were  collected  and  reduced,  a  major  milestone  in 
the  project  would  have  been  reached. 

To  perform  the  test,  an  MlAl  Abrams  tank  was  needed  to  fire  the  TPGID  arrununition.  A  Marm 
barrel,  a  device  constracted  for  precise  single-shot  testing  of  small  arms  ammunition,  was  used  to  perform 
the  7.62-mm  and  12.7-mm  shots.  The  Marm  barrel  consisted  of  a  towed  gun  carriage  and  two  barrel 
assemblies,  one  for  each  caliber. 

The  effort  to  field  a  calibrated  pressure  sensor  was  deemed  impractical  after  considering  the  number 
of  sensors  needed  to  cover  all  the  U.S.  Army  live-fire  gurmery  ranges.  This  made  any  parameter  not 
dependent  on  the  pressure  levels  very  attractive.  The  total  duration  was  one  parameter  that  seemed  to  be 
most  consistent,  and  discriminating.  The  same  seven  variables  examined  in  the  feasibility  test  were 
studied  in  the  data  collection  effort,  in  the  likelihood  that  the  total  duration  failed  as  a  discriminator  with 
inclusion  of  the  additional  data. 

The  same  two  recording  systems  used  during  the  feasibility  test  were  employed  for  the  data  collection 
test.  A  third  recorder  was  added.  This  was  a  digital  transient  data  recorder  card  mounted  in  a  PC.  It  had 
12-bit  resolution  in  the  ordinate,  with  a  250-kHz  sampling  rate.  This  was  deemed  sufficient  response  for 
the  7.62-mm  rounds.  The  resultant  plots  were  much  smoother  in  appearance,  and  the  accuracy  of  the  data 
was  enhanced. 

Five  B&K  high-fidelity  microphones  were  set  up  at  1  m,  3  m,  6  m,  9  m,  and  12  m,  respectively,  to 
the  left  of  the  bull’s-eye  in  fiont  of  the  target.  In  addition,  three  low-cost  ($150.00  each)  high-fidelity 
microphones  from  Countryman  were  set  up  next  to  the  B&K  microphones  at  3  m,  6  m,  and  9  m.  The 
Countryman  microphones  were  much  more  raggedized  and  rain  resistant  than  the  B&Ks. 
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A  new  data  collection  system  was  used.  The  system  consisted  of  a  486  Extended  Industry  Standard 
Architecture  (EISA)  computer,  two  four-channel  1-MHz  EISA  data  acquisition  cards,  and  a  graphical 
programming  software  tool  called  LabView  for  Windows.  The  system  emulated  an  eight-channel  high¬ 
speed  (250  kHz  per  channel)  digital  oscilloscope  with  the  programming  capabilities  of  a  PC.  The  acoustic 
signatures  of  the  microphones  were  automatically  captured  and  displayed  on  the  monitor,  scaled,  truncated, 
and  stored  on  the  hard  disk  after  each  tank  round  was  fired.  This  new  system  greatly  reduced  the  data 
processing  and  reduction  time.  The  same  four-channel  1-MHz  HP  digital  oscilloscope  was  used  as  a 
backup  system. 

Post-test  data  analysis  showed  that  the  Countryman  microphor^s  did  not  have  the  65-dB  dynamic 
range  necessary  to  capture  the  supersonic  shock  wave  of  large  projectiles.  For  small  rounds  such  as 
12.7  mm  and  the  7.62  mm,  the  Countryman  microphones  responded  wefi. 

A  second  series  of  live-fire  tests  was  also  conducted  at  H-Field  facility  at  the  APG,  EA.  Data  from 
these  tests  are  compiled  in  Appendix  D.  The  same  equipment  used  during  the  feasibility  test  was 
employed  during  the  data  collection  test.  The  target  was  positioned  on  top  of  the  same  berm  used  during 
the  feasibility  tests.  The  firing  platforms  were  moved  to  specific  firing  stations.  After  some  calibration 
rounds  were  fired,  the  testing  for  each  projectile  type  was  performed.  An  immediate  look  at  the  data  on 
the  field  indicated  that  the  signals  for  the  small-caliber  rounds  were  different  from  the  main  gun  roimds. 
This  observation  was  subsequently  verified  during  the  detailed  analysis. 

The  bulk  of  the  data  considered  was  from  the  transient  data  recorders.  The  resolution  in  the  abscissa 
and  ordinate  made  the  data  optimal  when  performing  data  reduction.  Figure  10  shows  the  typical  plot  of 
a  35-mm  DM18  projectile  raw  data,  fired  from  1,700  m  at  the  9-m  offset  distance.*  See  Appendix  E  for 
details.  Appendix  F  summarizes  the  microphone  positions  and  associated  shock  detection  times. 


*The  N-wave  appears  upside  down  because  the  signal  from  the  microphone  was  increased  by  an  inverting  amplifier. 
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Pressure  Signal 
(V) 


Figure  10.  Raw  data  of  35-mm  DM18  at  9-m  offset,  fired  from  1.700  m. 


5.  ANALYSIS  OF  LIVE-FIRE  DATA 


The  live-fire  data  was  analyzed  using  two  distinct  approaches.  The  first  was  an  engineering  approach, 
based  on  Whithams  formula,  that  described  a  mathematical  model  which  related  the  physical  properties 
and  events.  This  model  was  thought  to  be  adequate  for  predicting  occurrences  within  and  beyond  the  test 
matrix.  The  second  approach  used  advanced  statistical  analysis.  By  analyzing  the  data  on  the  variables, 
a  empirical  model  was  constructed  to  predict  within  the  test  matrix. 

5.1  Variables  Investigated.  The  seven  variables  investigated  were  as  follows:  positive  peak  pressure, 
negative  peak  pressure,  positive  phase  duration,  negative  phase  duration,  total  duration,  peak  net  pressure 
impulse,  and  peak  absolute  pressure  impulse.  It  was  hypothesized  that  the  positive  peak  pressure  and 
negative  peak  pressure  generated  from  the  main  gun  rounds  would  be  so  much  stronger  than  the  BFV 
roimds  that  merely  sensing  the  levels  would  provide  discrimination.  Subsequent  calculations  negated  this 
hypothesis.  A  small  round  close  to  the  sensor  would  appear  stronger  than  a  large  round  far  from  the 
sensor.  The  drop  in  peak  pressure  with  distance  was  estimated  to  be  roughly  to  the  power  of  l/5th. 
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The  predictions  of  time  durations  were  greatly  varied  over  round  type.  According  to  the  time  duration 
formula,  the  total  duration  is  strongly  related  to  the  diameter  of  the  projectile.  Since  the  rounds  all  have 
distinct  diameters,  and  the  other  variables  of  the  time  duration  formula  can  be  quantified,  time  duration 
should  be  a  pronoimced  discriminator.  Time  duration  does  not  require  any  calibration.  All  that  is  needed 
is  an  accurate,  high-speed  clock.  The  intervals  of  positive  and  negative  pressure  phase  time  durations  can 
also  be  used  to  identify  the  originator  of  an  N-Wave. 

There  was  much  speculation  that  the  environment  would  significantly  influence  the  variables 
mentioned.  The  temperature  of  the  ambient  air,  humidity,  air  density,  wind,  and  dust  were  thought  to  be 
capable  of  distorting  the  perfect  N-Wave  such  that  the  N-Wave  of  a  small  round  fired  on  a  cold  day  at 
sea  level  might  appear  the  same  as  a  large  projectile  N-Wave  fired  at  high  altitude  on  a  hot  day.  To 
attempt  to  evaluate  the  N-Wave,  a  measure  of  the  energy  imparted  to  the  atmosphere  was  used.  The  total 
duration  pressure  impulse  is  simply  the  integral  of  the  entire  N-Wave  pressure  curve.  The  value  of  the 
pressure  impulse!  at  some  point  in  time  along  the  N-Wave  is  the  net  value  of  energy  contributed  to  the 
local  environment.  As  the  projectile  flies  through  space,  the  KE  is  transferred  by  the  viscous  effects  of 
the  surrounding  fluid.  The  majority  of  the  energy  transfer  is  performed  in  the  nearly  isentropic 
compression  of  the  local  air,  thus  the  bow  and  tail  shock.  A  small  percentage  of  the  energy  is  consumed 
in  viscous  losses,  such  as  the  trailing  wake  eddies,  flow  separation,  friction  heating  of  the  projectile,  and 
the  losses  incurred  by  a  not  perfectly  isentropic  compression  of  the  air. 

The  total  duration  pressure  impulse  shows  the  entire  amount  of  net  energy  imparted  to  the  surrounding 
fluid.  The  negative  value  at  the  end  of  the  total  duration  pressure  impulse  curve  is  due  to  energy  loss  to 
the  surrounding  fluid.  The  projectile  passes  throu^  a  point  in  space  and  momentarily  displaces  the  air. 
The  energy  transferred  to  the  nearly  isentropic  bow  shock  radiates  from  the  bow.  The  air  moves  a  little, 
but  the  energy  in  the  shock  travels  at  the  local  speed  of  sound.  As  the  projectile  passes,  the  local  energy 
state  of  the  air  returns  to  the  original  level,  plus  any  energy  imparted  to  it  by  projectile  kinetic  losses.  The 
air  shocks  back  into  the  wake  as  the  wake  eddies  dissipate  energy  and  collapse.  This  marks  the  end  of 
the  energy  transference.  The  net  effect  is  a  negative  value.  The  value  should  be  approximately  equal  to 
the  losses  imparted  to  the  air,  which  result  in  local  heating  of  the  air. 

The  peak  value  of  the  total  duration  pressure  impulse  curve  demonstrates  the  total  duration  energy 
transferred  firom  the  travelling  projectile  to  the  local  air,  at  that  instant  in  time.  The  energy  transfer  of  the 
projectile  is  strictly  a  function  of  speed  and  projectile  shape.  Ambient  air  conditions  have  a  negligible 
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effect  on  this  value.  As  a  result,  a  projectile  could  be  discriminated  by  this  feature  without  regard  to  the 
environmental  conditions.  Another  method  of  using  this  variable  is  to  consider  the  absolute  contribution 
of  energy  transference. 

The  absolute  pressure  impulse  is  taken  when  considering  all  the  values  of  the  N-Wave  as  positive 
only.  In  this  manner,  the  absolute  total  duration  of  energy  transfer  is  considered.  This  works  such  that 
any  losses  in  the  projectile’s  flight  are  included  into  this  value  and  can  further  separate  the  groups  of 
projectiles  for  discrimination. 

Three  different  engineering  methods  of  analyzing  the  data  were  pursued.  The  first  method  was  the 
use  of  the  seven  variables  in  conjunction  with  the  time  duration  and  the  overpressure  formulations.  This 
method  would  provide  the  discriminant  algorithm  with  a  theoretically  derived  mathematical  model  that 
employed  thermodynamic  relationships.  This  model  could  be  used  to  discriminate  rounds  if  the  live-fire 
data  matched  the  calculated  results.  If  not,  a  modification  of  this  formulation  could  be  performed, 
resulting  in  quasi-empirical  formulation. 

The  other  two  methods  were  similar  to  each  other  in  process.  A  Fourier  transform  of  the  data  could 
produce  key  frequency  features  that  could  be  used  to  discriminate  the  projectiles.  Another  method  of 
transfer  analysis  is  called  Wavelet  analysis.  Using  this  method,  the  signal  is  transferred  into  the  domain 
of  a  discrete  function,  rather  than  the  continuous  function  of  the  frequency  domain.  By  picking  a  discrete 
function,  the  signal  is  transferred  by  taking  the  convolution  integral.  This  method  works  if  there  is  a 
specific  shape  to  each  N-Wave  when  normalized  in  the  abscissa  and  ordinate. 

The  results  from  the  live-fire  data  collection  efforts  were  N-Waves  that  were  extremely  similar  in 
shape.  The  Wavelet  theory  would  not  work  since  the  normalized  shape  of  the  N-Waves  were  so  closely 
matched.  A  Fourier  analysis  was  initially  performed  (Hulet  1992),  showing  extremely  promising  results 
(see  Figures  11-14).  The  peak  frequency  for  the  12.7  mm  was  3.58  kHz;  the  7.62  mm  was  4.15  kHz. 
Unfortunately  the  data  had  undergone  some  smoothing  process,  and  the  resulting  Fourier  analysis  yielded 
artificial  suppressing  frequency  peaks.  When  the  ARL  live-fire  data  were  transformed,  the  results  did  not 
yield  adequate  discriminants.  The  Fourier  transform  was  shown  to  relate  only  to  the  duration  intervals. 

Since  the  durations  were  strongly  related  to  a  rounds  specific  dimensions  and  the  transformation 
techniques  were  clearly  influenced  by  the  durations,  the  time  duration  formula  was  selected  as  the  best 
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Figure  11.  7.62-mm  projectile  at  1.500  m. 
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possible,  practical  method  of  discrimination.  The  need  to  calibrate  sensitive  pressure  gauges  for  a  mass 
production  item  was  impractical.  Therefore,  aU  the  variables  that  required  pressure  levels,  such  as  positive 
and  negative  pressure,  peak  value  of  total  duration  pressure  impulse,  total  duration,  and  absolute  pressure 
impulse,  were  disregarded.  These  data  were  retained  for  use  in  the  event  that  the  time  durations  could 
not  be  used  to  perform  a  high  rate  of  round  discrimination. 

The  total  time  duration  was  compared  to  the  data  and  found  to  compare  well.  The  resulting  computed 
durations  were  constantly  proportional  to  the  data.  This  indicates  that  the  earlier  assumption,  that  the 
formula  could  be  applied  to  large-caliber,  non-ogival-sht^  projectiles  was  good.  The  fairly  constant 
proportionality  for  each  projectile  resulted  in  a  stochastic  modification  of  the  duration  formula  with  a 
shape  constant. 

5.2  Feasibility  Study. 

5.2.1  Live-Fire  Data.  An  empirical  model  was  developed  fiom  the  feasibility  test  data  to  determine 
the  feasibility  of  round  discrimination.  The  following  sections  describe  the  data  analyzed  and  statistical 
tests  implemented.  The  analysis  begins  with  examination  of  the  1 ,000-m  data  and  ends  with  three  models 
developed  using  a  statistical  tool  known  as  discriminant  analysis.  The  number  of  roimds  in  the  analysis 
from  each  round  type  is  listed  in  Table  5.  As  described  earlier,  for  each  round  fired,  five  signatures  were 
collected  from  microphones  placed  1,  3,  6,  9,  and  12  m,  respectively,  from  the  center  of  the  target 


Table  5.  Number  of  Rounds  Analyzed  in  Feasibility  Study 


Distance 

(m) 

25  mm 

105  mm 

120  mm 

M793 

M910 

M490 

M724 

M831 

M865 

400 

4 

3 

3 

2 

3 

3 

1,000 

3 

3 

3 

3 

3 

3 

1,400 

3 

3 

- 

- 

— 

— 

1,700 

- 

- 

3 

3 

3 

2 

The  7th  ATC’s  original  minimum  requirement  was  to  distinguish  between  rounds  <30  mm  fiom  those 
>30  mm.  A  more  ambitious  task  was  undertaken  to  try  to  distinguish  among  all  the  round  types.  This 
goal  produced  six  categories  for  discrimination:  25-mm  KE,  25-rnm  HEI,  105-mm  KE,  105-mm  HEAT, 
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120-mm  KE,  and  120-mm  HEAT.  The  data  collected  Irom  the  signatures,  location  study,  and  Weibel 
radar  were  all  candidates  for  inclusion  in  the  models  for  discrimination. 

5.2.2  1,000-m  Data.  Preliminary  analyses  to  determine  the  feasibility  of  discriminating  different 
round  types  were  initially  performed  on  the  1,000-m  data.  Statistical  tests  were  performed  on  each 
variable,  and  at  each  microphone  separately.  These  tests  sought  to  determine  if  one  of  the  variables  could 
be  used  to  discriminate  the  six  round  types.  These  tests  also  sought  to  detennine  the  optimal  microphone 
distance  for  future  implementation. 

The  1-m  microphone  often  produced  saturated  signals.  The  actual  offset  distance  fell  into  four 
regions.  Region  1,  data  collected  from  both  the  1-m  and  3-m  microphones,  had  offset  distances  less  than 
4  m.  Region  2,  data  collected  from  the  6-m  microphone,  had  offset  distances  from  4  m  to  7  m.  Region 
3  has  offset  distances  >7  m  and  up  to  10  m.  Region  4  offset  distances  were  >10  m. 

The  Kruskal-Wallis  test  was  used  to  test  the  hypothesis  that  the  six  round  types  have  equivalent  mean 
values  in  comparison  to  the  alternative,  which  is  that  at  least  two  round  types  differ.  This  nonparametric 
test  performed  an  analysis  on  the  ranks  of  the  data  rather  than  the  actual  data  values.  The  ranks  started 
with  1  for  the  smallest  value  among  the  six  roimd  types.  The  largest  value  from  all  six  round  types 
received  a  rank  of  N,  the  number  of  test  data  from  the  six  populations.  The  advantage  of  the 
nonparametric  test  was  that  the  test  did  not  require  assumptions  of  normality  of  the  data  or  homogeneity 
of  variances  (HOV) — ^that  is,  equality  of  variance  among  the  groups.  Appendix  M  provides  an  example 
of  implementing  the  Kruskal-Wallis  test  and  provides  detailed  results  of  the  test  for  each  variable. 

The  results  of  the  Kruskal-Wallis  test  performed  on  each  variable  indicate  that  at  least  two  round  types 
differed  from  each  other  in  each  of  the  four  regions.  To  determine  exactly  which  round  types  differed, 
the  corresponding  multiple  comparison  test  was  implemented.  In  summary,  most  of  the  variables 
distinguished  between  the  small  and  large-caliber  rounds — that  is,  <30  mm  vs.  >30  mm — ^satisfying  the 
7th  ATC’s  minimum  requirement.  The  only  variable  that  did  not  was  negative  phase  duration  in  regions  3 
and  4. 

All  the  other  variables  in  each  region  distinguished  between  large  HEAT  and  large  KE  rounds.  The 
105-mm  and  120-mm  training  HEAT  rounds  were  not  statistically  different  fiom  each  other,  nor  were  the 
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105-mm  and  120-mm  training  KE  rounds.  But  the  105-mm  and  120-mm  training  HEAT  rounds  were 
statistically  different  from  the  105-mm  and  120-mm  training  KE  rounds. 

For  some  variables,  the  two  smaller-caliber  rounds  were  distinguishable  from  each  other.  These  results 
are  summarized  with  the  following  notation: 

25  K  25  H  (105  K  =  120  K)  (105  H  =  120  H). 

The  round  types  are  listed  in  order  from  lowest  average  rank  to  the  highest.  Parentheses  around  two 
round  types  indicate  there  is  no  significant  difference  between  the  rounds.  What  the  nonparametric 
analyses  indicated  was  that  many  of  the  variables  under  consideration  were  good  candidates  for 
discriminating  at  least  the  small-caliber  rounds  from  the  large-caliber  rounds  for  1,000-m  data.  The 
duration  and  pressure  impulse  values  were  good  for  discriminating  the  rounds  in  the  four  categories 
grouped  previously  (i.e.,  small  KE,  small  HEAT,  large  KE,  and  large  HEAT).  The  first  region  (offset 
distance  <4  m),  representing  data  collected  close  to  the  flight  of  the  round,  was  the  hardest  region  in 
which  to  discriminate.  The  facts  provided  a  guide  to  the  placement  of  the  microphones  for  future  tests. 

5.2.3  Analysis  of  Data  From  All  Ranges.  The  preliminary  analysis  of  the  data  at  the  1,000-m  range 
indicated  that  several  of  the  variables  were  independently  able  to  separate  the  roimd  types.  The  goal  was 
to  find  one  or  more  variables  that  could  separate  the  round  types  and  was  insensitive  to  range.  As 
mentioned  in  the  previous  analyses,  it  was  difficult  in  some  cases  to  separate  the  M724  from  the  M865 
rounds  and  difficult  to  separate  the  M490  and  the  M831  rounds  from  each  other.  Therefore,  in  the  next 
analyses,  the  discrimination  was  among  4  groups  (not  6):  M910,  M793,  M724/M865,  and  M490/M831. 

•  Visiud  Inspections.  Various  plots  of  the  signature  variables  of  interest  vs.  microphone  distance  were 
examined  to  determine  if  there  were  obvious  trends  in  the  data  that  separated  the  round  types,  independent 
of  the  range.  Positive  peak  pressure  was  plotted  in  Figure  15.  Some  separation  in  the  data  between  the 
large-caliber  rounds  and  the  small-caliber  rounds  was  observed.  However,  among  the  large-caliber  rounds, 
separating  the  HEAT  from  the  KE  rounds  over  all  the  distances  is  not  easy.  This  indicated  that  positive 
peak  pressure  is  sensitive  to  range  and  would  not  be  a  good  candidate  for  discrimination. 

Upon  examiriing  the  plots  of  total  duration  vs.  microphone  offset  distance  for  each  firing  position 
(range),  little  change  in  the  shape  of  the  curves  was  observed.  Figure  16  shows  all  the  ranges  together 
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Figure  15.  Positive  peak  pressure  vs.  microphone  distance  for  aU  firing  positions. 


Figure  16.  Total  duration  vs.  microphone  distance  for  all  firing  positions. 
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on  one  graph.  This  plot  suggested  that  total  duration  is  not  sensitive  to  range  and  is  a  good  candidate  for 
discriminating  round  types.  If  microphone  offset  distance  is  known,  the  type  of  round  can  be  determined 
based  on  the  magnitude  of  the  total  duration.  In  some  cases,  the  roimd  can  be  determined  independent 
of  the  microphone  offset  distance.  For  instance,  if  the  total  duration  were  >1,500  ps,  the  round  would  be 
a  large  HEAT  round  (M831  or  M490),  independent  of  the  microphone  offset  distance.  However,  if  the 
total  duration  were  600  ps,  one  could  not  be  certain  if  it  were  a  large  HEAT  round  with  a  microphone 
offset  distance  <2  m,  or  a  small  HEAT  round,  M793,  with  a  microphone  offset  distance  of  approximately 
12  m.  These  observations  led  to  conclude  that  total  duration  alone  carmot  be  use  to  discriminate  among 
the  rounds,  and  a  more  objective  and  quantifiable  way  to  combine  variables  and  separate  the  groups  was 
needed. 

5.2.4  Discriminant  Analysis  for  All  Ranges. 

•  General.  Discriminant  analysis  is  a  statistical  technique  which  uses  linear  combinations  of  variables 
to  determine  if  one  or  more  variables  will  categorize  the  populations  (round  types),  and  it  is  described  by 
Marriott  (1990)  in  the  following  way: 

Given  a  set  of  multivariate  observations  on  samples,  known  with  certainty  to  come  from  two  or 
more  populations,  the  problem  is  to  set  up  some  rule  which  will  allocate  further  individuals  to  the 
correct  population  of  origin  with  minimal  probability  of  misclassification.  This  problem  and 
sundry  elaboration  of  it  give  rise  to  discriminatory  analysis. 

In  linear  regression,  one  develops  a  linear  model  with  one  or  more  explanatory  (independent)  variables 
Xj’s  to  describe  the  response  (dependent)  variable  y.  A  linear  equation  is  of  the  form: 

y  =  Bq+BjXj+B2X2"*'B2X2+B^X4+ . . 

where  the  coefficients  Bj’s  are  linear.  The  x^’s  can  take  a  nonlinear  form  such  as  x  raised  to  a  power  or 
a  functional  form  as  the  natural  logarithm  or  exponent. 

In  the  application  of  categorizing  (or  classifying)  the  round  into  the  four  groups  (M910,  M793, 
M724/M865,  and  M490/M831),  round  type  is  not  a  measurable  response;  therefore,  univariate  regression 
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analysis  is  not  appropriate.  The  multivariate  tool,  discriminant  analysis,  allows  us  to  form  one  or  more 
functions,  where  the  value  of  the  function(s),  in  combination,  categorizes  the  round  types. 

fi  =  ao  +  ajxi  +  82X2  +  33X3  +  84X4  +  .... 

f2  =  bo  +  bjXi  +  b2X2  +  b3X3  +  b4X4  +  .... 

f3  =  Co  +  CjXi  +  C2X2  +  C3X3  +  C4X4  +  .... 

f„  =  zo  +  ZjXj  +  Z2X2  +  Z3X3  +  Z4X4  +  .... 

Only  variables  that  are  not  highly  correlated  with  one  another  and  contribute  significantly  are  selected 
into  the  model.  Forming  the  fimctions  is  performed  in  a  stepwise  fashion,  with  the  most  "significant 
variable"  that  most  separates  (and  therefore  discriminates)  among  the  round  types  entering  the  model  first. 
Given  that  this  first  variable  is  in  the  model,  the  next  variable,  which  is  not  highly  correlated  with  the 
variable(s)  already  in  the  model  and  would  contribute  most  from  the  remaining  significant  variables,  is 
entered  into  tire  model.  This  process  continues  xmtil  no  variable  can  satisfy  the  two  aforementioned 
conditions. 

Several  criteria  are  available  for  entering  significant  variables  into  the  model  (SPSS  Inc.  1990a). 
WUk’s  Lambda  was  used  in  this  analysis.  Wilk’s  Lambda  is  the  ratio  of  the  within-groups  sum  of  squares 
to  the  total  sum  of  squares  when  the  discriminant  scores  (the  fj’s)  are  the  dependent  variables,  and  the 
group  is  the  independent  variable. 

The  tolerance  for  correlation  of  variables  was  set  at  0.1.  Tolerance  is  the  amount  of  linear  association 
between  the  incoming  variable  and  the  variable(s)  already  in  the  model.  It  is  equal  to  l-R^^,  where 
is  the  squared  multiple  correlation  coefficient.  Variables  with  small  tolerance  values  (and  therefore  large 
r.2  values)  are  not  entered  into  the  model. 

These  analyses  were  performed  using  a  statistical  software  package  called  SPSS.  With  a  set  of 
functions  to  characterize  the  discrimination,  a  territorial  map  of  the  categories  of  at  most  two  dimensions 
is  provided.  This  procedure  is  explained  in  more  detail,  with  actual  data  from  the  test  program,  in  the  next 
paragraph. 
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•  Discriminatory  Functions  for  Live  Fire  Data  Set.  Wilk’s  Lambda  criterion  was  implemented  for  the 
model  development  of  the  live  fire  data  set.  Total  duration  was  the  first  variable  to  enter  the  model  and 
is  the  one  variable  that  best  separates  the  round  types  into  their  appropriate  groups.  As  described  earlier, 
total  duration  alone  carmot  guarantee  correct  classification  of  all  the  rounds.  The  stepwise  process 
continued  with  the  following  variables: 

1.  Total  Duration 

2.  Peak  of  Absolute  Pressure  Impulse 

3.  Velocity 

4.  Microphonic  Off-Set  Distance 

5.  Negative  Peak  Pressure 

6.  Positive  Peak  Pressure 

7.  Subtotal  Duration  (addition  of  positive  and  negative  phases). 

These  7  variables  were  able  to  correctly  classify  99.50%  of  the  round  types.  As  shown  in  Table  6.  There 
was  one  misclassification:  M793  classified  as  M910. 


Table  6.  aassification  Result  for  Original  Seven  Variables 


No.  of 
Cases 

Predicted  Group  Members 

hip 

Actual  Group 

1 

M910 

2 

M793 

3 

M724/M865 

4 

M490/M831 

1.  M910 

36 

36 

100% 

2.  M793 

36 

i 

I  : 

35 

97.2% 

3.  M724/M865 

61 

61 

100% 

4.  M490/M83 

69 

69 

100% 

If  one  looks  back  at  Figure  16  (total  duration  vs.  microphone  distance),  some  curvature  in  the  data 
is  seen.  The  natural  logarithm  (In)  of  total  duration  will  linearize  the  plot  as  shown  in  Figure  17. 
Repeating  the  discriminant  analysis  process,  the  natural  logarithm  of  total  duration  is  the  first  variable  to 
enter  the  model.  Continuing  the  analysis,  the  model  contains  the  following  variables: 
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1.  ln(Total  Duration) 

2.  Offset  Distance 

3.  Velocity 

4.  Total  Duration 

5.  Peak  of  Absolute  Pressure  Impulse 

6.  Negative  Peak  Pressure 

7.  Peak  of  Pressure  Impulse. 


Figure  17.  Initotal  duration)  vs.  microphone  distance. 


These  variables  produced  all  correct  classification  results.  Note  that  both  the  natural  logarithm  and 
raw  units  of  total  duration  entered  into  the  model.  Although  the  two  are  not  independent  of  one  another, 
they  pass  the  tolerance  criterion.  However,  the  pooled  within-groups  correlation  between  the  two  is  0.88. 
To  avoid  problems  with  multicoUinearity,  an  imstable  model  as  a  result  of  highly  correlated  explanatory 
variables,  it  was  not  advisable  to  keep  both  of  these  variables  in  the  model. 

From  the  seven  variables  listed  previously,  subsets  of  variables  were  formed  in  search  of  smaller 
models  that  may  be  equally  successful.  Three  reduced  models  will  be  explained  in  the  following  sections. 
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•  Reduced  Model  A.  From  the  set  of  variables  produced  in  the  aforementioned  model,  ln(total 
duration)  and  microphone  offset  distance  were  used  to  form  a  discriminating  model  that  requires  just  one 
function: 


f^  =  -84.059  +  12.774  ln(Total  Duration)  -  0.523  Microphone  Distance.  (15) 


When  the  frequency  of  the  fj -values  are  plotted  in  histogram  form,  as  shown  in  Figure  18,  four 
clusters  are  produced,  which  represent  the  four  groups  of  interest.  The  histograms  do  not  overlap, 
indicating  correct  classifications  for  all  the  groups  were  achieved.  The  values  on  the  abscissa  are  the 
f  j-values.  Below  that,  the  class  centriod  is  provided  in  two  forms.  The  relative  location  of  the  centroid 
compared  to  the  group  histogram  is  given  first.  Each  group  is  numbered  as  1  for  25-mm  KE,  2  for 
25-mm  HEI,  3  for  Large  KE,  and  4  for  Large  HEAT.  The  f-value  of  the  centriod  is  provided  next.  The 
limit  values  are  the  greatest  value  of  the  group  starting  from  the  left.  These  limits  provide  the  following 
functional  mapping: 


fi  <  -9.02  M910  (25  KE) 

-9.02  <fi  <  -2.39  M793  (25  HEI) 

-2.39  <fi  <  4.78  M724,  M865  (Large  KE) 

fl>4.78  M490,M831  (Large HEAT).  (16) 
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Figure  18.  Feasibility  phase  model  A:  histogram  of  group  classes. 
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•  Reduced  Model  B.  Referring  back  to  the  set  of  seven  variables  that  include  ln(total  duration),  other 
combinations  of  variables  were  considered  that  would  classify  aU  the  rounds  correctly.  Consider  flie 
variables  ln(total  duration)  and  velocity.  The  discriminant  analysis  requires  two  functions  as  shown  in 
equation  17,  resulting  in  98.03%  of  the  data  correctly  classified.  Of  the  misclassified  data,  one  M910 
round  was  classified  as  an  M793,  and  three  M793  roimds  were  predicted  as  M910’s  (see  Table  7). 

fj  =  ^2.451  +  5.994  ln(Total  Duration)  +  0.00129  Velocity 

f2  =  -2.800  +  0.529  ln(Total  Duration)  +  0.00650  Velocity.  (17) 


Table  7.  Classification  Result  for  ln(Total  Duration)  and  Velocity. 


Actual  Group 

No.  of 

Predicted  Group  Membership  || 

Cases 

BSSi 

2 

M793 

3 

M724/M865 

4 

M490/M831 

1.  M910 

36 

35 

97.2% 

t 

2.8% 

2.  M793 

36 

33 

91.7% 

3.  M724/M865 

61 

61 

100% 

4.  M490/M831 

69 

69 

100% 

•  Limitations  on  Model  B.  As  mentioned  earlier,  the  data  gathered  from  the  1-m  microphone  to  the 
target  was  not  reliable.  Therefore,  only  the  data  collected  firom  a  distance  >3  m  was  considered  in  the 
analysis.  With  the  same  variables  ln(total  duration)  and  velocity,  100%  correct  classification  was  achieved 
with  the  model. 


fi  =  -68.539  +  9.719  InCTotal  Duration)  +  0.000994  Velocity 

f2  =  -3.652  -  0.397  ln(Total  Duration)  +  0.00650  Velocity.  (18) 

The  data  plotted  in  (fl,  f2)  space  is  shown  in  a  scatterplot  in  Figure  19.  The  data  point  is  denoted 
with  the  group  number.  The  indicates  the  group  centroid.  Because  of  the  graphics  output  from  SPSS, 
overlapping  points  are  not  showa  The  number  of  points  indicated  on  the  scatterplot  may  not  add  up  to 
the  number  of  cases  actually  analyzed,  and  misclassified  cases  may  be  difficult  to  visualize.  Qassification 
results  are  provided  in  the  table  for  models  where  misclassification  occurs. 
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Figure  19.  Feasibility  phase  model  B:  scatterplot. 

The  territorial  map  in  Figure  20  shows  the  boundaries  for  each  group  in  (f^,  £2)  space.  These 
boundaries  are  constructed  such  that  any  point  in  a  group  in  the  fi-f2  plane  would  be  closer  to  that  group 
centroid  than  any  other  group  centroid.  This  is  useful  for  predicting  the  group  membership  of  an  unknown 
round. 
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Figure  20.  Feasibility  phase  model  B:  territorial  mat 
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The  process  of  selecting  subsets  of  variables  from  the  base  set  continues.  The  remaining  variables: 
peak  of  absolute  pressure  impulse,  negative  peak,  and  peak  of  pressure  impulse  are  difficult  to  obtain  in 
the  field  because  the  microphone  would  have  to  be  more  finely  calibrated. 

•  Reduced  Model  C.  A  third  and  quite  different  reduced  model  was  considered.  Total  duration  is  the 
singip.  variable  that  best  discriminates  among  the  groups  of  roimds.  If  this  variable  were  collected  from 
a  microphone  placed  on  each  side  of  the  target,  collectively  the  rounds  should  be  identified  uniquely. 
Although  our  test  array  of  five  microphones  was  on  the  same  side  of  the  target,  the  data  from  two 
microphones  were  carefully  selected  to  determine  if  discrimination  were  possible  with  nothing  more  than 
two  total  durations.  A  two-function  discrimination  model  was  developed  from  the  6-  and  12-m 
microphones: 

fj  =  -17.186  +  0.0115  Total  Duration  6  +  0.00331  Total  Duration  12 

fj  =  -0.0887  -  0.0317  Total  Duration  6  +  0.02651  Total  Duration  12.  (19) 

The  assumptions  for  an  optimal  model  is  normality  of  the  independent  variables  and  equal  covariance 
matrices  of  each  group.  Box’s  M  test  was  used  to  check  the  equality  of  covariance  matrices  based  on  the 
determinants  of  the  matrices.  Box’s  M  test  is  also  sensitive  to  the  normality  assumption.  Of  the  three 
models  developed  previously,  only  Model  C  did  not  pass  Box’s  M  test  for  equality  of  covariance  matrices. 
This  indicates  that  the  model  is  not  statistically  optimal,  and  problems  may  arise  in  predicting  rounds. 
In  spite  of  this,  the  model  predicted  all  correctly,  as  shown  in  Table  8,  and  the  scatterplot  and  territorial 
map  in  Figures  21  and  22,  respectively,  show  good  separation  among  the  groups. 

•  Verification  of  the  Feasibility  Models.  Each  of  the  three  models  described  in  the  previous  section 
represents  the  best  classification  rate  since  all  of  the  data  was  used  to  develop  the  model.  Often,  live  fire 
tests  are  conducted  to  verify  that  the  model  is  truly  adequate.  This  verification  was  expected  to  be 
performed  later  in  the  program  when  all  round  types  had  been  incorporated  into  the  final  model(s).  To 
check  the  adequacy  of  the  current  feasibility  models,  a  portion  of  the  data  was  taken  for  model 
development.  The  remaining  data  were  reserved  for  testing  the  model,  a  form  of  goodness-of-fiL 
Stratified  random  sampling  of  two-thirds  of  data  was  implemented  for  adjusting  the  coefficients  associated 
with  the  variables  in  each  of  the  three  models  listed  previously.  Stratification  was  performed  by  randomly 
selecting  one  round  from  each  round  type  at  each  range  to  be  reserved  for  testing  the  model  adequacy. 
The  other  rormds  were  used  for  the  model  development 
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Figure  21.  Feasibflitv  phase  model  C:  scatterplot. 
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Figure  22.  Feasibility  phase  model  C:  territorial  map. 
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All  three  resulting  models  performed  well.  Each  adjusted  model  correctly  predicted  aU  of  the  rounds 
used  to  develop  the  model,  and  each  model  correctly  predicted  each  of  the  reserved  rounds.  Additionally, 
a  set  of  theoretical  data  was  created  in  a  blind  study  to  test  the  adequacy  of  Model  C  with  two  total 
durations.  Twelve  theoretical  total  duration  values  were  computed  for  the  various  round  types  as  if  the 
microphones  collecting  acoustical  signatures  were  placed  on  each  side  of  the  target  These  values  were 
placed  in  the  discriminatory  function.  The  resulting  discriminant  values  correctly  predicted  the  group 
membership  for  each  round  type,  as  seen  in  Table  9. 


Table  8.  Qassification  Result  for  Two  Total  Durations 


Actual  Group 

No.  of 

Predicted  Group  Membership 

Cases 

2 

M793 

3 

M724/M865 

4 

M490/M831 

1.  M910 

9 

9 

100% 

2.  M793 

9 

9 

100% 

3.  M724/M865 

15 

15 

100% 

4.  M490/M831 

18 

18 

100% 

Table  9.  Verification  Results 
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•  Feasibility  Test  Conclusion.  Acoustical  signatures,  particularly  total  duration  of  the  N-Wave  curve, 
can  be  used  to  discriminate  between  the  following  round  types: 

small-  and  large-caliber  rounds 

-  large  (105-  or  120-mm)  HE  and  large  KE  training  rounds 

-  25-mm  HEI  and  KE  training  rounds. 

With  the  use  of  discriminant  analysis,  three  models  that  correctly  classify  100%  of  the  rounds  have 
been  developed.  These  models  contain  the  following  variables: 

-  total  duration  and  microphone  offset  distance 
total  duration  and  velocity  (>3  m) 

-  two  total  durations. 

5.2.5  Location  and  Velocity  Data.  During  the  feasibility  test,  the  six  microphones  were  set  up  in 
two  triangular  arrays  as  shown  in  Figure  23  (see  Appendix  C  for  exact  microphone  locations).  The  time 
differences,  microphone  and  target  locations,  and  meteorological  data  were  processed  through  the  Stallings 
algorithm  discussed  in  3.2.3.  The  data  were  first  reduced  for  velocity  prediction.  The  Stallings  algorithm 
predicted  the  velocity,  on  the  average,  within  2%  of  the  actual  velocity  verified  by  radar.  The  locations 
were,  on  the  average  20  cm  from  the  actual  video  scored  locations  (see  Appendix  G).  The  error  was 
mostly  due  to  microphone  placement.  The  projectile  hits  did  not  go  through  the  center  of  the  microphone 
array.  The  microphones  were  purposely  placed  low  to  the  ground  to  keep  them  from  being  hit  and  to  see 
what  errors  would  be  generated. 

Six  M831  projectiles,  fired  from  the  1,000  m  position,  were  analyzed  to  study  the  effects  of  input  data 
on  the  accuracy  of  the  predicted  location.  Table  10  shows  how  a  change  in  each  input  value  of 
meteorological  data,  fall  angle,  and  gun  position  affects  the  projectile  location. 

The  error  associated  with  the  microphone  positions  depends  on  the  contribution  of  diat  microphone 
when  determining  the  location.  If  the  projectile  flies  close  to  a  microphone,  the  algorithm  weights  that 
microphone’s  data  such  that  it  provides  a  larger  contribution  to  the  output  location.  Generally,  the  average 
survey  error  of  all  the  microphones  will  be  close  to  the  average  scoring  error.  The  survey  error  was 
typically  around  1  mm  for  each  microphone  with  the  use  of  a  precision  theodolite  and  an  experienced 
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Input  Values  Change  for  Each  Additional  Change  in 

Location 

(mm) 


Meteorological  Data 


Temperature  1°  3 

Wind  1  m/s  4 

Humidity  Any  Humidity  3  max 


Ballistics 


FaU  Angle  1  mil  3.6 


Gun  Position 


X  1  m  5 

Y  1  m  0.5 

Z  1  m  4 
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surveyor.  Errors  associated  with  microphone  times  also  depended  on  that  microphone  s  contribution. 
There  is  more  of  an  effect  on  location  for  slower  projectiles  since  their  shock  waves  are  more  oblique. 
The  errors  associated  with  meteorological  data  and  fall  angle  are  considered  relatively  small.  The  major 
downfall  with  the  Stallings  technique  is  the  dependence  on  knowing  the  gun’s  position.  Elevation 
changes,  Z,  are  not  that  great  on  a  training  range.  The  major  errors  come  from  not  knowing  an  accurate 
range  and  deflection  of  the  gun  position  relative  to  the  target.  The  Stallings  technique  was  eliminated  as 
a  candidate  for  RDS  location  algorithm  because  it  was  slow  and  had  large  potential  gun  position  errors. 

5.3  Collection  Phase. 

5.3.1  Live-Fire  Data.  The  next  phase  of  testing  began  with  the  purpose  of  expanding  the  existing 
data  set  to  include:  7.62-mm  ball  (M59)  and  tracer  (M62)  rounds,  12.7-mm  ball  (M33)  and  tracer  (M17) 
rounds,  and  three  35-mm  TPGID  (old  HEAT  DM18,  new  HEAT  DM68,  and  KE  simulator)  rounds.  This 
phase  is  designated  as  the  collection  phase.  Randomization  of  rounds  for  testing  took  place  as  much  as 
possible  with  the  tight  test  schedule  and  small  budget 

Since  changing  the  firing  platform  for  each  round  at  each  range  was  time-consuming,  several  rounds 
of  the  same  size  were  fired  at  one  range  before  moving  the  platform  to  another  range.  The  goal  was  to 
divide  the  sample  in  half  so  that  two  groups  of  rounds  of  the  same  type  were  fired  on  two  different  days 
to  average  out  any  effects  due  to  weather  conditions.  This  was  not  always  accomplished.  The  DM18  and 
DM68  round  were  late  arriving.  The  DM68  and  the  TPGID  KE  simulators  had  not  been  type  classified. 
The  7.62-mm  and  12.7-mm  rounds  were  fired  first.  Ball  and  tracer  rounds  were  completely  randomized. 
When  the  35-mm  KE  was  type  classified,  it  was  randomized  with  the  DM18.  The  DM68  was  tested  with 
the  DM18  on  the  last  two  days.  However,  the  duration  data  for  the  last  two  days  was  not  available  for 
the  analyses.  Table  14  lists  the  number  of  rounds  reduced  that  were  included  firom  each  type  in  the 
analysis. 

The  ball  and  tracer  rounds  are  listed  separately  in  Table  11,  but  there  was  little  interest  in 
discriminating  between  them  in  this  study.  There  were  some  differences  in  total  duration  and  velocity. 
The  average  velocity  of  the  tracer  rounds  at  each  range  was  greater  than  the  ball  rounds  for  both  7.62  mm 
and  12.7  mm.  The  average  total  duration  of  the  tracer  rounds  was  approximately  13  ps  greater  than  the 
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Table  11.  Number  of  Rounds  Analyzed  in  Collection  Phase 


Distance 

7.62  mm 

12.7  mm 

35  mm 

(m) 

M59 

M62 

M33 

M17 

DM68 

KE 

Simulator 

3 

4 

5 

4 

3 

3 

3 

4 

— 

— 

— 

— 

800 

3 

3 

— 

— 

— 

— 

1,000 

— 

— 

4 

4 

4 

3 

1,700 

— 

— 

— 

— 

3 

3 

ball  roimds  for  both  7.62  mm  and  12.7  mm  with  the  exception  of  the  7.62  mm  at  800  m.  For  this 
combination,  the  average  ball  round  was  230  ps  greater  than  the  tracer  round.  More  investigation  is 
required  to  determine  what  caused  this  large  difference  at  800  m.  See  Appendix  E  for  tables  of  total 
duration  and  velocity  values.  For  the  remaining  analyses,  the  ball  and  tracer  rounds  of  7.62  mm  will  be 
treated  as  one  group  as  well  as  the  12.7-mm  rounds. 

5.3.2  Discriminate  Analysis  for  Empirical  Model  Formation.  The  entire  data  set,  from  both  the 
feasibility  and  collection  phases,  had  a  total  of  8  groups:  7.62  mm,  12.7  mm,  M910,  M793,  35-mm 
TPGID  KE  simulator,  DM18  and  DM68,  M724  and  M865,  and  M490  and  M831.  A  cursory  analysis 
performed  on  the  data  before  the  complete  data  set  was  reduced  indicated  that  prediction  would  be  good 
overall.  However,  within  the  TPGID  groups,  only  66%  of  the  rounds  could  be  classified  correctly. 

A  total  of  34%  of  the  DM18/DM68  rounds  were  predicted  as  TPGID  KE  simulator  rounds.  Likewise, 
34%  of  the  TPGID  KE  simulator  rounds  were  predicted  as  DM18/DM68  rounds.  An  encouraging  point 
was  that  fiiese  rounds  were  stiU  classified  as  35-mm  TPGID.  Upon  examination  of  these  two  round  types, 
we  found  that  the  velocities  of  KE  training  rounds  at  each  range  were  on  the  average  250  m/s  slower  than 
the  120-mm  KE  training  rounds.  To  continue  the  analysis,  the  35-mm  TPGID  rounds  were  combined  into 
one  group  resulting  in  seven  groups  for  discrimination. 

•  Model  I:  Total  Duration,  Microphone  Distance,  and  Velocity.  Based  on  the  three  models  developed 
during  the  feasibility  phase,  discriminant  functions  were  formed  first  with  variables  ln(total  duration)  and 
microphone  distance.  Unlike  the  one-function  model  formed  for  the  four  groups  in  that  phase,  the  model 
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for  seven  groups  required  two  functions.  However,  misclassification  was  high,  particularly  between  the 
M910  group  and  12.7  mm.  The  misclassification  rate  was  17%. 

The  next  attempt  was  to  form  a  model  with  variables  ln(total  duration)  and  velocity,  as  in  Model  B. 
This  model  produced  88%  correct  classification  of  the  groups.  Only  the  M490/M831  rounds  were  without 
misclassification.  TPGID  was  the  most  difficult  of  all  the  groups  to  classify  correctly. 

Finally,  all  three  variables — ^ln(total  duration),  microphone  offset  distance,  and  velocity — ^weie  included 
to  produce  a  two-function  model  (see  equation  20).  This  resulted  in  a  good  classification  rate  of  97.85% 
as  shown  in  Table  12.  From  the  scatterplot  of  the  data  in  (f^,  f.^  space  (Figure  24),  one  can  see  why 
several  7.62-mm  rounds  (group  1)  were  classified  as  12.7  mm  (group  2);  the  functional  values  were  closer 
to  the  centroid  of  group  2.  Likewise,  several  M910  rounds  (group  3)  were  misclassified  as  12.7  mm 
(group  2)  since  their  values  were  closer  to  that  centroid.  The  misclassified  rounds  could  be  easily 
identified  by  overlajdng  the  territorial  map  (Figure  25)  on  the  scatteiplot  (Figure  24).  (The  pattern 
numbers  in  the  territorial  map  form  the  border  lines  between  the  groups.) 

fj  =  -87.510  +  13.594  ln(Total  Duration) 

+  0.00413  Velocity  -  0.529  Microphone  Offset  Distance 
f2  =  6.592  -  1.928  ln(Total  Duration) 

+  0.00650  Velocity  +  0.0671  Microphone  Offset  Distance.  (20) 


This  model  was  tested  with  data  collected  during  the  verification/validation  phase.  The  groups’ 
centroids  listed  in  Table  13  are  necessary  for  determining  the  group  membership  of  the  unknown  rounds 
fired.  The  boundary  conditions  are  set  at  fi(-20,  20)  and  f2(-4,  4).  Any  rounds  producing  values  of  fj 
and  f2  outside  these  boundaries  were  classified  as  an  unknown  roimd. 

•  Model  II:  Three  Variables  at  6  m.  As  discussed  earlier,  the  field  equipment  array  had  two  high- 
fidelity  microphones  placed  6  m  from  the  center  of  the  target.  The  same  variables  used  in  the  previous 
model  were  considered,  but  the  coefficients  were  based  on  the  6-m  microphone  data  only.  The 
classification  rate  for  this  model  (see  equation  21)  was  close  to  the  Model  I  with  97.12%  correct  as  shown 
in  Table  14.  The  centroids,  scatterplot,  and  territorial  map  are  provided  (see  Table  15  and  Figures  26  and 
27,  respectively).  The  boundary  conditions  are  fj(-20,  20),  f2(-4.  4),  and  f3(-4,4).  This  model  was 
denoted  as  Empirical  Model  II  and  tested  in  the  verification/validation  phase. 
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fj  =  -95.122  +  14.681  ln(Total  Duration) 

+  0.00432  Velocity  -  0.5122  Microphone  Offset  Distance 
£2  =  5.512  -  1.878  ln(Total  Duration) 

+  0.00618  Velocity  +  0.258  Microphone  Offset  Distance 
f3  =  -9.133  -  0.202  ln(Total  EHiration) 

-  0.00041  Velocity  +  1.987  Microphone  Offset  Distance.  (21) 


Table  12.  Qassification  Results  for  ln(Total  Duration),  Microphone  Offset  Distance  and  Velocity 


Actual  Group 

No.  of 

Predicted  Group  Membership  || 

Cases 

1 

7.62  mm 

2 

12.7  mm 

3 

M910 

4 

M793 

5 

TPGID 

6 

M724/ 

M865 

7 

M490/ 

M831 

1.  7.62  mm 

59 

54 

91.5% 

•5 

'8.5% 

2.  12.7  mm 

65 

65 

100.0% 

3.  M910 

29 

3  ■ 
m3% 

26 

89.7% 

4.  M793 

29 

29 

100.0% 

5.  TPGID 

29 

29 

100.0% 

6.  M724/M865 

56 

56 

100.0% 

7.  M490/M831 

65 

65 

100.0% 

NOTE:  Cases  correctly  classified:  97.85%, 


Table  13.  Discriminant  Functions  Group  Centroids 


Group 

^2 

M62/M59 

-13.170 

-0.621 

M17/M33 

-8.676 

-0.609 

M910 

-6.346 

2.030 

M793 

-1.686 

-0.990 

DM18/DM68/rPL-T 

1.800 

0.581 

M724/M865 

9.119 

2.207 

M490/M831 

14.447 

-1.810 
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Figxire  24.  Model  I:  scatterplot 
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Figure  25.  Model  I:  territorial  map. 
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Table  14.  Model  n  aassification  Results 


NOTE:  Cases  correcdy  classified:  97.12%. 


Table  15.  Model  II  Group  of  Centroids 


Group 

fi 

h 

7.62 

-14.276 

-0.582 

-0.191 

12.7 

-9.323 

-0.583 

0.011 

M910 

-6.790 

1.826 

0.145 

M793 

-1.598 

-0.893 

1.060 

TPGID 

2.035 

0.474 

-0.474 

M724A1865 

9.811 

2.089 

M490/M831 

15.726 

-1.697 

-0.100 

•  Model  III:  Two  Total  Durations.  The  final  empirical  model,  equation  22,  includes  two  total 
duration  variables  as  constmcted  in  the  feasibility  phase  (Model  C).  The  data  from  the  9-m  microphone 
were  the  most  discriminating,  and  the  3-m  microphone  data  were  the  second  most  discriminating. 
Although  an  effort  was  made  to  use  the  data  measured  at  6  m,  the  data  at  3  m  and  9  m  were  the  least 
correlated  as  weU  as  the  most  discriminating.  The  correlation  between  the  total  duration  values  at  9  m 
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Figure  26.  Model  II:  scatterolot. 
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Figure  27.  Model  II:  territorial  map. 
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and  6  m  is  0.94.  This  high  correlation  between  explanatory  variables  can  cause  a  problem  known  as 
multicollinearity  in  which  the  coefficients  are  very  sensitive.  The  goal  is  to  obtain  coefficients  that  are 
robust  so  that  the  model  wiU  be  valid  under  various  conditions. 

=  -16.158  -  0.000608  Total  Duration  (3-m  Mic) 

+  0.0199  Total  Duration  (9-m  Mic) 
f2  =  0.200  +  0.0316  Total  Duration  (3-m  Mic) 

-  0.0233  Total  Duration  (9-m  Mic).  (22) 

The  correct  classification  rate  for  this  model  was  lower  than  the  other  two  models  at  89.69%  (see 
Table  16).  The  difficulty  in  discriminating  occurred  for  rounds  35  mm  and  smaller.  The  scatterplot  and 
territorial  map  (Figures  28  and  29)  respectively  show  this  from  the  location  of  the  group  centroids.  The 
centroids  of  three  smallest  rounds — iJ&l  mm,  12.7  mm,  and  M910  are  very  close  as  are  the  centroids  of 
M793  and  TPGID.  The  boundary  values  for  the  discriminant  functions  below  are  fj(-15,25)  and  f2(-4,4). 
The  group  centroids  are  provided  in  Table  17. 

The  three  models  developed  in  this  collection  phase  were  tested  against  a  new  set  of  data  during  the 
verificationA^alidation  test  along  with  the  brassboard  system.  Those  results  are  given  in  section  9. 

5.3.3  Location  and  Velocity  Data  Reduction.  The  microphones  were  placed  in  the  H  configuration, 
two  lines  of  three  microphones  each  (Figure  30).  For  this  study,  the  intersection  of  Hypertwlas  of 
Ambiguity  method  was  used  to  calculate  offset  distance.  The  predicted  locations  were  on  average  24  cm 
away  firom  the  actual  video  scored  locations,  and  the  predicted  velocities  were  within  3%  (see 
Appendix  H).  Errors  in  location  and  velocity  are  mostly  due  to  the  projectile’s  incoming  angle.  The 
projectile  must  travel  parallel  to  the  centerline  connecting  the  two  center  microphones  for  the  calculation 
to  be  accurate.  As  this  incoming  angle  increases,  the  calculated  errors  increase  quadratically. 

6.  RDS  ALGORITHM  DEVELOPMENT 

6.1  Theoretical  Formulation  Technique.  The  total  duration  formula  described  in  equation  2  was  based 
on  test  data  with  small-caliber  projectiles  fired  from  M14  and  AR15  rifles  (Kalb,  unpublished).  The 
constant  value  of  1.82  was  used  in  the  fonnula  for  smaU-caliber  ogival  bodies.  A  new  variable  was 
created  for  this  constant  and  named  the  acoustic  shape  constant,  K.  During  phase  I  and  n  of  the  data 
collection  testing,  each  K  was  solved  for  by  rearranging  Schochko’s  (1966)  formula  for  K  and  plugging 
in  the  known  values. 
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Table  16.  Model  III  Qassification  Results 


Actual  Group 

No.  of 

Predicted  Group  Membership 

- 

Cases 

1 

7.62  mm 

2 

12.7  mm 

3 

M910 

4 

M793 

5 

TPGID 

6 

M724/ 

M865 

7 

M490/ 

M831 

1.  7.62  mm 

17 

15 

88.2% 

2 

11.8% 

2.  12.7  mm 

15 

12 

80.0% 

3 

20.0% 

3.  M910 

9 

1.  "t.:" 

7 

77.8% 

4.  M793 

9 

6 

66.7% 

3 

33.3% 

5.  TPGID 

18 

18 

100.0% 

6.  M724/M865 

14 

100.0% 

7.  M490/M831 

15 

100.0% 

NOTE:  Cases  correctly  classified:  89.69%. 
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Figure  28.  Model  III:  scatterplot. 
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Figure  29.  Model  III:  territorial  map. 


Table  17.  Model  HI  Group  of  Centroids 


Group 

h 

7.62 

-10.574 

0.198 

12.7 

-8.561 

0.515 

M910 

-8.036 

-0.132 

M793 

-3.415 

-0.581 

TPGID 

-1.444 

0.037 

M724/M865 

6.787 

-0.878 

M490/M831 

22.815 

0.464 

(23) 


An  averaged  K  value  was  then  determined  over  the  velocities  and  offset  distances  tested.  Since  each 
projectile  has  unique  physical  properties  (diameter  and  length)  and  K  value,  equation  23  was  rearranged 
to  equate  the  projectile  physical  dependents  on  one  side  with  the  measured  values  on  the  other.  The 
projectile  dependents  were  grouped  together  and  labeled  the  projectile  shape  value  (PSV).  The  formula 
becomes 

PSV  =  _2_K  =  .  (24) 

l1/4  (O)^^^M 
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LM6 


Figure  30.  Data  collection  test  microphone  array  and  target  setup. 


Table  18  contains  the  empirical  PSV  and  category  for  each  projectile  type  tested.  The  measured 
values  were  input  into  equation  3  after  the  raw  analysis.  Next,  the  measured  were  compared  to  the 
empirical  values.  The  system  identiftes  the  projectile  t)^  by  calculating  the  smallest  measured  empirical 
PSV  differential. 


After  the  data  were  analyzed  and  the  measured  values  were  input  into  equation  24,  the  next  step  was 
to  check  to  see  how  these  values  compared  with  tabular  values.  This  was  accomplished  by  finding  the 
difference  of  the  measured  PSV  fiom  microphone  1  and  the  theoretical  PSV  and  selecting  the  projectile 
type  that  best  agrees. 


PSVl  =  ^_K  = 

L1/4 


(T  Mic  l)a(M^-l)^^ 
(O  Mic 


(25) 
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Table  18.  Average  Projectile  Shape  Values 


Gun 

Caliber 

(rrun) 

Projectile  Type 

Diameter 

(mm) 

Length 

(m) 

Average 

K 

Average 

PSV 

Projectile 

Category 

(mm) 

7.62 

M59  (baU) 

7.82 

0.032 

2.08 

0.038 

7.62 

7.62 

M62  (tracer) 

7.82 

0.034 

2.16 

0.039 

7.62 

12.7 

M33  (baU) 

12.95 

0.059 

2.11 

0.056 

12.7 

12.7 

M17  (tracer) 

12.95 

0.059 

2.21 

0.058 

12.7 

25 

M910  (KE-TP) 

16.20 

0.076 

1.99 

0.062 

25  KE 

25 

M793  (HEI-TP) 

25.00 

0.115 

2.18 

0.093 

25  KE 

35 

TPL  (KE-TPGID) 

35.00 

0.180 

2.04 

0.110 

35 

35 

DM18  (HE-TPGID) 

35.00 

0.180 

2.12 

0.114 

35 

35 

DM68  (HE-TPGID) 

35.00 

0.180 

35 

105 

M724  (KE-TP) 

63.50 

0.256 

1.72 

0.154 

105/120  KE 

120 

M865  (KE-TP) 

38.00 

0.468 

3.71 

0.170 

105/120  KE 

105 

M490  (HEAT-TP) 

105.00 

0.654 

2.50 

0.291 

105/120  KE 

120 

M831  (HEAT-TP) 

120.00 

0.844 

2.37 

0.296 

105/120  KE 

PSVl  =  PSV  for  HMl 
Remainder!  =  abs  (PSVl-  PSV). 


The  remainders  were  ranked  from  the  smallest  to  the  largest.  In  the  same  manner,  the  PSV  for 
microphone  2  is  solved  for 


PSV2  =  _£_K  =  CT  Mic  2)a(M^-l)^^ 
L  (O  Mic  2)^/^M 


(26) 


PSV2  =  PSV  for  HM2 
Remainder2  =  abs(PSV2-  PSV). 

The  remainders  are  ranked  from  smallest  to  largest. 

If  both  rankings  agree,  that  is  the  projectile  type  predicted.  If  not,  there  are  three  possible  cases  where 
the  ranks  are  different.  In  case  1,  PSVl  is  closer  and  greater  than  round  type  A,  and  PSV2  is  closer  and 
lesser  than  rormd  type  B.  This  yields  two  positive  remainders  for  round  type  A  and  two  negative 
remainders  for  round  type  B.  Round  type  B  is  selected.  In  case  2,  PSVl  is  closer  and  greater  than  round 
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type  A,  and  PSV2  is  closer  and  lesser  than  round  type  B.  Round  type  A  is  selected.  In  case  3,  PSVl  is 
closer  and  lesser  than  round  type  A,  and  PSV2  is  closer  and  greater  than  round  type  B.  Take  the  average 
of  PSVl  and  PSV2  and  reevaluate  the  remainder  to  find  the  round  type. 

Figure  31  shows  the  projectile  terminal  velocity  plotted  against  PSV.  It  can  be  seen  that  some  PSV’s 
from  different  round  types  overlap  at  extreme  conditions.  AH  of  the  projectiles  fired  have  well-defined 
velocities  at  known  ranges  (see  Table  2).  If  a  PSV  is  computed  and  there  are  two  possibfiities,  the  correct 
round  type  can  be  selected  if  referenced  to  the  velocity  axis.  For  instance,  the  12.7-mm  and  25-mm  KE 
rounds  have  similar  PSV’s.  If  the  terminal  velocity  of  the  projectile  fired  was  greater  than  850  m/s,  the 
12.7-nun  round  could  be  ruled  out  as  a  possible  choice  since  it  cannot  achieve  that  velocity.  This  model 
was  called  the  theoretical  model. 

6.2  Hybrid. 

6.2.1  Discriminatory  Analysis  for  PSV.  The  quasi-empirical  formulation  for  PSV  is  a  function  of 
the  same  significant  variables  used  in  the  empirical  model,  total  duration,  microphone  offset  distance,  and 
velocity.  The  theoretical  and  the  empirical  results  were  combined  to  form  a  "hybrid"  model.  The  PSV 
will  be  used  in  the  discriminant  analysis. 

From  the  data  collected  in  the  feasibility  phase  and  data  collection  phase,  the  projectile  shape  was 
tested  to  determine  if  it  could  serve  as  a  discriminator  of  the  seven  groups.  Discriminant  analysis  was 
implemented  with  the  model  based  on  the  data  collected  from  the  6-m  microphone.  This  distance  was 
chosen  since  this  was  the  same  location  for  the  high-fidelity  microphones  as  planned  for  the  live-fire 
validation  test.  The  classification  of  the  remaining  data  was  performed  for  a  sensitivity  check  of  the  model. 

The  model  with  projectile  shape  alone  resulted  in  8  misclassifications  and  94  correct  classifications 
for  a  success  rate  of  92.16%.  The  shape  variable  was  added  to  the  database  of  variables  to  determine  if 
it  were  a  superior  discriminant  variable  when  compared  to  the  other  variables.  Using  the  discriminant 
analysis  models,  PSV  was  found  to  improve  the  performance  of  the  classifications.  PSV’s  were  first  to 
enter  the  model  with  ln(total  duration)  following.  The  correct  classification  rate  dropped  to  85.29%. 
Although  total  duration  was  allowed  to  enter  the  model  based  on  the  tolerance  criterion,  the  pooled  group 
correlation  between  shape  and  ln(total  duration)  was  0.814. 
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Velocity  became  the  variable  to  enter  the  model  after  PSV.  This  was  intuitively  appealing,  since  PSV 
was  basically  replacing  ln(total  duration)  in  the  earlier  model  fonned  with  ln(total  duration)  and  velocity. 
If  we  compare  plots  of  ln(total  duration)  vs.  velocity  (Figure  31)  with  PSV  vs.  velocity  (Figure  32),  we 
notice  that  the  relative  positions  of  the  round  groups  are  similar.  The  plot  with  PSV  forms  tighter 
groupings. 

The  two-function  model  with  variables  PSV  and  velocity,  equation  27,  provided  a  99.02%  correct 
classification  rate  from  the  6-m  data,  with  one  25-mm  KE  misclassified  as  a  12.7-mm  round.  The 
discriminatory  functions  Fj  and  F2  plotted  in  Figure  33  resembles  the  data  in  the  original  umts  of  PSV 
vs.  velocity  Figure  32. 


Fj  =  -22.757  +  144.233  PSV  +  0,00500  Velocity 

F2  =  -3.904  -  10.733  PSV  +  0.00632  Velocity.  (27) 


6.2.2  Model  Implementation  for  Prediction.  A  model  was  developed  that  separated  the  seven  types 
of  rounds  with  good  success.  The  trae  test  of  the  adequacy  of  the  model  would  come  when  additional 
rounds  were  fired  in  the  verificationAalidation  phase  (section  9).  How  is  the  two-function  model 
implemented  once  the  data  are  collected?  First,  the  centroids  of  each  group  were  calculated  for  both 
functions  as  provided  in  Table  19. 


Table  19.  Group  Centroids 


Group 

Fl 

F2 

7.62 

-14.765 

-1.296 

12.7 

-11.495 

-0.909 

M910 

-8.781 

1.716 

M793 

-5.610 

-0.787 

TPGID 

-1.338 

0.848 

M724/M865 

7.800 

2.678 

M490/M831 

25.093 

-1.586 

Second,  boundary  conditions  for  Fj  and  F2  were  determined  to  guard  against  signals  outside  of  the 
target  area  of  interest.  This  is  not  really  necessary  since  the  round  location  method  would  not  allow  the 
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Shape  Value  versus  Velocity 


Figure  32.  PSV  value  vs.  velocity. 
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target  to  drop  if  the  round  did  not  hit  the  target  But  to  implement  the  round  discrimination  method 
independently  of  the  round  location  method,  boundary  conditions  are  necessary.  The  boundary  area  is  a 
simple  rectangular  area  in  Fj  and  F2  space.  The  minimum  and  maximum  values  for  each  function 
is:  Fi  (-17.0,  29.0)  and  F2  (-4.5,  5.5). 

When  an  unknown  round  is  fired  at  the  target,  the  total  duration,  microphone  offset  distance,  and 
velocity  will  be  collected,  and  the  PSV  wiU  be  computed.  Then  values  of  Fj  and  F2  are  computed  for 
this  round.  The  distance  is  calculated  between  the  unknovm  roimd  and  each  group  centroid.  The 
minimum  distance  will  dictate  the  classification  of  the  unknovm  round. 

Because  two  microphones  will  be  placed  on  each  side  of  the  target  in  the  verification/validation  test 
to  collect  total  duration  data,  two  pair  of  discriminant  scores  will  be  computed,  one  for  the  left  microphone 
and  one  for  the  right  It  is  possible  that  the  minimum  distance  of  the  discriminant  score  for  the  left 
microphone  classifies  the  unknown  round  into  a  different  group  as  the  right  In  that  case,  the  final  group 
assignment  would  be  the  smallest  of  those  two  minimum  distances. 

7.  SOFTWARE  DEVELOPMENT 

A  microprocessor  was  determined  to  be  the  most  effective  way  to  process  the  input  data  using  the 
model,  perform  decisions,  and  execute  commands.  The  RDS  needed  to  take  the  data  developed  by  the 
sensors  after  a  roimd  was  sensed.  The  signals  for  location  and  acoustic  signature  were  to  be  processed, 
and  the  values  of  that  process  utilized  by  the  prediction  model.  Once  a  prediction  was  made  as  to  round 
type  and  location  of  hit,  commands  based  on  a  decision  of  "kill,"  "hit  but  no  kiU,"  and  "no  hit"  would 
command  die  RETS  system.  A  486  PC/104  embedded  computer  was  selected  as  the  microprocessor  due 
to  the  ease  of  programming,  ability  to  use  peripherals  (disk  drives,  monitors,  and  keyboards),  and 
accessibility  to  remote  sites.  Additional  hardware  consisted  of  power  supplies,  a  floppy  drive,  a  hard 
drive,  and  three  electronic  interface  boards.  The  boards  were  designed  and  fabricated  at  ARL.  They  were 
the  location  interface  board,  the  sensor  interface  board,  and  the  A/D  board.  To  control  the  computer, 
special  software  was  generated. 

The  software  for  the  RDS  system  was  developed  in  Microsoft  C  to  run  under  an  IBM  PC  or 
compatible  platform.  The  I/O  devices  used  to  collect  the  signal  wave  forms  and  counter  values  were 
memory-mapped  devices  which  enabled  a  high  rate  of  data  transfer.  Two  data  files  were  created  to  run 
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the  program.  The  first  file  is  called  "locate.dat"  and  contained  microphone  and  target  coordinates  needed 
to  calculate  round  location.  The  second  file  was  called  "shapes.dat,"  and  it  contained  "shape  values 
needed  to  categorize  rounds  for  identification. 

The  program  began  with  a  series  of  questions  to  initialize  the  system.  The  first  question  asked 
whether  or  not  results  will  be  transmitted  back  to  a  remote  location.  The  second  question  asked  whether 
or  not  raw  data  and  results  will  be  saved  onto  the  hard  disk.  If  data  and  results  are  to  be  saved,  a  file 
name  was  specified.  Finally,  the  program  asked  what  type  of  target  was  being  used  for  the  training 
exercise. 

The  system  waited  for  a  hardware  trigger  signal  from  any  one  of  the  location  microphones.  The 
trigger  was  caused  by  the  shock  fiom  a  projectile  flying  by  any  of  the  location  microphones.  While  the 
system  polled  the  trigger  line,  a  2,000-point  circular  buffer  was  being  continuously  filled  with  data 
collected  through  the  high-fidelity  microphones  and  into  the  dual-chaimel  analog-to-digital  board  at  a  rate 
of  250  kHz  per  channel. 

When  triggered,  1,000  points  of  shock  wave  data  were  collected  from  each  high-fidelity  microphone. 
A  total  of  12,000  points  of  shock  wave  data  were  collected  from  each  high-fidelity  microphone  when  the 
circular  buffer  data  were  added.  The  circular  buffer  contained  pretrigger  data  which  ensured  the  entire 
wave  form  of  the  shock  wave  was  preserved  in  case  the  shock  wave  reached  the  high-fidelity  microphones 
before  the  location  microphones  which  triggered  the  collection. 

Counter  values  were  next  read  in  from  the  quartz  counter  board.  There  were  five  counters,  and  each 
measured  the  relative  delay  of  the  shock  wave  as  it  passed  between  adjacent  location  microphones. 
Figure  34  shows  how  the  counters  were  set  up  in  relation  to  the  location  microphones.  Each  counter  was 
controlled  by  a  pair  of  location  microphones;  Coimter  0  by  LMl  and  LM2,  Counter  1  by  LM2  and  LM3, 
Counter  2  by  LM2  and  LM5,  Counter  3  by  LM4  and  LM5,  and  Counter  4  by  LM5  and  LM6.  Counters 
started  counting  when  either  one  of  its  controlling  microphones  encountered  the  shock  wave  and  stopped 
when  its  other  controlling  microphone  encountered  the  shock  wave  (both  controlling  microphones  could 
start  or  stop  the  counter).  Knowing  the  rate  at  which  the  counters  counted  and  the  counter  values  when 
they  stopped,  the  time  delays  of  the  shock  wave  between  microphones  are  then  determined.  In  addition 
to  the  counters,  four  data  bits  were  read  to  determine  which  of  the  microphones  in  each  pair  started  each 
counter  counting.  The  sign  of  the  counters  was  adjusted  accordingly  finom  these  four  data  bits  to  show 
if  the  round  came  fi'om  the  left  or  the  right. 
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Using  the  counter  values,  the  location  of  the  round  on  the  target,  the  offset  distances  of  the  round  from 
the  high-fidelity  microphones,  and  the  round  velocity  were  calculated.  The  round  velocity  was  calculated 
first  using  Coimter  2.  The  location  and  offset  distances  were  calculated  using  velocity  and  Coxmters  0, 
1,  3,  and  4, 

Using  the  shock  wave  length,  round  velocity,  and  round  offset  distance  from  the  high-fidelity 
microphones,  the  round  type  was  determined.  To  calculate  the  shock  wave  length,  the  wave  form  data 
was  first  passed  through  a  five-point  moving  average  smoothing  filter  to  eliminate  random  noise.  Next, 
the  data  was  truncated  to  500  points  to  eliminate  the  data  points  before  and  after  the  shock  wave.  The 
shock  wave  length  was  found  by  performing  a  point-by-point  differentiation  of  the  remaining  500-point 
wave  form,  picking  out  the  two  peaks.  The  separation  between  the  two  peaks  corresponded  to  the  shock 
wave  length.  Figures  35a  through  35d  graphically  illustrate  the  process  of  finding  the  shock  wave  length. 

After  the  round  location  and  type  had  been  detemiined,  the  program  checked  if  the  round  hit  inside 
the  kill  area  and  if  the  round  type  could  kill  the  target  type.  If  both  of  these  conditions  were  true,  the 
program  would  signal  the  RETS  to  lower  the  target.  The  RDS  would  accept  a  signal  from  the  RETS 
vibration  sensor,  which  indicates  a  roimd  striking  the  target  The  round  type,  location,  and  velocity  was 
radioed  back  to  a  remote  site  via  a  modem  and  hand-held  radio.  The  system  was  reset  to  the  initial 
values. 
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8.  BRASSBOARD 


8.1  Theory  of  Operation.  The  RDS  brassboard  was  a  stand-alone  system  that  directly  interfaced  with 
the  RETS.  The  RDS  sent  a  kill  signal  to  the  RETS  when  the  system  discriminates  the  round.  It  is  a 
modular  PC-based  system,  and  it  has  its  own  radio  communication  link  to  the  range  control  tower.  The 
RDS  brassboard  was  made  of  off-flie-shelf  components  except  for  the  two  electronic  interface  boards.  The 
system  operates  on  llO-V/60-Hz  or  220-V/50-Hz  AC  power. 

Prior  to  running  the  RDS,  software  initialization  was  required  with  each  power  up.  Temperature, 
target  type,  microphone  positions  relative  to  the  target,  and  kill  zone  information  must  be  entered  by  the 
user.  There  was  a  storage  option  where  all  raw  data  could  be  stored  on  die  hard  disk  for  post-test 
analysis.  After  initialization,  the  system  established  radio  communication  with  the  tower.  When  this 
connection  was  established,  the  RDS  was  set.  As  a  round  flew  through  the  sensor  array,  the  RDS  was 
triggered  by  shock  wave  impinging  on  one  of  the  location  microphones.  The  two  high-fidelity 
microphones  captured  the  shock  wave.  The  RDS  computer  calculated  the  durations,  the  Cartesian  position 
of  the  round,  and  the  velocity  of  the  round. 

With  the  durations,  the  microphone  offset  distances  (derived  fi'om  the  location  information),  and  the 
velocity,  the  RDS  identified  a  round  type.  If  that  round  satisfied  the  kill  criteria  for  that  particular  target 
type  and  it  was  within  the  kill  zone,  the  RDS  sent  a  "target  down"  signal  to  the  ECU  of  the  RETS.  The 
RETS  lowered  the  target  as  if  it  received  a  target  down  signal  frcan  the  vibration  sensors.  The  round 
identification,  location,  and  velocity  information  was  transmitted  back  to  the  tower.  Another  PC  in  the 
tower  displayed  a  picture  identification  of  the  round,  an  "X"  on  a  tank  silhouette  target,  and  the  velocity 
(see  Figure  37).  The  RDS  transmitted  results  to  the  tower  each  time  the  system  was  triggered,  regardless 
if  the  kiU  criteria  was  satisfied.  It  took  2  s  for  the  RDS  to  process  and  transmit  data,  but  it  took  10-15  s 
if  the  storage  option  was  enabled.  During  the  processing  time,  any  other  incoming  round  would  be 
ignored.  The  RDS  brassboard  logged  the  vibration  input.  It  was  not  part  of  an  algorithmic  decision.  The 
RDS  was  PC  compatible.  This  enabled  the  designers  to  use  DOS  and  Windows  application  software  and 
hardware  to  develop  the  discrimination  algorithms  (Figure  36)  and  the  communication  system. 
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Figure  36.  Algorithm  flow  chart. 


8.2  Construction. 

8.2.1  Location  Interface  Board.  The  location  interface  board  (LIB)  was  a  6-in  x  6-in  printed  circuit 
board  (PCB)  that  incorporated  a  similar  circuit  designed  by  Dr.  Joel  Kalb  of  the  Human  Research  and 
Engineering  Directorate  (HRED)  for  his  round  location  system.  The  RDS’s  version  had  six  electret 
microphone  inputs  instead  of  four.  The  additional  inputs  were  the  right  front  and  left  front  microphones. 
The  microphone  front-end  electronics,  threshold  comparator,  debouncing  latches,  and  time  difference  gates 
were  identical  to  Dr.  Kalb’s  (Kalb,  unpublished). 
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Figure  37.  Tank  silhouette. 

The  LIB  provided  power  to  the  six-electrode  microphones  and  measured  the  time  delays  between 
location  microphones.  The  board  had  outputs  firom  five  time  difference  gates  and  four  direction  registers. 
The  gates  measured  the  time  delays  between  pairs  of  microphones:  front  left  and  front  center,  front  right 
and  front  center,  back  left  and  back  center,  back  right  and  back  center,  and  front  center  and  back  center. 
Each  gate  sent  a  pulse  when  one  of  the  two  microphones  was  triggered,  and  it  sent  another  pulse  when 
the  second  microphone  was  triggered.  Each  direction  register  associated  a  negative  or  a  positive  sign  to 
each  time  difference  gate  except  for  the  two  center  microphones.  It  was  assumed  that  the  shock  wave 
would  impinge  upon  the  front  row  of  the  array  before  the  back  row.  The  outputs  of  the  LIB  fed  directly 
into  the  PC/104  computer  through  the  counters  and  I/O’s  mini-module.  The  time  duration  and  direction 
of  fire  calculations  were  done  by  software. 

In  addition,  the  LIB  had  a  NPN  transistor  circuit  for  the  target  down  signal.  When  the  discrimination 
algorithm  determined  that  the  target  should  be  lowered,  the  computer  sent  a  5-ms  pulse  to  the  circuit  The 
circuit  acted  as  a  driver  and  sent  an  inverted  pulse  (HIT  SIG)  to  the  ECU  of  the  RETS  to  lower  the  target. 
Refer  to  LIB  circuit  diagrams  for  further  details  (Figure  38). 
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Figure  38.  Location  interface  board  circuits. 

8.2.2  Sensors  Interface  Board.  The  sensors  interface  board  (SIB)  was  also  a  6-in  x  6-in  PCB.  It 
contained  circuits  for  the  high-fidelity  microphones  (B&K’s  and  Countrymans),  the  vibration  sensors,  and 
additional  electronics  for  future  use  such  as  the  switches  and  LED’s.  Only  the  circuits  for  the  B&K 
microphones  and  the  vibration  sensors  were  used.  The  RDS  brassboard  used  two  high-fidelity  B&K 
microphones  to  capture  the  supersonic  shock  waves.  The  interfacing  circuit  was  a  simple  voltage  divider 
with  gain  of  1/4  and  impedance  buffer.  The  two  outputs  were  fed  directly  into  the  A/D’s. 

The  SIB  had  two  vibration  sensor  inputs.  The  vibration  circuit  was  identical  to  the  one  in  the  RETS 
ECU  except  for  the  hit  sensor  failure  circuit  (TM9-6920-442-14&P).  As  soon  as  one  of  the  two  vibration 
sensors  was  triggered,  the  circuit  produced  a  5-V,  1-10-ms  output,  signifying  a  tagged  hit  This  HIT 
signal  was  fed  into  the  computer  by  a  register  on  the  A/D  board.  Refer  to  SIB  circuit  diagrams  for  further 
details  (Figure  39). 
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Figure  39.  Sensor  interface  board  circuits. 


8.2.3  A/D  Board.  The  A/D  board  was  a  4-in  x  4-in  PCB  with  a  PC/104  bus  connector  that  stacks 
onto  the  PC/104  computer.  It  had  two  12-bit  analog-to-digital  converters  (ADC)  with  track  and  hold.  The 
ADC’s,  part  no.  MAX  120  from  Maxim  Integrated  Products,  ran  on  a  4-MHz  clock  and  sampled  at  a 
250-kHz  rate.  The  A/D  board  had  encoding  and  decoding  logic  that  enabled  the  data  bus  from  the  ADC’s 
to  interface  directty  with  the  PC/104  bus.  The  digitized  data  from  the  ADC’s  were  continuously  being 
dumped  on  the  bus.  The  computer  took  the  data  and  stored  them  in  a  circular  buffer.  When  one  of  the 
location  microphone  was  triggered  by  the  shock  wave,  the  BIGOR  gate  from  the  LIB  sent  a  pulse  to  the 
computer  through  a  register  on  the  A/D  board  to  collect  10,000  more  data  points.  After  the  10,000th 
point,  the  computer  stopped  the  ADC’s  from  dumping  anymore  data  onto  the  bus.  The  ADC’s  remained 
in  this  state  until  the  computer  finished  with  the  discrimination  algorithm.  The  A/D  board  also  accepted 
the  HIT  signal  from  the  SIB.  The  software  logged  the  HIT  signal,  but  it  did  not  use  the  signal  in  any 
decision  process.  Refer  to  the  A/D  board  circuit  diagrams  for  fumre  details  (Figure  40). 
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9.  RDS  INTEGRATION  WITH  RETS. 


9.1  Electronics  Packaging.  The  RDS  brassboard  system  architecture  is  shown  in  Figure  41.  The 
system  was  housed  in  a  24-in  x  20-in  x  12-in  stainless  steel  box  with  a  sealed  lid.  The  box  will  protect 
the  electronics  from  rain  and  fragment  debris.  On  the  outside  of  the  box,  there  were  six  location 
microphone  input  connectors  (MICl  to  MIC6),  two  high-fidelity  microphone  input  connectors  (B&Kl  and 
B&K2),  two  vibration  sensor  input  connectors  (VIBl  and  VIB2),  one  target  signal  down  output  cormector 
(HIT  SIG),  a  parallel  port,  a  serial  port,  a  VGA  monitor  input,  and  a  keyboard  input  Eight  40-ft-long, 
shielded,  twisted  pair  cables  connected  the  microphones  to  the  RDS  unit  A  20-ft  cable  connects  the  HIT 
SIG  output  to  the  RETS’s  ECU. 


The  core  of  the  RDS  brassboard  was  the  PC/104  embedded  computer.  The  computer  was  composed 
of  four  4-in  x  4-in  stackable  modules:  a  486-25SX  CPU  core  module  with  a  math  coprocessor,  a 
hard/floppy  disk  controller  mini -module,  a  SVGA  mini-module  and  a  counter,  and  a  I/O  s  mini-module. 
The  embedded  computer  only  required  a  single  +5-V,  750-mA  power  supply  for  operation.  The  PC/104 
embedded  architecture  was  chosen  because  it  is  compact  and  modular.  During  the  development  phase, 
a  more  sophisticated  computer  system  could  be  assembled  by  stacking  more  mini-modules.  A  mini¬ 
module  can  be  a  flash  memory  card,  a  data  acquisition  card,  an  I/O  card,  or  a  commumcation  card.  The 
RDS  brassboard  only  used  four  modules.  Only  the  486-25SX  core  module  and  the  counters  and  I/O’s 
mini-module  were  needed  for  operation  after  system  initialization. 
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9.2  ECU  Connection.  The  RDS  brassboard  interfaces  with  the  RETS  with  no  modification.  Due  to 
time,  budget  and  proprietary  RETS  hardware  constraints  it  was  impossible  to  incorporate  the  RDS  outputs 
into  the  RETS’s  communication  system  without  working  and  consulting  with  Paramax  System  Corp.,  the 
main  RETS  contractor.  Instead,  a  separate  radio  communication  subsystem  for  the  RDS  was  designed  and 
used.  In  the  future,  however,  the  RDS  outputs  would  be  incorporated  into  the  RETS’s  CCU,  and  the 
results  would  be  displayed  in  the  control  tower. 

The  RDS  interfaced  with  the  RETS  through  the  ECU  of  the  RETS.  It  hooked  up  to  the  RETS  by 
disconnecting  all  vibration  sensors  from  the  ECU.  Then,  cormecting  two  of  those  vibration  sensors  to  the 
RDS  vibration  sensor  inputs,  VIBl  and  VIB2,  and  connecting  the  RDS’s  target  down  signal  output,  HIT 
SIG,  to  one  of  the  hit  sensor  inputs  on  the  ECU,  the  hardware  connection  was  then  complete. 

9.3  Radio  Commnniratinn  to  Tower.  In  addition  to  the  program  written  for  the  RDS  hardware,  a 
second  program  was  written  for  the  RDS  system  to  interpret  results  sent  to  a  remote  site  from  the  target 
area.  As  mentioned  earlier,  the  target  area  and  remote  site  communicated  through  a  modem  and  hand-held 
radio  at  each  site.  The  RDS  computer  communicated  with  the  modem  through  a  serial  port.  The  modem 
was  a  Kantronics  modem  running  at  9,600  baud  while  the  radio  was  a  Motorola  GP300  VHP  radio  with 
an  operating  range  of  146  MHz  to  174  MHz.  The  radio  frequency  used  was  dependent  on  the  usable 
frequency  bands  at  each  test  site. 

10.  VERMCATIONA^ALIDATION  TEST 


10.1  Purpose.  The  primary  purpose  of  the  test  was  to  verify  the  performance  of  the  RDS  brassboard 
and  validate  the  system’s  operational  capabilities  within  the  spectrum  of  live-fire  conditions  expected  on 
a  U.S.  Army  training  range.  A  secondary  propose  was  to  allow  private  contractors  the  opportunity  to 
verify/validate  their  own  round  identification  systems  and/or  collect  data  that  may  assist  in  the  future 
development  of  RDS  for  use  on  U.S.  Army  training  ranges. 

The  verification  test  was  conducted  from  23  Jun-1  Jul  1993  at  the  high  velocity  range,  APG,  MD. 
A  RETS  lifter  was  used.  Two  6-ft  x  6-ft  x  2-ft  concrete  blocks  were  placed  to  the  left  of  the  target  to 
protect  the  RDS  brassboard  and  radio  communication  system.  The  RDS  was  connected  to  the  RETS  ECU. 
The  7th  ATC  training  environment  was  replicated  as  best  as  possible.  A  berm  was  built  around  the  lifter 
to  simulate  the  target  pit.  Discrimination  microphones  were  attached  to  1-m-long  stakes.  Location  array 
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microphones  were  attached  to  stakes  and  located  20  cm  above  ground  level.  Support  data  for  the  test 
included  Weibel  radar  velocity  data,  video  scoring  on  the  target,  local  meteorological  data,  and  survey 
data.  The  test  plan  called  for  all  training  projectile  types  fired  in  the  feasibility  and  data  collection  tests 
to  be  fired.  The  verification/validation  test  firing  matrix  and  test  results  are  presented  in  Appendix  K. 
Average  velocity  errors  of  3%  and  average  location  errors  of  9  cm  were  observed. 

10.2  Statistical  Summary  of  Data  at  APG.  Five  models  were  ccmstmcted  in  sections  4  and  5,  and 
each  was  evaluated  against  its  ability  to  predict  the  rounds  fired  during  the  validation  test.  The  five 
models  include  PSV  formulation  with  equation  25  (Whitham’s  theoretical  model),  the  three  empirical 
models,  and  the  Hybrid  model. 

The  three  empirical  models  and  the  Hybrid  model  that  used  discriminant  analysis  were  not  statistically 
optimal  since  the  assumption  of  equality  covariances  for  all  the  groups  was  not  met.  This  indicated  that 
the  estimate  of  pooled  covariance  matrix  was  a  biased  estimate  of  the  true  but  unknown  covariance  matrix. 
It  was  not  surprising  that  Model  HI,  with  two  total  durations  as  variables  in  the  model,  did  not  have  equal 
covariance  matrices,  since  the  covariances  were  not  equal  when  only  four  groups  were  considered  during 
the  feasibility  phase.  Attempting  to  achieve  this  for  seven  groups  to  discriminate  makes  it  even  harder 
for  all  the  models  considered.  Model  H  came  closest  to  having  equal  convariances.  However,  it  should 
be  noted  that  due  to  these  inequalities,  some  of  the  coefficients  may  be  unstable.  Nevertheless,  using 
pooled  estimates  of  the  covariance  matrix  for  each  analysis  resulted  in  robust  algorithms  for  most  of  the 
models,  based  on  the  results  from  this  phase  of  testing. 

Table  20  summarizes  the  number  of  rounds  correctly  classified  by  each  model.  Four  of  the  models 
performed  comparably  at  around  90%;  however,  empirical  Model  HI  performed  the  worse.  One  problem 
with  Model  HI  stemmed  from  the  microphone  array.  Two  total  durations  were  selected  from  the 
five  microphones  collecting  the  signatures.  This  change  in  configuration  made  the  coefficients  for  the 
algorithm  imreliable.  The  model  resulted  in  poor  prediction  especially  for  the  25-nun  rounds.  The  model 
predicted  a  12.7-mm  round  in  each  case. 

Model  I  used  data  from  all  the  microphones.  This  fact  made  it  appear  that  there  were  five  times  more 
independent  shots  than  there  actually  were.  The  five  duration  values  obtained  from  one  shot  were  not 
independent  outcomes,  and  the  velocities  paired  with  each  of  the  five  durations  were  the  same. 
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Table  20.  Model  Comparison— Number  of  Correct  aassifications 


Round 

No.  of 
Rounds 

Theoretical 

Empirical 

I 

II 

III 

Hybrid 

7.62  mm 

3 

3 

1 

1 

1 

3 

12.7  mm 

9 

9 

9 

9 

6 

9 

M910 

9 

9 

8 

9 

0 

9 

M793 

5 

5 

4 

3 

4 

TPGID 

15 

15 

15 

15 

14 

15 

M724/M865 

10 

8 

9 

9 

4 

9 

M490/M831 

11 

6 

10 

9 

5 

7 

Total 

62 

55 

56 

55 

30 

56 

Percent  Correct 

89% 

90% 

89% 

48% 

90% 

Model  II  has  the  greatest  appeal  of  the  three  empirical  models.  The  algorithm  was  developed  from 
the  6-m  microphone  which  is  suited  for  the  final  microi^one  array.  One  data  set  of  duration,  microphone 
distance,  and  velocity  corresponds  to  one  round  tested. 

The  theoretical  model  and  the  Hybrid  model  are  favored,  since  they  have  a  theoretical  basis  and 
provide  prediction  results  equal  to  or  greater  than  the  other  models.  Although  equation  25  was  derived 
from  smaller-caliber  rounds,  the  model  performed  quite  weU,  as  did  the  Hybrid.  Both  models  (as  well 
as  Models  I  and  II)  have  redundancy  incoiporated  in  them.  They  only  require  one  high-fidelity 
microphone  for  discrimination.  The  other  microphone  acts  as  a  check  and  as  a  backup  system  if  the  other 
microphone  malfimctions  during  training  exercises. 

Unusually  high  duration  signatures  measured  for  the  large-caliber  rounds  became  a  problem  during 
testing.  These  readings  were  at  least  50%  higher  than  expected,  and  they  did  not  occur  consistently  on 
the  same  microphone.  This  affected  Model  III  more  than  the  other  models  because  both  microphones  are 
required  for  round  classification.  It  was  hypothesized  that  the  large  durations  were  due  to  reflections  off 
the  concrete  block  protecting  the  instrumentation.  If  this  were  tme,  one  would  expect  the  effect  to  be 
greatest  for  the  microphone  closest  to  the  concrete  waU;  however,  this  did  not  explain  the  high  signal 
received  from  the  microphone  on  the  far  side  when  the  microphone  closest  to  the  wall  recorded  values 
as  expected.  The  source  of  the  reflection  is  unlikely  to  be  from  the  ground.  The  problem  was  not 
resolved  in  this  phase  of  the  program. 


11.  DEMONSTRATION  RESULTS  AT  7TH  ATC 

The  RDS  brassboard  and  a  backup  system  were  transported  to  Grafenwoehr,  Germany,  for  a  live  fire 
demonstration.  On  14  My  1993,  die  RDS  was  set  up  in  front  of  target  V-16  on  Range  1 17.  Target  V-16 
was  located  approximately  1,500  m  fiom  the  initial  firing  line.  This  target  was  chosen  in  order  to  get  the 
best  engagement  ranges  for  both  small-caliber  and  large-caliber  weapons.  The  microphones  and  tatget 
locations  were  surveyed  in  the  same  manner  as  in  previous  tests  (see  Figure  30).  Rve  different  projectile 
types  were  fired  fiom  an  MlAl,  an  MlAl  outfitted  with  TPGID,  and  an  M2.  Each  vehicle  was  fired 
from  a  random  location  on  the  range.  An  orange  kill  zone  was  painted  on  the  target.  Thermal  blankets 
were  placed  in  the  same  area  (see  Figure  42).  If  projectiles  were  determined  to  be  in  that  area,  the  RETS 
would  be  sent  information  to  lower  the  target.  During  the  pre-test,  9  out  of  12  projectiles  were  identified 
correctly  by  the  theoretical  model,  and  the  RETS  system  lowered  die  target  when  signaled  by  the  RDS. 
The  stored  data  was  analyzed  on-site  to  determine  why  the  three  projectiles  were  not  identified.  On-site 
data  reduction  revealed  that  acoustic  reflections  were  being  ct^tured  within  the  prescribed  window  being 
analyzed.  These  reflections  were  assumed  to  be  only  ground  reflections. 


Figure  42.  Kill  zone. 


These  reflections  sometimes  caused  an  error  in  the  total  duration  pulse  period.  A  longer  determination 
of  total  duration  pulse  period  caused  the  algorithm  to  incorrecdy  identify  the  projectile.  The 
predemonstration  test  matrix  and  results  are  located  in  Appendix  L. 
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In  an  effort  to  separate  any  ground  reflection  from  the  N-Wave,  the  target  was  raised  approximately 
0.5  m  for  the  demonstration  test  on  15  July  1993.  The  store  mode  was  turned  off  so  the  RDS  would 
compute  real-time  data  reduction.  The  MlAl  with  TPGID  fired  10  projectiles.  Only  3  of  10  TPGID 
rounds  were  identified  as  TPGID’s.  Similar  results  were  recorded  with  the  other  projectile  types  fired. 
Problems  were  also  encountered  with  the  RETS  interface.  The  RETS  shorted  during  the  demonstration 
and  did  not  respond  to  the  RDS  communication.  At  that  time,  it  was  not  known  why  the  RDS  perforaied 
so  poorly.  It  was  thought  that  ground  reflections  seen  from  the  previous  day  caused  the  problems. 

On  the  following  day,  the  RDS  was  set  up  again  with  the  store  mode  activated.  Many  projectiles  were 
fired  to  gather  data  in  the  field.  Due  to  rain,  most  of  the  data  were  not  useful.  Upon  our  return  to  APG, 
the  pretest  data  were  analyzed.  It  was  determined  that  reflections  were  coming  from  the  ground  and  from 
the  top  of  the  microphone  stake.  When  the  projectiles  were  fired  around  the  microphone  height  or  higher, 
the  reflections  were  very  close  in  time  to  the  original  N-Wave  or  even  superimposed  on  the  N-wave  (see 
Figure  43).  These  reflections  were  mathematically  determined  to  initiate  fi-om  the  flat  surface  on  top  of 
the  microphone  support  stake  (see  Figure  44). 

120  mm  HEAT  Round:  N-Wave  and  Reflection 


78 


Figure  44.  Microphone  support  stake. 

The  microphone  location  relative  to  the  stake  and  ground  surface  was  measured  during  this  test.  The 
shock’s  speed,  projectile’s  speed,  and,  therefore,  the  Mach  angle  were  also  measured.  The  location  of  the 
projectile  was  hand  scored  at  the  conclusion  of  the  test.  The  pressure  time  history  of  the  N-Waves  was 
recorded.  Therefore,  the  distance  that  the  shock  wave  traveled  over  time  could  be  calculated.  This  value 
is  represented  as  DISTl  and  DIST2  for  high-fidelity  microphones  1  and  2.  This  distance  was  compared 
with  the  actual  distance  the  shock  wave  would  have  traveled  if  reflected  finm  the  ground  or  the  stake. 
The  stake  reflection  represents  the  actual  distance  the  shock  wave  traveled  from  the  projectile  to  the 
microphone  using  the  stake  as  a  reflecting  surface.  By  comparing  the  calculated  reflection  distance  with 
those  due  to  the  ground  or  stake  it  can  be  shown  in  Table  21  that  at  least  five  of  the  projectiles  had  their 
N- Waves  reflected  off  of  the  stake.  Only  those  projectiles  fired  above  the  microphones  had  their  shock 
waves  reflected  off  of  the  stake.  Only  three  of  the  five  caused  algorithm  problems  when  computing 
N-Wave  duratioru 

As  a  consequence  of  raising  the  target  in  the  demonstration  test,  the  projectiles  were  fired  higher 
relative  to  the  ground  level.  This  separated  the  time  between  N-Wave  and  ground  reflections  but  caused 
the  problem  of  stake  reflections.  They  were  even  more  prevalent  than  during  the  pretest  and  resulted  in 
the  poor  demonstration  results.  The  problem  could  have  been  resolved  by  using  a  tripod  to  hold  the 
microphone  or  angling  the  top  of  the  stakes. 
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Table  21.  7th  ATC  Demonstration  Results  Analysis 


Round 

No. 

Mach  Angle 
(rad) 

n 

Time 

diff2 

(ps) 

Distl 

(m) 

Dist2 

(m) 

Hand  Scored 
Locations 

Ground 

refl 

(m) 

Stake 

refl 

(m) 

X 

(m) 

Z 

(m) 

1 

0.27 

4,092 

1.3 

0.552 

0.229 

1.2 

2 

0.37 

1,000 

1,000 

0.3 

0.3 

0.305 

0.686 

1.7 

0.3 

3 

0.48 

516 

476 

0.2 

0.1 

4 

0.61 

1,004 

1,316 

0.3 

0.4 

5 

0.58 

6 

0.38 

2,364 

0.7 

0.641 

1.842 

2.8 

1.2 

7 

0.38 

1,828 

1,616 

0.6 

0.5 

0.787 

1.334 

2.3 

0.7 

8 

0.31 

4,112 

3,900 

1.3 

1.2 

1.029 

0.203 

1.0 

9 

0.31 

3,724 

4,428 

1.2 

1.4 

0.457 

0.210 

1.2 

10 

0.66 

928 

0.2 

1.422 

0.734 

1.7 

0.3 

11 

0.56 

2,940 

2,228 

0.8 

0.6 

1.168 

1.524 

2.5 

0.5 

12 

0.60 

6,176 

1.7 

12.  CONCLUSIONS/RECOMMENDATIONS/FUTURE  PLANS 

The  goal  of  this  project  was  to  determine  if  round  discrimination  were  feasible  using  acoustic 
signatures.  If  it  were,  the  secondary  goal  was  to  build  a  RDS  to  perform  on  a  live-fiie  gunnery  range. 
The  first  goal  was  achieved  without  any  doubt.  The  models  created,  both  theoretical  and  empirical, 
worked  very  well.  The  sensors  used  were  excellent  laboratory  instruments.  The  data  recording  and 
analysis  equipment  was  well  chosen  for  this  endeavor  and  produced  an  excellent  database.  The  extent  of 
the  database  was  more  than  sufficient  for  the  analyses.  The  intellectual  contribution  from  outside  agencies 
was  tremendous  and  directly  impacted  the  success  of  the  program. 
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The  secondary  goal  was  only  partially  achieved.  Hampered  by  technical  problems,  the  final 
demonstration  at  the  7th  ATC  was  not  the  great  success  desired.  Although  the  RDS  performed  extremely 
well  on  the  last  days  of  the  validationA^erification  test,  the  weather  and  unforeseen  technical  problems 
precluded  a  favorable  result  at  Grafenwoehr.  This  can  be  most  likely  attributed  to  the  urgent  nature  of 
having  to  demonstrate  the  RDS  before  the  planned  development  date.  The  models  were  incorporated  into 
a  powerful,  flexible  system  that  allowed  for  easy  modification.  The  instrumentation  was  extremely  reliable 
laboratory  equipment  but  was  severely  affected  by  the  outdoor  environment.  Because  of  the  time  devoted 
to  hardware  performance  after  the  verification  test,  little  attention  was  given  to  the  reflection  jAienomenon 
discovered  during  the  verification  test  Ultimately,  this  problem  manifested  itself  into  a  significant  major 
technical  problem  during  the  final  demonstration. 

A  retrospective  analysis  sheds  new  light  on  these  problems.  Apparently,  the  shock  incident  to  the 
wood  microphone  stand  was  reflected.  The  reflection  was  actually  returning  to  the  microphone  before  the 
N-Wave  duration  was  completed.  In  differentiating  the  signal,  additional  peaks  were  encountered,  and 
incorrect  durations  were  calculated.  This  is  the  primary  reason  smaller  rounds  were  being  calculated  as 
large  tank  training  rounds.  The  RETS  equipment  failed  to  perfonn  due  to  the  electronic  short  in  the 
RETS,  affecting  the  testing  on  the  third  day.  If  these  small  problems  are  allowed  to  be  worked  out,  there 
is  a  high  degree  of  confidence  that  a  woiking  RDS  could  be  used  on  a  live-fire  gunnery  range. 

A  test  plan  has  been  developed  to  continue  the  RDS  development.  Future  modifications  include 
employing  different  microphone  setups.  A  possible  setup  involves  discrimination  microi^ones  that  reside 
on  ground  level  so  only  one  N-Wave  signature  will  be  analyzed.  Improvements  to  the  data  acquisition 
and  reduction  program  will  be  made.  The  N-Wave  period  can  be  better  determined  by  using  maximum 
and  minimum  peaks,  zero  crossings,  as  well  as  the  slopes.  Ruggedizing  the  microphone  hardware  will 
greatly  enhance  the  survivability  in  the  rain  and  snow.  These  concepts  wiU  be  field  tested  on  a  U.S.  Army 
training  range,  such  as  Fort  Knox  or  Fort  Stewart,  that  closely  resembles  the  terrain  of  Grafenwoehr.  The 
data  win  be  analj^ed  and  correlated  with  the  N-Wave  signatures  that  already  reside  in  the  data  base. 

For  a  future  demonstration  of  the  proof-of-concept,  the  same  RDS  brassboard  should  be  used  with  a 
new  high-fidelity  four-microphone  system.  Two  low-cost  high-fidelity  microphones  with  high  sensitivity 
should  be  used  to  process  the  N-Waves  of  small  roimds,  and  two  pressure  transducers  with  high-dynamic 
range  should  be  used  to  process  the  N-Waves  of  large  tank  rounds.  Both  microphones  and  transducers 
must  be  environment-resistant  and  mggedized.  The  current  RDS  hardware  must  be  modified  to 
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accommodate  two  additional  transducer  inputs.  The  signal  processing  algorithm  must  be  modified  to 
isolate  the  N-Wave  from  its  reflection.  The  discrimination  algorithm  should  incorporate  durations  from 
the  transducers.  Several  more  field  tests  are  needed  to  characterize  the  responses  of  the  pressure 
transducers.  There  is  a  high  degree  of  confidence  that  these  improvements  to  the  hardware  and  the 
algorithm  will  allow  the  RDS  to  perform  at  its  potential  and  permit  total  achievement  of  the  goals  of  this 
project. 
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APPENDIX  A: 

FEASIBILITY  TEST  RANGE  DATA 
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APPENDIX  B: 

PROJECTILE  CHARACTERISTICS  DATA  FROM  FEASIBILITY  TEST  PHASE 
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7611  735 


1^1 

617.0 

1133.01 

203637 

155|  M831|  3|  7775.91  2651.01  4100.0|  6561  369  |  429|  14541  238031  4154.7  192.802 


gj  4045.^  283li|  758!  589|  401]  1748 1  1451.6|  3105.8j  809.08 
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I  2581.1 

1537.7^ 

2059.4 

1282.91 

1547.6 

r  129!  717| 

393 

2554.01  274533 


203637 


7499.2I 

2651.01 

4000.01 

645!  3721 

413 

- ^ 

- 1 

61 

4045'Sl 

2427.2^ 

283131 

752i  6541 

1  333| 

1  1739 

[  1476.1 

1  3103.0 

1  809.08 

9521  742  291 


121  1893.8}  1547.61  1669.81  135  !  857|  250]  2142|  1048.6|  2166.0| 


157!  M831|  3!  8167.8 


274533 


203637 


192.802 


61  3640.9!  1618.21  26293|  7821  457!  548|  1787|  14233|  2927.8!  809.08 


274533 


1  2581.1 

I  1016.0 

i  2059.4 

r  9921  4501 

1  2158.6 

1  1282.9 

1  1669^ 
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400  meter  Data 


shot  i  Round  t  mic 

PosP 

NegPl 

Neg  P2 

A1  1  A2 

A3 

Atot 

Peak 

ABS 

cal 

i  type 

Pa 

Pa 

Pa 

msec  i  msec 

msec 

msec 

mPa-s 

mPa-s 

Pa/V 

3 

M724 

1 

909.8 

3200.0 

- 

107 

587 

- 

694 

52.6 

1646.0 

36539 

6 

2857.3 

1632.7 

- 

455 

483 

- 

938 

442.4 

925.1 

816.^7 

9 

1332.6 

1243.81 

> 

509 

544 

- 

1053 

3653 

723.0 

177.68 

HIH 

HP^ 

■S 

957.5 

- 

5421  592 

- 

1134 

586.7 

854.877 

j 

1 

4 

M724 

T 

913.5 

3197.2 

_ 

1231  5^ 

- 

700 

57.0 

1644.0 

36539 

_ 1 

6 

3061.3 

1836.8 

- 

4591  484 

• 

943 

4533 

942.0 

816.357 

- ^ 

1 

9 

1332.6 

1243.8 

• 

5061  554 

- 

1060 

3623 

731.4 

177.68 

1 

”l2 

1293.5 

1017.3 

- 

5431  604 

1147 

3013 

6772 

854.877 

. 

- 1 

6 

M490 

T 

36539 

— 

6 

5306.3 

3061.3 

- 

7311  897 

- 

1628 

1329.0 

2575.0 

816.357 

1 

9 

2785.3 

1915.5 

- 

853!  9^ 

- 

1839 

1132.8 

1990.8] 

177.68 

!■ 

Us 

1658.5 

1547.3 

- 

8971  1084 

1981 

8102 

16293 

854.877 

j 

- ! 

i 

7 

M490 

1 

! 

i 

36539 

t 

1 

6 

5102.2 

2857.3" 

- 

7271  1020 

- 

1747 

12673 

1  2523.0 

816.357 

j 

r~9 

3282.9 

3000.0 

- 

8391  980 

- 

1819 

12367 

1  2086.7 

177.68 

1 

1 

HIHl 

12 

1658.5 

1547.3 

- 

9061  1063 

- 

19691 

854.877 

1  ^ 

! 

8 

M490 

1 

i 

36539 

i 

6 

3061.^ 

- 

805I  ^ 

* 

1679 

1425.1 

2776.0 

816.357 

! 

9 

3260.6 

2872.2 

-| 

9271  968 

1 

1895 

1456.6 

2330.6 

177,68 

m 

1658.5 

1709.8 

- 

9731  1076 

- 

2049 

856.0 

1745.0 

854.877 

i 

1 

1  ■ 

i 

j  T 

, 

_ : _ 

- - 

- - 
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121  488.7 i  325.81  427.61  319!  23l|  1061  ffi6|  85.0i  165.41  203637 
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1000  meter  Data 


Atot 

Peak 

msec 

mPa-s 

1000  meter  Data 


shot  iRound  mic  1  Pos  P  I  Neg  Pi  j  NegP2  |  A1  |  A2  j  A3 


71|M910|  3|  976.51  -976.5 


— i 

9 

1  459.8 

1  ^.0! 

196 

!  214 1 

-I 

i  I 

S  121 

i  376.0 

-414.0! 

i - ^ 

2081 

238| 

-I 

6| 

i  635.8 

-555.41 

- 

159 

1 

1 

12! 

1  401.61 

j  -339.0i 

195| 

2471 

-! 

4421 

37.6| 

ABS 

cal 

mPa-s 

Pa/V 

1  74.5| 

177.3 

781.18 

367| 

1  5171 

118.4! 

1  730.78 

1  4101 

1  47.41 

95.^ 

1  333,15 

1  446i 

i  41.0 

1  93.3 

257.43 

1  -1380.01 

i _ H 

113] 

1621 

i _ d 

1  275| 

66.6| 

9 

!  439.8 

-353.1! 

i  q 

1961 

1  214 1 

i _ d 

!  4101 

1  46.5 

j  -693.31  -j 

147] 

1  182 

; - ij 

329| 

51.0 1 

61  781.91  -453.1] 


!  439.8 

-393.1! 

!  -! 

2001 

213| 

- 

I  4131 

1  46.9| 

82|M793|  1|  1284.01  -1037.01  -1421.0]  1711 


-504.21 

i  -599.2] 

2511 

1  2041 

115| 

1 

1  i 

9i 

1  499.7 

-353.1! 

1  -416.41 

3001 

1  2321 

lOOj 

-337i! 

!  -370.7] 

318] 

265] 

74| 

83 

M793I 

11 

1  1287.0 

-1119.01  -1421.01 

1541 

!  109] 

147 

1  4101 

1  135.5 

!  i 

61 

1  730.7] 

-504.21  -591.91 

262 1 

1  187| 

1  121: 

1  570 

1  106.5 

1  _22^ 

1  109 

1  626 

1  80.6| 

I  I  121  481.41  -339.81  -388.71  3001  272|  91!  663| 


84|M7931  11  1242.5!  -1036.61  -1421.0|  167 


61  679.6| _ -555.4!  -643.1i  259]  195]  116|  570} 


!  j  91  539.71  -419.2!  -479.7i  241}  99|  627] 


i 

1  i 

121 

1  499.41 

1  -337.21  -370.7! 

1  312, 

1  276 

1  77 

1  665 

1  71.9 

198.21  686.49 


117.11  730.78 


-306.3:  -1 

2111 

1  2291 

- 

1  4401 

39.8 

79.11 

1  257.431 

1  5701 

1  88.01 

207.5| 

1  730.78 

i  6321 

1  8Z5| 

1617 

1  333.15 

1  657] 

1  73.6] 

r  150.7] 

1  257.43 

435 

128.6 

1  343.0] 

1  686.49 

1  570 

1  93.2 

2097] 

1  730.78 
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1700  meter  Data 


Intentionally  left  blank. 
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APPENDIX  C: 

MICROPHONE  DATA  FROM  FEASIBILITY  TEST  PHASE 


107 


Intentionally  left  blank. 
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FEASIBILITY  TEST  HIGH  FIDELITY  MICROPHONE  POSITIONS: 


All  coordinates  are  in  meters. 
22SEP92 


Im  microphone 

x:  1700.104 

y:  0.157 

z:  5.694 

6m  microphone 

x;  1697.448 

y:  5.284 

z:  5.282 

12m  microphone 

x:  1696.780 

y:  11.641 

z:  5.109 

23SEP92 

Im  microphone 

x:  1700.104 

y:  0.157 

z:  5.523 

3m  microphone 

x:  1697.835 

y:  1.702 

z:  4.072 

6m  microphone 

x:  1697.448 

y:  5.284 

z:  5.118 

9m  microphone 

x:  1697.044 

y:  8.485 

z:  5.054 

12m  microphone 

x:  1696.780 

y:  11.641 

z;  4.944 

24SEP92 

Im  microphone 

x:  1700.104 

y:  0.157 

z:  5.548 

3m  microphone 

x:  1697.835 

y:  1.702 

z:  5.109 

6m  microphone 

x:  1697.448 

y:  5.284 

z:  5.139 

9m  microphone 

x:  1697.044 

y;  8.485 

z:  5.078 

12m  microphone 

x:1696.780 

y;  11.641 

z:  5.023 

25SEP92 

Im  microphone 

x:  1700.104 

y:  0.157 

z:  5.551 

3m  microphone 

x:  1697.835 

y:  1.702 

z:  5.112 

6m  microphone 

x:  1697.448 

y:  5.284 

z:  5.115 

9m  microphone 

x:  1697.044 

y;  8.485 

z:  5.075 

12m  microphone 

x:1696.780 

y:  11.641 

z:  5.039 

28SEP92 

Im  microphone 

x;  1700.104 

y;  0.157 

z:  5.547 

3m  microphone 

x:  1697.835 

y:  1.702 

z:  5.096 

6m  microphone 

x:  1697.448 

y:  5.284 

z:  5.139 

9m  microphone 

x;  1697.044 

y:  8.485 

z:  5.066 

12m  microphone 

x:  1696.780 

y:  11.641 

z:  5.014 

29SEP92 

Im  microphone 

x:  1700.104 

y:  0.157 

z:  5.541 

3m  microphone 

x:  1697.835 

y:  1.702 

z:  5.090 

6m  microphone 

x:  1697.448 

y:  5.284 

z:  5.121 

9m  microphone 

x;  1697.044 

y:  8.485 

z:  5.087 

12m  microphone 

x:  1696.780 

y:  11.641 

z:  5.008 

30SEP92 

Im  microphone 

x:  1700.104 

y:  0.157 

z;  5.547 

3m  microphone 

x:  1697.835 

y:  1.702 

z;  5.108 

6m  microphone 

x:  1697.448 

y:  5.284 

z:  5.136 

9m  microphone 

x:  1697.044 

y:  8.485 

z:  5.084 

12m  microphone 

x:  1696.780 

y:  11.641 

z:  5.026 
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FEASIBILITY  TEST  ROUND  LOCATION  MICROPHONE  POSITIONS! 


All  coordinates  are  in  meters. 
Array  1 


21SEP92 


sensor  1 

x:  1697.846 

y:  -4.458 

z:  3.385 

sensor  2 

x;  1697.151 

y:  0.128 

z:  4.691 

sensor  3 

x:  1696.569 

y:  4.278 

z:  5.226 

sensor  4 

x:  1699.552 

y;  -5.346 

z:  5.976 

sensor  5 

x:  1699.348 

y;  0.009 

z:  4.949 

sensor  6 

x:  1699.222 

y:  5.276 

z:  5.977 

Array  2 

28SEP92 

sensor  1 

x:  1697.852 

y:  -4.444 

z:  3.382 

sensor  2 

x:  1697.108 

y:  0.125 

z:  4.468 

sensor  3 

x:  1696.554 

y:  4.277 

z:  5.222 

sensor  4 

x:  1699.520 

y:  -5.323 

z:  5.969 

sensor  5 

x:  1699.475 

y:  0.056 

z:  4.901 

sensor  6 

x:  1699.238 

y;  5.273 

z:  5.972 

Change  dining  testing 
sensor  1 

Ax:  0.006 

Ay:  0.014 

Az:  -0.003 

sensor  2 

Ax:  -0.043 

Ay:  -0.003 

Az:  -0.223 

sensor  3 

Ax:  -0.015 

Ay:  -0.001 

Az:  -0.004 

sensor  4 

Ax:  -0.032 

Ay:  0.023 

Az:  0.002 

sensor  5 

Ax:  0.127 

Ay:  0.047 

Az:  -0.048 

sensor  6 

Ax:  0.016 

Ay:  -0.003 

Az:  -0.025 
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FEASIBILITY  TEST  Gl)N  POSITIONS: 


All  coordinates  are  in  meters. 
400  meter  Range 
22SEP92 


M60  x:  1299.993 
MlAl  x:  1299.666 

y:  24.205 
y:  24.539 

z:  3.735 
z:  3.558 

23SEP92 

M2  x:  1299.901 

y:  23.715 

z;  4.023 

1000  meter  Ranee 

28SEP92 

MlAl  x:  700.059 
M60  x:  700.017 

M2  x:  699.711 

y:  6.717 
y:9.122 
y:  9.160 

z:2.335 
x:  4.282 
x:  4.617 

1700  meter  Range 

29SEP92 

M60  x:  0.024 

MlAl  x;  0.033 

y:  1.739 
y:  2.299 

z:  4.007 
z:  3.723 

1400  meter  Range 

29SEP92 

M2  x:  299.572 

y;  4.877 

z:  4.510 

400  meter  Range 

30SEP92 

MlAl  x:  1299.639 
M2  x;  1300.034 

y:  24.031 
y:  24.603 

z:  3.489 
x:  4.012 
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85.03 

176.74  105.102 

Shotl 

Round  1 

Ch 

PosPI 

Neg  Pi! 

Neg  P2I 

A1| 

A2i 

A3 

ATot 

Peak 

ABS 

Cal 

Type 

F^l 

Pa 

F^i 

msec! 

mseci 

msec 

msec 

mPa-sec 

mPa-sec 

Pa/Volt 

75 

DM18 

T 

3501.951 

-3503.31 

140 

288! 

388 

516 

175.83 

1186.95 

700.68 

i 

2 

2411.68i 

-1853.21 

252 

!  3521 

476 

588 

244.21 

563.95 

554.026 

900.7681 


603.687 

483.884“ 


3505.51 


2772.621 


1139.65! 


780.66! 

-708.15 

356 

513.74! 

-509.85 

404 

[■  i 

1  ! 

1 

1 

3505.5 

-2454.4 

1 

2251.8 

-1621.9 

1 

1  748.71] 

-658.961 

364 

-900.74 


-573.97I 

-479.37! 


•3499.9I 


-2178.1! 


-996.891 


155.46 


119.12 


111.21 


244.29 


221.53 


319.416 


105.102 


548 

286.78 

645.11 

554.026 

684 

172.84 

380.291 

350.769 

756 

144.84 

299.88 

319.416 

816 

123.08 

249.02 

105.102 

- ’ 

— 

. 

552 

293.44 

725.2 

700.68 

1 - 1 

! 

608 

235.02 

525.04 

554.026 

867.45] 

-807.82 

]  332 

788| 

132.24 

1  286.73] 

319.416 

3505,51 


2772.61 


-426.191 


-3499.91 


-2470.il 


-1038.91 


793.75 

-703.351 

352 

519.44 

-524.791 

416 

313.76 


177.85 


146.31 


714.9^ 


394.76 


299.39 


554.026 


350.769 


319.416 


848 

127.89 

260.381 

105.102 

79l  DM181 


3505.5! 

2139.371 


864.651 


-3046.5! 

-1665.11 


-874.47! 


335.45 

839.39 

700.68 

239.62 

538.97 

1  554.026 

4! 

705.91 

-645.54 

372 

si 

501.65 

-524.77 

416 

_ 1_ 

3505 

-3500 

]  76 

2 

2772.6 

1 

-2211.9 

1  228 

1120.7 


872| 

5261 


82ITPL-TI  1  1  3503.31 

!  2770.11 


I  3 


904.581 


8!  526.5 


-2598.il 

-1332.3] 


-990.21 


-746.21 

-524.81 


-3503.31 

-2526.5! 


-755.741 


-526.51 


796 

133.93 

279.19 

319.416 

860i 

i  122.87 

253.84 

105.102 

554.026 


684 

173.92 

394 

350.769 

764 

151.85 

321.3 

319.416 

258.83  105.102 


181.721 


295.34 


995.49 


715.531 


393.72 


700.68 


554.026 


350.769 


i  556 

r — 

716 

144.46 

312.97 

319.416 

127.781 


259.6  105.102 


831  TPL-Ti 


2772.6 

1 

-2506.7 

i  212 

1147.37 

1 

-1109.1 

1  284 

464 

148.541 

1116.15 

700.68 

500 

326.49 

721.92 

554.68 

826.011 


1  526.031 

1  -524.771 

i  384 

319.416 


3503.3| 

2672.6! 


1052.41 


-846.48 


-458.15 


-2360.5 

-2114.2 


-1065.3 


808.761 

i  -331.24 

-753.28 

I  344 

525.51 1 

1  -525.41 

-525.41 

i  380 

228|  3161  448 
232!  3321  436 


3001  4201 _ 

4801 


I  532 

294.83 

685.26 

700.68 

544 

266.72 

610.92 

554.026 

374.4 

299.85' 


350.769 

319.416 


1  I  3505] 


21  2772.621 


31  1060.731 


-2988.41 


-2414.71 


684.471 


382.3 


554.026 


350.769 


793.751 

-761.81 

352! 

748 

143.52 

302.34 

319.416 

1  ROfi  77I 

-524.81 

4001  i 

816| 

128.3 

256.35 

105.102 
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Shotl  Round!  Chi  PosPi  Neg  PI 


Neq  P2I 


F^I 


•2469.91 


msec!  msecl 


Peak! 


mPa-sec 


307.64! 


257.71 


178.96 


152.14 


122.51 


ABSl 


mPa^sec 


703.97 


701.731 


316.88 


237.491 


Cai 


Pa/Volt 


700.68 


554.026 


319.416 


105.102 


1025.6 

700.68 

678.77 

554.026 

381.85 

350.769 

323.82 

319.416 

237.651  105.102 


1328.2 

700.68 

646.8 

554.026 

391.58 

350.769 

293.7 

319.416 

232.61 

105.102 

700.68 


554.026 


350.769 


700.68 


554.026 


350,769 


319.416 


105.102 


1  81 

1  485.5 

_ 

94!  TPL-T  1 1  3505.5 


I  2|  2511.9 


•2492 


•1900 


•1006 

•658.9 


j  220 

548 

1  236 

572 

257; 

590 

171.2 

380.8 

136.7 

281.8 

113.2 

226.8 

554.026 


350.769 


319.416 


105.102 


708 

776 

- 1 

1 

1  8 

!  439 

t 

t 

766.71 


59 


348.21 


274.al 


216.1 


350.769 


319.416 


105.102 


APPENDIX  F: 

MICROPHONE  DATA  FROM  COLLECTION  TEST  PHASE 
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Intentionally  left  blank. 


136 


All  coordinates  are  in  meters. 
30  Mar  93 


3-m  hub 

x: -1.512 

y:  1697.509 

z:  1.388 

6-m  hub 

x:  -5.022 

y:  1697.238 

z:  1.506 

9-m  hub 

x:-8.270 

y:  1697.001 

z:  1.450 

12-m  hub 

x:-l  1.432 

y:  1696.818 

z:  1.380 

Microphone  Elevations,  top  of  microphone: 


Hub  Reference 

31  Mar  93 

5  Apr  93 

6  Apr  93 

7  Apr  93 

3-m  hub 

2.723024 

2.723024 

2.726072 

2.738264 

6-m  hub 

2.86236 

2.853216 

2.86236 

2.856264 

9-m  hub 

2.80636 

2.80636 

2.812456 

2.818552 

12-m  hub 

2.7516 

2.711976 

2.699784 

2.6754 
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DATA  COLI.ErTTON  ROUND  I.OCATION  MICROPHONE  POSITIONS; 


All  coordinates  in  meters. 
30  MARCH93 


sensor  1 

x:  -.0242 

y:  1697.945 

z;1.516 

sensor  2 

x:  2.798 

y:  1698.901 

z:  1.685 

sensor  3 

x:  -0.192 

y:  1698.927 

z:  1.692 

sensor  4 

x:-3.230 

y:  1698.966 

z:  1.692 

sensor  5 

x:  -1.796 

y:  1698.902 

z:  1.692 

sensor  6 

x;  -2.205 

y:  1698.915 

z:  1.725 

8  APR  93 

10®  Rotation,  Counter  Clockwise  about  sensor  3 


sensor  1 

x:  -0.006 

y:  1697.976 

z:  1.654 

sensor  2 

x:  2.772 

y:  1699.465 

z:1.719 

sensor  3 

x:  -0.192 

y:  1698.927 

z:  1.692 

sensor  4 

x:-3.150 

y:  1698.448 

z:  1.607 

sensor  5 

x:  2.953 

y:  1698.460 

z:  1.698 

sensor  6 

x;  -2.974 

y:  1697.498 

z:  1.666 
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DATA  COLI.ECTTON  GUN  POSITIONS; 


All  coordinates  in  meters. 


References 


1700  M  Position 

x:0 

y:0 

z: 

0 

1000  M  Position 

x:  -7.181 

y:  700.252 

z: 

0.244 

800  M  Position 

x:  -8.530 

y;  900.249 

z: 

0.801 

600  M  Position 

x: -11.340 

y:  1100.137 

z: 

-0.175 

400  M  Position 

x:  -23.602 

y:  1299.785 

z: 

-0.534 

30  MAR  93 

7.62  mm  Muzzle 

x:  22.583 

y:  1299.923 

z: 

0.676 

.5  APR  93 

7.62  Muzzle 

x:  -6.296 

y:  900.175 

z: 

1.990 

6  APR  93 

.50  cal  Muzzle 

x:  -5.230 

y:  700.283 

z: 

1.490 

7.62  Muzzle 

x:  -22.074 

y:  1299.934 

z: 

0.737 

7  APR  93 

35  mm  Muzzle 

x:  -5.487 

y:  698.197 

z: 

2.051 

35  mm  Muzzle 

x:  -22.813 

y:  1298.062 

z: 

1.370 

35  mm  Muzzle 

x; -1.967 

y:  1.882 

z: 

1.778 

8  APR  93 

0.50  cal  Muzzle 

x:  -22.126 

y:  1299.823 

z: 

0.736 

35  nun  Muzzle 

x:  -22840 

y:  1297.942 

z: 

1.376 

9  APR  93 

35  mm  Muzzle 

x:  -1.900 

y: -1.877 

z: 

1.819 

35  mm  Muzzle 

x:  -8.713 

y:  698.340 

z: 

2.063 
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APPENDIX  H: 

ANGLE  CORRECTION  TECHNIQUE 
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A  method  to  theoretically  conect  the  velocity  is  to  employ  a  system  of  a  family  of  hyperbolas,  and 
solve  for  the  speed.  This  is  an  iterative  method,  since  more  than  one  of  the  variables  are  unknown. 
Referring  to  Figure  H-1,  the  following  is  defined: 


where 

Dgj  =  Distance  between  mic  1  and  PI. 
h  =  Elevation  of  projectile  above  ground  plane. 

M  =  Mach  Number. 

P  =  Azimuth  Deviation. 

yj  =  Equation  for  Shock  Front. 

y^  =  Equation  for  lead  Hyperbola. 

y^  =  Equation  for  succeeding  Hyperbola. 

y^  =  Equation  for  line  going  through  mic  1  and  mic  2. 

c  =  trae  distance  for  correct  velocity  calculatioa 


Figure  H-1.  Hyperbola  created  by  cone  intersecting  ground  plane. 

First  one  must  consider  that  the  conic  section  that  intersects  a  plane  that  is  parallel  to  the  conix  axis  will 
form  a  hyperbola  of  the  following  fonn: 
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The  angle  of  the  shock  cone  is  0  =  sin  h/M.  The  slope  of  the  shock  cone  can  be  written  as 

lllTlx_^oo  ■ 

Taking  the  derivative  of  the  hyperbolic  equation  yields 


dyh  _  dCbU’^) 


lim  _ 
x-»oo  dx 


The  distance  from  the  origin  to  where  the  hyperbola  intersects  the  x-axis  is  where  x^— a.  This  is  shown 


y  =  b 

a 


_-l  =  0, 
a 


2 

=  a. 
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This  is  the  point  where  the  shock  cone  generated  by  the  projectile  intersects  the  microphone,  as  depicted 
in  Figure  H-2. 


Figure  H-2.  Shock  cone  geometry. 
Using  Figure  H-2,  we  are  able  to  say  that 

£.  .  -I. 

h 


(33) 


This  implies  a  solution  for  the  magnitude  of  a  plus  solving  for  d, 

a  =  h^  (M^  -  1), 

_L  =  -  1, 

b  =  VaCM^  -  1)  =  h(M^  -  1).  (34) 

By  knowing  the  length  of  D,  since  it  intersects  the  trajectory  at  PI,  we  know  the  y  value  at  mic  1. 

yi  =  Dg-p  sin  p.  (35) 
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The  value  of  x  also  needs  to  be  solved  and  can  be  done  by  substituting  the  result  of  y  into  the  solution 
for  X  of  the  intersecting  hyperbolic  equation. 

Now  the  point  of  mic  1  is  known  in  Cartesian  coordinates.  Substituting  cot  P  for  into  the  equation 
that  describes  the  line  of  mic  1  and  mic  2,  we  can  define  the  intercept,  f2. 


yh  =  ^ 


N  a 


.-1 . 


Xu  =  a 


r  Y 

yh 


Xi  = 


r  Y 

Xh 


(36) 


f  =  ys  -^s  cot  p. 


(37) 


or 


f  =  yi  -  Xip , 

y*  =  yi  +  (^s  -  cot  p. 


(38) 


To  find  the  Cartesian  values  for  mic  2,  we  repeat  the  substitution,  but  for  the  second  hyperbola. 


y2  =  (DgT+lm)  sin  p , 
y2  =  yi  +(^2  cot  p , 

X2  =  (Y2  -Yj)  tan  P  +  Xj.  (39) 

Kjiowing  the  values  of  X2  and  y2  is  the  key  to  finding  c.  We  write  the  equation  for  the  hyperbola 
intersecting  mic  2  at  Dt  =  0.  This  is  identical  to  the  first  hyperbola  except  it  is  displaced  along  the  x-axis 
by  value  c.  Transforming  x’=x+c. 
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yh2  =  b 


(40) 


(Xh2  -  C)^ 
a 


-  1. 


Solving  the  equation  for  c, 

y2  =  (D  +  Im)  sin  p, 

X2  =  (y2  -  Yi)  tan  P  +  xj. 


- 

=  X, 

a 

yz 

+ 

r>l _ 

N 

In 

J 

(41) 


With  c  determined,  the  corrected  velocity  should  be  calculated  as 

V  =  _L.  (42) 

T2 

For  this  method  to  work,  an  initial  guess  must  be  made  for  the  Mach  number.  On  the  first  guess,  we 
assume  c  is  1  m.  For  an  increase  in  azimuth  deviation,  a,  the  true  velocity  length,  c,  is  decreasing.  As 
weU,  the  further  to  the  left  or  right  of  the  centerline,  c  becomes  even  smaller.  With  any  azimuth  deviation, 
the  predicted  velocity  wiU  always  be  higher.  The  resulting  c,  based  on  the  higher  Mach  number,  win  be 
calculated  to  slightly  less  than  reality.  The  resultant  updated  velocity  will  be  more  accurate. 

It  can  be  shown  geometrically  that  an  iterative  process  of  updating  the  velocity  and  Mach  number  and 
calculating  a  new  value  for  c  wiU  eventually  converge  to  the  true  velocity  rapidly.  For  the  purpose  of  the 
RDS,  one,  maybe  two,  iterations  would  be  aU  that  was  necessary.  The  true  test  of  this  method  lay  in 
processing  the  data  obtained  from  the  data  collection  test. 

Unfortunately,  the  method  did  not  work.  It  was  determined  that  the  fault  did  not  lay  in  the  theory  but 
in  the  assumptions.  It  was  assumed  that  the  azimuth  deviation  would  be  known.  The  method  of  finding 
the  azimuth  deviation  was  not  accurate.  The  solution  for  the  true  velocity  is  strongly  related  to  that  angle. 
As  a  result,  this  method  could  not  be  used. 
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APPENDIX  I: 

LOCATION  DATA  FROM  DATA  COLLECTION  TEST 
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RD.  IpREDICTED  impact  location  on  target _ Error  Analysis  Error  Analysis  PREDICTED  Error  Ana 

fjO.  X  1  Y  Z  delX  delZ  miss  distance  Velocity  _ del  V< 

Consec.  (m)  (m)  (m) _ m _ m _ m _ m/s _ 


ON  cs  o  S 

ON  Cn  ONilfi  r-»!ON  OO  On  GO 

O  CO  ^  ^  ^ 

d  d  o  did  oio  o  o  o 


Tf  lO  On  nO  !  I  lO  CO  nO  O  nO  ^  ^ 

D  ^  K  ^5  TT  t  VO  I  cs  00  CO  CO  Q  Q  DlDiQlQIQ 

CO  m  CO  CN }  CO  I  CO  cn  cs  on  co  ^  Z  2 

^  ^  tt  ^  CO 


r^iQlQ  Q 
2  2 


QOOOQOOOioOlOOOOOOOOGO 

^  ^  ^  ^ 


S  ^  ON  ONiON  ON  ON  ON  ON  ON  Q  Q  Q  QjQ  Q  Q  ON  Q  Q  ON  Q 


s  s  s 


^  ^  ^  ^ 


S  2  2  2  2!2  2  2 


f-^lcN  cO^  tn  volt^  oOlONi;:;;;;  ;^ 


ocioNO^lcvico'^trivot^ao 
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RP.  PREDICTED  IMPACT  LCXZATION  ON  TARGET  Error  Analysis  Error  Analysis  PREDICTED  Error  Analysis  EA 

NO.  X  Y  Z  deIX  del  Z  miss  distance  Velocity _ del  Vel  del  V 

Consec.  (m)  (m)  (m)  m  1  m  m  m/s  m/s  % 


QIQIQ  O  Q  Q  Q  S  o  S  S  S  oo  go  o  ^  ^  S  2 

5  CO  CO  VO  ^  VO  CN  ON  40  ;7  p  vo  in  10  10 

^  ^  ^  ^  ^  ^  CO  vd  CO  ^  Tf  rf  rr  rf  Tt 


vo  co  in  in  co  QQ  -- 
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d  d  fj*  9  9 


P  o  O  *7 
•7  d  'T  d 
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«  Z  o  I  i  I 


K  c5  ^  si^jg  s  s  s  s 


RD.  PREDICTED  IMPACT  LOCATION  ON  TARGET  Error  Analysis_  Error  Analysis  PREDICTED  Error  Analysis  EA  Array 

NO.  X _ Y _ Z _ del  X  del  Z  miss  distance  Velocity  delVel  delV  Config 

Consec. _ (n^ _ (m) _ (m) _ m _ m _ m _ m/s  m/s  % 

85  I  1700.823  |  -0,244  |  4.2^  o]i95  -0.027  |  0.197  1018  |  2  1 ,2,3,4 
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<N  M  9  9  cN  ^  9  m  9 


ON 

Tf 

1  c;  CO 

VO  i 

in 

ON 

nj 

ON 

tv 

CN 

m 

t  in 

1  CO 

o 

NO 

ON  ! 

tN 

OC 

00 

f— t 

ON 

CN 

CO 

i  OC 

1  IN 

IN 

NO 

NO  ! 

NO 

NO 

NO 

GO  1 

00 

ON 

t>s 

IN 

CN 

CN 

3 

NO 

CN 

CN 

CN 

CN 

ON 

ON  ! 

NO 

NO 

NO 

in 

ON 

o 

1  CO  t  1 

in 

NO 

NO 

ts 

ON  i 

i  NO  1 

i 

ON 

in 

00 
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00 

o 
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00 
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1 
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d 
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APPENDIX  J: 

VERMCATION/VALIDATION  TEST  DATA 
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1158.08 

1165.50 

922.08 

923.36 

1030.13 

1016.00 

834.90 

848.54 

849.62 

718.26 

870.70 

356.47 

848.18 

723.33 

724.77 

718.65 


1430.62 

951.25 

143471 

1454.55 

845.85 

1300.39 

123774 

1233.05 

702.74 

743.08 

773.69 

763.94 


669.46 
683.64 
730.99 
603.59 
533.12 
594.27 
‘364.10 
591.19 
597.55 

383.47 
364.37 

797.45" 

805.96 

125876 

1249.61 

863.93 

848.54 

863.93 


1154.36 

1150.06 

1134.46 


-0.349 

-0.039 

-0.042 

1.004 

-0.093 

-1.147 

0.556 

0.107 

-0.747 

0.347 

0.916 

-0.156 

-0.215 

0.363 


-0.216 

0715 

-0.548 

-0.708 

0.255 

0.446 

-0.489 

-0.624 

0.366 

-0.403 

0.599 

0.311 


0.445 

0.437 

0.762 

-0.432 

0.68 

-2.295 

-0.879 

-0.408 

-0,519 

-0.051 

0.647 


PREDICTED  ROUND  TYP 


Statistical 


Formulation 


PREDICTED 


HIT 

TARGET 

VELOCITY 

TARGET 

WENT 

ERROR 

KZONE 

DOWN 

% 

LOCATION 

ERROR 

(m) 


35  TPGID 

TPGID  KE 

6.405 

6.505 

35  TPGID 

TPGID  KE 

6.028 

6.307 

35  TPGID 

TPGID  HEAT 

6.475 

5.706 

35  TPGID 

TPGID  HEAT 

6.257 

6.100 

35  TPGID 

TPGID  KE 

6.293 

6-172 

35  TPGID 

TPGID  KE 

5,256 

7.150 

35  TPGID 

TPGID  HEAT 

6.297 

6.070 

35  TPGID 

TPGID  HEAT 

7.336 

5.089 

35  TPGID 

TPGID  HEAT 

5.764 

6.721 

35  TPGID 

TPGID  heat! 

6.090 

6.231 

35  TPGID 

TPGID  HEAT 

6.802 

5,269 

12.7 

TPGID  HEAT 

6.013 

7.098 

35  TPGID 

TPGID  HEAT 

5.205 

6.930 

35  TPGID 

TPGID  HEAT 

6.445 

6.118 

35  TPGID 

TPGID  HEAT 

6.270 

5.743 

35  TPGID 

LARGE  HE 

5,732 

6.349 

i-ir  -Lrr^ 

'5m 

^  71^ 

ARGE  HEA 
ARGE  HEA 
35TPGID 


LARGE  HE 
TPGID  KE 
TPGID  KE 
TPGID  HEAT 
LARGE  KE 
LARGE  HE 
TPGID  KE 
TPGID  HEAT 
LARGE  HE 
LARGE  HE 


7.62  OR  12.7 

5.614 

7.62  OR  12.7 

5.725 

7.62  OR  12.7 

5.584 

7.62  OR  12.7 

6.469 

7.62  OR  12.7 

6.436 

7.62  OR  12.7 

9.205 

7.62  OR  12.7 

7.959 

7.62  OR  12.7 

6.853 

7.62  OR  12.7 

6.588 

7.62  OR  12.7 

6.156 

7.62  OR  12.7 

5.536 

-1.670  1.000 

-2.109  1-501 

0.025  1.488 


7.62  OR  12.7 
7.62  OR  12.7 
7.62  OR  12.7 
7.62  OR  12.7 
7.62  OR  12.7 
7.62  OR  12.7 
7.62  OR  12.7 


7.62  OR  12.7  7.678 

7.62  OR  12.7  8.089 

7.62  OR  12.7  5.977 


1122,10 


629,76 

566.35 

1046.65 

1064.69 

1056.84 


955.96 

1444.60 

1426.14 


0.030 

1.474 

25  KE 

7.62  OR  12.7 

5.974 

6.003 

1.185 

0.943 

25  HE 

7.62  OR  12.7 

4.861 

7.163 

0.451 

1.504 

25  HE 

7.62  OR  12.7 

5.553 

6.424 

4- 

4- 

0.240 

2.592 

25  KE 

7.62  OR  12.7 

5.872 

6.332 

4- 

4- 

-0.770 

1.324 

25  KE 

7.62  OR  12.7 

6.768 

5.209 

-0.262 

3.091 

25  KE 

7.62  OR  12.7 

6.471 

5,965 

4- 

4- 

4> 

4- 

-0.325 

1.689 

LARGE  HEAl 

LARGE  HE 

6.329 

5.662 

4- 

4- 

4> 

* 

-0.124 

1.596 

LARGE  HEAl 

LARGE  HE 

6.118 

5.865 

* 

4- 

-0.201 

1.744 

LARGE  KE 

7.62  OR  12.7 

6.209 

5.786 

4- 

♦ 
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RD. 

# 


GUN 

POS. 

(m) 


ROUND 

TYPE 

(mm) 


DATE 


TIME 

H:M:S 


TERMINAL 

VELOCITY 

(m/sec) 


VIDEO  SCORING  LOCATION 


HORIZ 

(m) 


VERT 

(m) 


OFFSET 

(m) 


(I 


1 

800 

35  TPGID  KE 

6/23/93 

13:23:55 

1142.18 

2 

35TPGIDKE 

6/23/93 

13:29:54 

1156.61 

0.118 

3.012 

6.042 

6.276 

3 

35  OLD  HEAT 

6/23/93 

14:42:50 

901.42 

4).420 

2.584 

6.490 

5.664 

4 

35  NEW  HEAT 

6/23/93 

15:38:59 

908.96 

^.099 

3.056 

6.263 

6.075 

5 

1200 

35  TPGID  KE 

6/23/93 

16:iai7 

NA 

-^.090 

3.273 

6.308 

6.139 

6 

•• 

35  TPGID  KE 

6/23/93 

16:13:50 

1013,41 

0.951 

3.143 

5.270 

7.118 

7 

35  OLD  HEAT 

6/23/93 

16:21:56 

825.50 

^.142 

3.056 

6.305 

6.033 

8 

*• 

35  NEW  HEAT 

6/23/93 

16:25:51 

825.05 

-1.218 

3.164 

7.379 

5.028 

9 

1700 

35  TPGID  KE 

6/23/93 

17:09:16 

843.22 

0.465 

3.288 

5.777 

6.679 

10 

•• 

35  NEW  HEAT 

6/23/93 

17:13:16 

702.55 

0.049 

2.910 

6.087 

6.187 

11 

•• 

35  TPGID  KE 

6/23/93 

17:21:54 

839.64 

-0.810 

2.141 

6.829 

5.217 

12 

*• 

35  TPGID  KE 

6/23/93 

17:24:44 

NA 

6.218 

6.213 

13 

35  TPGID  KE 

6/23/93 

17:29:05 

848.95 

0.847 

2.388 

5.208 

6.891 

14 

35  NEW  HEAT 

6/23/93 

17:31:24 

707.04 

-0.220 

3.397 

6.466 

6.049 

15 

•• 

35  NEW  HEAT 

6/23/93 

17:40:09 

705.15 

-0.290 

1.699 

6.290 

5.711 

16 

35  NEW  HEAT 

6/23/93 

17:43:28 

708.16 

0.292 

2.148 

5.732 

6.315 

17 

800 

105  KE 

6/24/93 

10:58:24 

1382.3 

-0.285 

2.010 

6.296 

5.729 

18 

105  HEAT 

6/24/93 

11:44:24 

943.8 

0.146 

1.965 

5.863 

6.157 

19 

105  KE 

6/24/93 

11:51:11 

1382.9 

-0.593 

2.032 

6.605 

5.423 

20 

•• 

105  KE 

6/24/93 

11:54:37 

1387.2 

-0.787 

1.869 

6.791 

5.220 

21 

1200 

105  HEAT 

6/24/93 

13:30:44 

839.7 

0.199 

1.840 

5.805 

6.203 

22 

105  KE 

6/24/93 

13:35:08 

1300.2 

0.402 

1.847 

5.602 

6.406 

23 

1800 

105  KE 

6/24/93 

14:32:40 

1203.0 

6.218 

6.213 

24 

105  KE 

6/24/93 

15:23:19 

1196.7 

-0.619 

3.657 

6.922 

5.757 

25 

105  HEAT 

6/24/93 

15:29:07 

694.6 

0.411 

2.608 

5.674 

6.488 

26 

105  HEAT 

6/24/93 

15:39:13 

725.3 

6.218 

6.213 

27 

105  HEAT 

6/24/93 

15:45:08 

762.9 

6.218 

6.213 

28 

105  HEAT 

6/24/93 

15:48:55 

754.1 

0.340 

3.709  I 

6.029 

6.678 

29 

400 

12.7  n 

6/25/93 

12:47:07 

NA 

NA  ! 

NA 

NA 

NA 

!  30 

12.7 

6/25/93 

12:50:22 

NA 

NA 

NA 

NA 

NA 

31 

12.7 

6/25/93 

12:56:36 

NA 

NA 

NA 

NA 

NA 

32 

- 

7.62 

6/25/93 

13:10:13 

NA 

NA 

NA 

NA 

NA 

33 

800 

12.7 

6/25/93 

13:51:09 

NA 

NA 

NA 

NA 

NA 

34 

•• 

12.7 

6/25/93 

13:55:19 

NA 

NA 

NA 

NA 

NA 

33 

1. 

12,7 

6/25/93 

13:59:27 

NA 

NA 

NA 

NA 

NA 

36 

« 

12.7 

6/25/93 

14:01:44 

NA 

NA 

NA 

NA 

NA 

37 

12.7 

6/25/93 

14:03:15 

NA 

NA 

NA 

NA 

NA 

38 

M 

7.62 

6/25/93 

14:10.50 

NA 

NA 

NA 

NA 

NA 

39 

7.62 

6/25/93 

14:16:25 

NA 

NA 

NA 

NA 

NA 

40 

7.62 

6/25/93 

14:19:12 

NA 

NA 

NA 

NA 

NA 

42 

600 

25  HE 

7/1/93 

10:54:00 

NA 

NA 

NA 

NA 

NA 

43 

25  HE 

7/1/93 

10:55:00 

NA 

NA 

NA 

NA 

NA 

44 

M 

25  KE 

7/1/93 

11:15:00 

NA 

NA 

NA 

NA 

NA 

45 

25  KE 

7/1/93 

11:18:00 

NA 

NA 

NA 

NA 

NA 

46 

1700 

25  KE 

7/1/93 

11:30:00 

NA 

NA 

NA 

NA 

NA 

47 

25  KE 

7/1/93 

11:33:00 

NA 

NA 

NA 

NA 

NA 

48 

25  KE 

7/1/93 

11:43:00 

NA 

NA 

NA 

NA 

NA 

49 

•• 

25  HE 

7/1/93 

11:45:00 

NA 

NA 

NA 

NA 

NA 

50 

•• 

25  HE 

7/1/93 

11:46:00 

NA 

NA 

NA 

NA 

NA 

31 

25  HE 

7/1/93 

11:48:00 

NA 

NA 

NA 

NA 

NA 

52 

1000 

25  KE 

7/1/93 

14:15:00 

NA 

NA 

NA 

NA 

NA 

53 

25  KE 

7/1/93 

14:2a00 

NA 

NA 

NA 

NA 

NA 

54 

•• 

25  KE 

7/1/93 

14:22:00 

NA 

NA 

NA 

NA 

NA 
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25  KE 
25  HE 
25  HE 
25  HE 
25  KE 
25  KE 
25  KE 
25  HE 
120  HE 
120  KE 
120  KE 


7/1/93 

14:27:00 

7/1/93 

14:29:00 

7/1/93 

14:53:00 

7/1/93 

14:55:00 

7/1/93 

14:56:00 

7/1/93 

14:58:00 

7/1/93 

15:01:00 

7/1/93 

15:04:00 

7/1/93 

16:28:00 

7/1/93 

16:36:00 

7/1/93 

16:38:00 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 


NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA 

NA  ' 

NA 

NA 

NA 

NA 

NA 

NA 
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APPENDIX  K: 

STATISTICAL  ANALYSIS  OF  1,000-m  DATA 
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Intentionally  left  blank. 
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Implementing  Kmskal- Wallis  Test:  An  Example 

The  Kruskal-WaUis  nonparametric  test  was  implemented  in  the  feasibility  test  phase  for  the  1,000-m 
data.  An  example  of  the  test  procedure  is  provided  in  this  appendix.  The  null  hypothesis  (the  hypothesis 
that  is  to  be  rejected)  is  that  the  six  round  types  have  identical  distributions  (and  therefore  equal  means) 
vs.  the  alternative  that  at  least  two  round  types  differ. 

The  advantage  of  a  nonparametric  test  over  a  parametric  test  is  that  the  test  does  not  require  the 
assumptions  of  normality  of  the  data  or  homogeneity  of  variances  among  the  groups.  This  test  performs 
analysis  of  variance  on  the  ranks  of  the  data.  The  ranks  start  with  1  for  the  smallest  value  among  all  the 
six  round  types.  If  more  than  one  response  value  is  the  same  (that  is,  there  are  ties),  then  these  tied  values 
would  be  assigned  the  average  rank.  This  test  was  performed  on  each  variable  separately  and  at  four 
regions  corresponding  to  the  data  collected  from  the  microphones.  The  data  from  the  1-  and  3-m 
microphones  were  grouped  into  one  region  since  the  1-m  microphone  had  incomplete  data.  The  range  of 
offset  distances  in  meters  for  each  region  is  as  foUows: 

Region  1  ^ 

4<  Region  2  <7 
7<  Region  3  <10 
10<  Region  4. 

An  example  of  the  analysis  provided  for  this  test  is  for  the  positive  peak  pressures  collected  at  the 
microphone  stationed  12  m  from  die  center  of  the  target,  region  4. 


Table  K-1.  Example  of  Positive  Peak  Pressure  From  Region  4 


25  mm 

105  mm 

120  mm 

HE 

KE 

HEAT 

KE 

HEAT 

KE 

(Pa) 

(Pa) 

(Pa) 

(Pa) 

(Pa) 

(Pa) 

1495.70 

1109.00 

1508.30 

1088.00 

1733.00 

1109.50 

1800.00 

1109  JO 

404.20 

1735.00 

1112.00 

1987.00 

1145.00 
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Table  K-2.  Order  Values  and  Ranks 


Round 

(mm) 

Ordered  Values 
(Pa) 

Rank 

Rank  with  Ties 

376.00 

1 

1 

401.60 

2 

2 

25  KE 

404.20 

3 

3 

25  HEI 

481.40 

4 

4 

25  HEI 

483.98 

5 

5 

25  HEI 

499.41 

6 

6 

120  KE 

1088.00 

7 

7 

105  KE 

1109.00 

8 

8 

105  KE 

1109,50 

9 

93 

120  KE 

imso 

10 

93 

105  KE 

1112.00 

11 

11 

120  KE 

1145.00 

12 

12 

105  HEAT 

1495.70 

13 

13 

120  HEAT 

1508.30 

14 

14 

105  HEAT 

1733.00 

15 

15 

105  HEAT 

1735.00 

16 

16 

120  HEAT 

1800.00 

17 

17 

120  HEAT 

1987.00 

18 

18 

All  the  positive  peak  pressure  values  from  Table  K-1  are  ordered  and  then  ranked.  Table  K-2  lists 
the  ordered  values  and  the  corresponding  rank.  The  two  values  with  ranks  9  and  10  are  the  same  and 
therefore  should  be  given  the  same  rank.  The  fourth  column  provides  the  rank  correction  by  assigning 
the  average  of  the  two  tied  ranks,  9.5.  These  values  become  the  data  for  which  the  analysis  will  be  based 
(see  Table  K-3). 
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Table  K-3.  Ranks  of  Positive  Peak  Pressure 


25-mm 

25-mm 

105-mm 

105-mm 

120-mm 

120-mm 

M793 

M910 

M490 

M724 

M831 

M865 

4 

1 

8 

7 

5 

2 

9.5 

9.5 

6 

3 

16 

11 

12 

sum  of 
ranks 

15 

6 

44 

28.5 

49 

28.5 

The  test  statistic,  T,  for  the  Kruskal-Wallis  test  is  computed  the  following  way  (Conover  1980): 


T  = 


k=6 


i=l 


SR  2  N(N+1)2 


n; 


where 

SR  =  sum  of  the  ranks  as  given  in  Table  K-3,  and 


S^  = 


N-1 


f  N 

J^Ranki^  - 
Vi=l 


2  N(N+1)^ 


The  test  statistic,  is  compared  to  the  chi-square  distribution  with  5  (6  groups  - 1)  degrees  of  freedom 
at  the  a  =  0.05  level.  (The  chi-square  test  provided  an  approximate  quantile  since  the  exact  distribution 
was  not  available.) 

The  calculated  test  statistic  of  15.794  is  significant  at  the  alpha  =  0.01  level.  This  indicates  at  least 
two  groups  differ  in  their  average  rank.  To  determine  which  ones  differ,  pairwise  comparisons  were 
perfonned.  Two  groups,  i  and  j,  would  be  statistically  different  if  the  following  inequality  is  satisfied: 


SRj 

SRj 

>  t 

[52  N-l-Tl 

1 

1 

N-k 
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The  "t"  is  the  t-distribution  with  N  -  k  degrees  of  freedom  at  the  1  -  (x/2  quantile.  For  this  data  set,  if  two 
average  ranks  differ  more  than  3.01,  then  the  two  groups  are  significantly  different  at  the  a  =  0.05  level. 
The  ordered  average  rank  for  each  group  is  listed  in  Table  K-4. 


Table  K-4.  Ordered  by  Average  Rank 


25-mm  M910 

2.00 

25-mm  M793 

5.00 

105-mm  M724 

9.50 

120-mm  M865 

9.50 

105-mm  M490 

14.67 

120-mm  M831 

16.33 

To  summarize  the  results,  the  following  notation  is  used  for  this  example  and  for  the  other  N-Wave 
variables  in  Table  K-5.  The  round  types  are  listed  in  order  from  lowest  average  rank  to  the  highest. 
Parentheses  around  two  or  more  rotmd  types  indicate  there  is  no  significant  difference  between  (or  among) 
the  rounds.  The  short-hand  notation  to  summarize  the  results  for  positive  peak  pressures  from  region  4 
is: 


(25  K  =  25  H)  (105  H  =  120  H)  (105  K  =  120  K). 

To  conserve  space,  H  and  K  are  abbreviated  for  HEI  or  HEAT  and  KE,  respectively.  The  result  for 
example  teUs  us: 

-  25-mm  KE  is  not  significantly  different  from  25-mm  HEI  in  average  rank. 

-  25-mm  KE  and  HE  are  significantly  different  from  both  round  types  of  size  105  and  120  mm. 

-  105-rmn  HE  is  not  significantly  different  from  120-mm  HE. 

-  105-mm  KE  is  not  significantly  different  from  120-mm  KE. 

-  The  large  HEAT  rounds  (105  mm  and  120  mm)  are  significantly  different  from  the  large  KE  rounds 
(105  mm  and  120  mm). 
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Table  K-5.  Multiple  Comparison  Results 


Variable 

Region 

Notational  Summary 

(25  K  =  25  H)  ^  (105  K  =  120  K)  ^  (105  H  =  120  H) 

Positive  Peak 
Pressure 

(25  H  =  25  K)  (105  K  =  120  K)  ;£  (105  H  =  120  H) 

25  K  ^  25  H  (105  K  =  120  K)  9^  (105  H  =  120  H) 

(25  K  =  25  H)  (105  K  =  120  K)  (105  H  =  120  H) 

■1 

(25  K  =  25  H)  (105  H  =  105  K  =  120  H  =  120  K) 

Negative  Peak 
Pressure 

(25  H  =  25  K)  (105  K  =  [105  H)  =  120  H)  =  120  K] 

(25  K  =  25  H)  (105  K  =  105  H  =  120  H  =  120  K) 

wM 

(25  H  =  25  K)  ^  (105  K  =  120  K  =  105  H  =  120  H) 

H 

25  H  25  K  ?£  (105  K  =  120  K)  (105  H  =  120  H) 

Positive  Phase 
Duration 

25  H  25  K  105  K  120  K  ^  105  H  9^120  H 

25  H  25  K  105  K  9i:  120  K  *  105  H  120  H 

■■ 

25  H  25  K  105  K  120  K  105  H  120  H 

H 

(25  H  =  25  K)  (105  H  =  120  K  =  105  K  =  120  H) 

Negative  Riase 
Duration 

(25  H  =  25  K)  (105  K  =  [120  K  =  105  H)  =  120  H] 

(25  K  =  [25  H)  =  [105  K]  =  105  H  =  120  K  =  120  H} 

■■ 

(25  K  =  [25  H)  =  [105  K]  =  105  H  =  120  K  =  120  H) 

H 

(25  K  =  25  H)  (105  K  =  120  K)  ^  120  H  =  105  H) 

25  K  ^  25  H  105  K  ^  120  K  (105  H  =  120  H) 

Total  Duration 

mm 

25  K  25  H  (105  K  =  120  K)  *  (105  H  =  120  H) 

■■ 

25  K  25  H  ^105  K  5*  120  K  ;6  105  H  ^  120  H 

H 

25  K  25  H  (105  K  =  120  K)  9^  (105  H  =  120  H) 

Peak  Pressure 
Impulse 

25  K  ^  25  H  105  K  120  K  *  105  H  ^  120  H 

25  K  ^  25  H  ?6105  K  120  K  105  H  ^  120  H 

■■ 

25  K  ^  25  H  105  K  120  K  ^  105  H  120  H 

1 

25  K  25  H  ^  (105  K  =  120  K)  ^  (105  H  =  120  H) 

Peak  of  Absolute 

25  K  ^  25  H  105  K  120  K  ;tl05  H  ^120  H 

mm 

25  K  ^  25  H  105  K  120  K  105  H  120  H 

Pressure  Impulse 

mM 

25  K  25  H  ^  105  K  120  K  ^  105  H  120  H 

n 

25  K  25  H  (105  K  =  120  K)  (105  H  =  120  H) 

Negative  Phase 

25  K  25  H  105  K  120  K  ^  (120  H  =  105  H) 

25  K  25  H  ?£  105  K  91120  K  9t  105  H  9^  120  H 

Pressure  Impiilse 

25  K  ^  25  H  105  K  120  K  9t  (120  H  =  105  H) 

Most  of  the  variables  at  each  region  can  distinguish  between  large-  and  small-caliber  rounds  (25  mm 
vs.  105  mm  or  120  mm).  The  only  variable  that  cannot  is  negative  phase  duration  in  regions  3  and  4. 
The  results  of  the  pairwise  multiple  comparison  tests  for  each  region  on  each  variable  indicate  that  we  can 


not  discriminate  among  all  six  round  types.  A  general  conservative  result  is  the  relationship  as  seen  for 
the  positive  peak  pressure: 

(25  K  =  25  H)  ;£  (105  K  =  120  K)  *  (105  H  =  120). 

For  some  variables  in  some  of  the  regions,  the  two  small-caliber  rounds  were  distinguishable: 

25  K  ^  25  H  (105  K  =  120  K)  *  (105  H  =  120  H). 

For  the  most  part,  there  is  good  agreement  of  the  observed  values  with  the  statistical  nonparametric 
results.  However,  on  occasion,  there  are  a  few  results  that  do  not  appear  intuitive.  Take  for  instance, 
negative  phase  duration  plotted  against  offset  distance  in  Figure  K-1.  In  regions  three  and  four,  the 
25-mm  HE  and  KE  roimds  are  not  significantly  different  from  each  other  by  the  Kruskal-Wallis  test,  as 
indicated  by  the  parentheses  ( )  around  25  K  =  25  H,  in  Table  K-5.  Furthermore,  all  of  the  large-caliber 
rounds  are  not  significantly  different  from  each  other,  as  indicated  by  the  braces  {}.  These  two  results 
should  be  obvious  from  the  plot.  What  is  not  inmitive  is  the  nonparametric  test  indicating  that  25-mm 
HEI  and  105-mm  KE  are  not  significantly  different.  From  the  plot,  one  would  think  that  they  should  be 
different.  This  result  shows  that  by  analyzing  the  data  using  the  ranks,  large  differences  in  magnitude  of 
the  response  variable  in  the  original  units  may  be  masked. 

Conversely,  the  nonparametric  test  can  be  very  powerful  in  detecting  small  differences  in  magmtude 
if  the  ranks  of  the  data  within  a  round  type  are  close  together.  This  is  illustrated  in  Rgure  K-2  for  the 
peak  net  pressure  impulse.  The  smallest  rarrks  are  the  25  KE  rounds  followed  by  the  25  HEI.  Then,  in 
order  for  regions  1-4,  the  ranking  is  as  follows:  105  KE,  120  KE,  105  HEAT,  and  120  HEAT. 

What  we  have  shown  in  diese  analyses  is  that  many  of  the  variables  under  consideration  are  good 
candidates  to  discriminate  among  at  least  the  small-caliber  rounds  from  the  large-caliber  rounds.  In  some 
regions  for  some  variables,  all  the  rounds  are  significantly  different  from  each  other.  For  the  most  part, 
the  first  region,  representing  data  collected  closes  to  the  flight  of  the  rotmd  (<4  m),  is  the  hardest  region 
to  discriminate.  This  result  is  a  guide  to  the  placement  of  the  microphones  in  fumre  tests. 
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1 000  m  Firing  Position 


(ed)  uoj^ejnp  osegd  OAiieSou 
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1000  m  Firing  Position 


❖  +  □  X  <  * 


(s-BdUJ)  asinduj!  ajnssajd  lau  >jBad 
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APPENDIX  L: 

PREDEMONSTRATION  TEST  DATA 
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GUN 

POS. 

(m) 

ROUND 

TYPE 

(nun) 

DATE 

TIME 

H;M:S 

TERMINAL 

VELOCITY 

(m/s) 

HAND  SCORED  LOCATION 

OFFSET  DIS 
(m) 

TANCE 

(m) 

HORIZ 

(m) 

VERT 

(m) 

1 

2 

1410 

120(M865) 

7/14/93 

12:24:25 

NA 

mmm 

1.229 

120(M831) 

7/14/93 

12:26:56 

NA 

BSH 

1.686 

12.7  (M17) 

7/14/93 

12:31:16 

NA 

7.62  (M62) 

7/14/93 

12:32:27 

NA 

7.62  (M62) 

7/14/93 

12:32:41 

NA 

35  (DM18) 

7/14/93 

12:44:43 

NA 

0.641 

2.842 

^^Bip 

35  (DM18) 

7/14/93 

12:46:15 

NA 

0,787 

2.334 

1000 

25  (M910) 

7/14/93 

12:54:27 

NA 

1.029 

1.203 

1000 

25  (M910) 

7/14/93 

12:56:39 

NA 

0.457 

1.210 

1000 

25  (M793) 

7/14/93 

12:57:13 

NA 

1.422 

1.734 

1000 

25  (M793) 

7/14/93 

12:57:55 

NA 

1.168 

2.524 

1410 

35  (DM18) 

7/14/93 

12:58:49 

NA 

V 
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LIST  OF  SYMBOLS 


a 

A 

c 

D 

Di 

F(y) 

h 

K 

L 

M 

O 

P 

Po 

Rc 

Rl 

Rr 

T 

Ti 

T2 

T3 

T4 

T5 

V 

Vp 

X 

y 

yf 

yh 

yh2 

YlM2 

^LM5 

^LMi 

ys 

a 

P 

Y 

e 


ambient  speed  of  soimd,  meters/second 
microphone  separation  distance,  1  meter 
tme  distance  for  correct  velocity  calculation 
projectile  diameter,  meters 
arrival  time  difference,  i  =  1-4 

is  a  fimction  used  to  solve  the  differential  equation  determining  the  flow  past  the  body  of 
revolution 

elevation  of  projectile  above  ground  plane 
shape  constant,  1.82 
projectile  length,  meters 
location  microphone,  i  =  1-6 
Mach  number 

offset  distance,  perpendicular  distance  from  projectile  trajectory  to  mic,  meters 

overpressure 

ambient  pressure 

in-plane  distance  from  center  microphone  to  projectile  location,  meters 

in-plane  distance  from  left  microphone  to  projectile  location,  meters 

in-plane  distance  from  right  microphone  to  projectile  location,  meters 

total  duration  pulse  period,  seconds 

time  difference  between  microphone  2  and  1 

time  difference  between  microphone  3  and  2 

time  difference  between  microphone  5  and  2 

time  difference  between  microphone  5  and  4 

time  difference  between  microphone  6  and  5 

projectile  velocity,  meters/second 

predicted  projectile  velocity,  meters/second 

distance  along  the  axis  from  the  bow 

describes  a  characteristic  curve  relationship  between  x  and  O 

equation  for  shock  front 

equation  for  lead  hyperbola 

equation  for  succeeding  hyperbola 

y  position  of  microphone  2,  meters 

y  position  of  microphone  5,  meters 

y  position  of  microphone  i 

equation  for  line  going  through  mic  1  and  mic  2 

shock  angle,  ° 

azimuth  deviation 

ratio  of  specific  heats  of  air  (1.4) 

angle  from  microphone  array  axis  to 
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ATTN  TECH  LB 
STRBE  N  HEBERLEIN 
STRBE  NA  WEAVER 
STRBE  BLORE 
STRBE  CFLO 
STRBE  NE  3  CYS 
STRBE  ND  SPITZER 
STRBE  NDM  DILLON 
FT  BELVOm  VA  22060-5606 

14  DIR 

ATTN  AMSRL  SS  FF  J  GERBER 

TPHAM 

MFONG 

KTRAN 

AMSRL  WT  N  J  GWALTNEY 
J  INGRAM 
J  MCGARRITY 
AMSRL  WT  NA  R  KEHS 
AMSRL  WT  NB  M  ABE 
AMSRL  WT  ND  J  MILETTA 
AMSRL  WT  NF  L  JASPER 
AMSRL  WT  NG  T  OLDHAM 
AMSRL  WT  HD  J  CORRIGAN 
AMSRL  SL  C  J  HUGHES 
US  ARMY  RESEARCH  LABORATORY 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1145 
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3  DIR 

ATTN  TECH  LB 

AMSRL  MA  PA  W  HASKELL 

ERICAS 

US  ARMY  RESEARCH  LABORATORY 
WATERTOWN  MA  02172-0001 

1  CDR 

ATTN  STEW  NED  J  MEASON 
US  ARMY  WSMR 
WSMR  NM  88002 

2  CDR 

ATTN  TECH  LB 
AMFTA  RSA  T  MEITZLER 
US  ARMY  TACOM 
WARREN  MI  48397-5000 

3  CDR 

ATTN  TECH  LB 
GLONG 
D  BACH 

US  ARMY  NUC  AND  CHEML  AGCY 
7150  HELLER  LOOP  STE  101 
SPRINFIELD  VA  22150-3198 

2  COMMANDER 
ATTN  TECH  LB 
M  OGORZALEK 

US  ARMY  CONCEPTS  ANALYSIS  AGCY 
8120  WOODMONT  AVE 
BETHESDA  MD  20814 

1  CDR 

ATTN  TECH  LB 
US  ARMY  RSRCH  OFC 
PO  BOX  12211 
RSRCH  TRI  PK  NC  27709 

3  CDR 

ATTN  T  REEDER 
G  GOODWIN 
TECH  LIB 
US  ARMY  NGIC 
220  7TH  ST  NE 

CHARLOTTESVILLE  VA  22901-5396 


1  CDR 

ATTN  TECH  LB 

US  ARMY  DUGWAY  PROVING  GROUND 
DUGWAY  UT  84022 

2  US  ARMY  CORPS  OF  ENGRS 
ATTN  TECH  LB 

C  HERRING 

CONST  ENGRG  RSRCH  LAB 
PO  BOX  9005 

CHAMPAIGN  EL  61826-9005 

8  US  ARMY  ENGR  WATERY  EXPERMT  STAT 

ATTN  TECH  LB 
J STOUT 
J  INGRAM 
RDINAN 
G  MCMAHON 
PKING 
C  JOACHIM 
B  ARMSTRONG 
PO  BOX  631 

VICKSBURG  MS  39180-0631 

1  CMDT 

ATTN  TECH  LB 

US  ARMY  ENGINEER  SCHOOL 

FT  LEONARD  WOOD  MO  65473 

1  CDR 

ATTN  TECH  LB 

US  ARMY  COLD  REGIONS  R&E  LAB 
PO  BOX  282 
HANOVER  NH  03755 

2  CDR 

ATTN  STEYT  MT  EA  C  HASTON 
STEYT  MT  AT  A  HOOPER 
US  ARMY  YUMA  PROVING  GROUND 
YUMA  AZ  85365 

1  CDR 

ATTN  R  FRANSEEN 
US  AMCOM 

FIELD  ASST  IN  SCI  AND  TECHLGY  PROG 
FT  BELVOm  VA  22060-5606 
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1  US  ARMY  CORPS  OF  ENGINEERS 
ATTN  D  NEBUDA 

OMAHA  DISTRICT 
MAIL  CODE  CEMRO  ED  SH 
215  N  I7TH  ST 
OMAHA  NE  681024978 

8  OFFICER  IN  CHARGE 
ATTN  TECH  LIB 
FWARNOCK 
R  FERGUSON 
KCHIEN 
J  COLLINS 
CODE  R  45  R  PERSH 
M  RUPPALT 
CODE  R31  A  LE 

NAVAL  SURFACE  WARFARE  CTR 
10901  NEW  HAMPSHIRE  RD 
SILVER  SPRING  MD  20903-5640 

4  CDR 

ATTN  TECH  LIB 
L  FONTENOT 
J  BROWN 

CODE  G72  J  POWERS 
NAVAL  SURFACE  WARFARE  CTR 
17320  DAHLGREN  RD 
DAHLGREN  VA  22448-5000 

2  CDR 

ATTN  TECH  LB 
M  ORR 
CODE  1101 

NAVAL  SURFACE  WARFARE  CTR 
PHILADELPHIA  PA  19112-5083 

2  CDR 

ATTN  TECH  LB 
WILLIAM  A  SCHMIDT 
NAVAL  RSRCH  LAB 
WASHINGTON  DC  20375 

1  CDR 

ATTN  TECH  LB 

NAVAL  SURFACE  WARFARE  CTR 
SILVER  SPRING  MD  20903-5000 


1  CDR 

ATTN  TECH  LB 

NAVAL  SURFACE  WARFARE  CTR 
CHINA  LAKE  CA  93555-6001 

1  OFHCER  IN  CHARGE 
ATTN  TECH  LB 
NAVAL  EOD  FACILITY 
INDIAN  HEAD  MD  20640 

3  CDR 

ATTN  R  DENTON 
MAJCUTCHALL 
TECH  LB 

NAVAL  COASTAL  SYSTEMS  CTR 
PANAMA  CITY  FL  32407 

1  NAVAL  WARARE  ASSESSMENT  CTR 
ATTN  B  HULET 

PO  BOX  5000 
NWAC 

CORONA  CA  91718-5000 

2  DIR 

ATTN  TECH  LB 
J  BRIONES 

NUCLEAR  EFFECTS  DIRECTORATE 
WSMR  NM  88002 

2  CDR 

ATTN  D  VAUGHN 
TECH  LB 

MARINE  CORPS  RD&E  CMND 
QUANTICO  VA  22134-5080 

1  CDR 

ATTN  AETT  SA  MR  JEFF  SMITH 
7TH  ARMY  TRAINING  COMMAND 
HQ  7TH  ATC  UNIT  #28130 
APO  AE  09114 

1  CDR 

ATTN  AMXLS  SA  MR  SCOTT  KOHNKE 
US  NATIONAL  TRAINING  CTR 
BLDG  502 
FT  IRWIN  CA  92310 
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1  CDR 

ATTN  AMSTl  RLO  NTC  HECTOR  LOPEZ 

US  ARMY  STRICOM 

STRICOM  LNO 

BLDG  130  PO  BOX  10332 

FT  IRWIN  CA  92310 

1  CDR 

ATTN  AFZF  CS  SA  MR  ROY  J  HOLLEY 
HQ  III  CORPS 
FT  HOOD  TX  76544-5056 

1  DIR 

ATTN  AMSRL  SS  SL  JIM  CHOPAK 
US  ARMY  RESEARCH  LABORATORY 
2800  POWDER  MILL  RD 
ADELPHI  MD  20783-1197 

2  US  DEPT  OF  ENERGY 
ATTN  TECH  LB 

KK  22  K  SISSON 
WASHINGTON  DC  20585 

1  US  MILITARY  ACADEMY 
ATTN  CAPT  KEITH  E  MATTHEWS 
DEPT  OF  MATH  SCI 

WEST  POINT  NY  10996 

2  USIA  WORLD  NET 
ATTN  TECH  LB 
JRYAN 

RM  2410  PATRICK  HENRY  BLDG 
601  D  ST  NW 
WASHINGTON  DC  20547 

1  PHILLIPS  LAB 
ATTN  TECH  LB 
KIRTLAND  AFB  NM  87118-6008 

5  FIELD  COMMAND  DNA 
ATTN  LCDR  Z  MYERS 
CPT  M  SCOTT 
J  RENICK 

E  MARTINEZ  2  CYS 
KIRTLAND  AFB  NM  87115-5000 


3  WLMNME 
ATTN  G  PARSONS 
J  FOSTER 

TECH  LIB 

EGLIN  AFB  FL  32542-5000 

1  DIR 

NASA  SCIENCE  TECHL  INTNAT  FACLTY 

PO  BOX  8754 

BWI  AIRPORT  MD  21240 

2  NOAA  OCEANIC  AND 
ATMOSPHERIC  RSRCH 
ATTN  CPT  SMART 

W  CALLENDER 
1335  E  W  HWY 
SSMC3  MAILCODE  R  PDC 
SILVER  SPRING  MD  20910 

1  OB  CSD  CRB 

ATTN  A  M  JONES 
RM  1413  OHB 
WASHINGTON  DC  20505 

1  DB 

ATTN  TECH  LB 
IDAHO  NAT  LAB 
PO  BOX  1625 
IDAHO  FALLS  ID  83415 

1  DB 

ATTN  TECH  LB 
LOS  ALAMOS  NAT  LAB 
PO  BOX  1663 
LOS  ALAMOS  NM  87545 

4  DB 

ATTN  TECH  LB 
R  OSTENSEN 
S  SNYDER 
M  SAGARTZ 
SANDIA  NAT  LABS 
PO  BOX  5800 

ALBUQUERQUE  NM  87185 


1  AFOSR 

ATTN  TECH  LB 
BOLLING  AFB  DC  20332 
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2  DIR 

ATTN  TECH  LIB 
AKUHL 

LAWRENCE  LIVERMORE  NAT  LAB 
PO  BOX  808 
LIVERMORE  CA  94550 

1  AAI  CORPORATION 
ATTN  TECH  LB 
PO  BOX  126 

HUNT  VALLEY  MD  21030-0126 

2  ABERDEEN  RESEARCH  CENTER 
ATTN  J  KEEFER 

N  ETHRIDGE 
PO  BOX  548 
ABERDEEN  MD  21001 

1  AEROSPACE  CORPORATION 
ATTN  TECH  LB 
PO  BOX  92957 
LOS  ANGELES  CA  90009 

1  ALCAN  POWDERS  AND  CHEMICALS 
ATTN  TECH  LB 
PO  BOX  290 
ELIZABETH  NJ  07207 

3  ALLBD  CONTRACTORS  INC 
ATTN  TECH  LB 

T  CRAWFORD 
A  SIMPSON 
204  E  PRESTON  ST 
BALTIMORE  MD  21202 

1  ALUMINUM  COMPANY  OF  AMERICA 
ATTN  TECH  LB 
1501  ALCOA  BLDG 
PITTSBURG  PA  15219 

1  AMPAL 

ATTN  TECH  LB 
PO  BOX  31 

FLEMINGTON  NJ  08822 


2  ANSER  MISSILE  DIVISION 
ATTN  TECH  LB 
RLEGINUS 

CRYSTAL  GATEWAY  3 
1215  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22202 

1  APPLIED  RESEARCH  ASSOCIATES  INC 
ATTN  R  FLOREY 

2750  EISENHOWER  AVE  STE  104 
ALEXANDRIA  VA  22314 

2  APPLIED  RESEARCH  ASSOOATES  INC 
ATTN  R  GUICE 

RHEYMAN 

5941  S  MIDDLEFIELD  RD 
LITTLETON  CO  80123 

1  APPLIED  RESEARCH  ASSOQATES  INC 
ATTN  J  DRAKE 
3202  WISCONSIN  AVE 
VICKSBURG  MS  39180 

1  ATLANTIC  RESEARCH  CORPORATION 
ATTN  TECH  LB 
5390  CHEROKEE  AVE 
ALEXANDRIA  VA  22314 

1  ATOMIZED  METAL  POWDERS  INC 

ATTN  TECH  LB 
25  E  39TH  ST 
NEW  YORK  NY  10016 

1  BATTELE 

TWSTIAC 
505  KING  AVE 
COLUMBUS  OH  43202-2093 

1  BDM  CORPORATION 
ATTN  TECH  LB 

7915  JONES  BRANCH  DR 
MCLEAN  VA  22102 

2  BOEmG  HELICOPTERS 
ATTN  TECH  LB 

J  COSGROVE 

PO  BOX  16858  MS  P30  07 

PHILADELPHIA  PA  19142-0858 
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2  BOEING  MILITARY  AIRPLANE  COMPANY 
ATTN  TECH  LIB 
R  LORENZ 
PO  BOX  7730 
WICHITA  KS  67277-7730 

1  BOOZ  ALLEN  &  HAMILTON  INC 
ATTN  TECH  LIB 

CRTSTAL  SQ  2  STE  1100 
1725  JEFFERSON  DAVIS  HWY 
ARLINGTON  VA  22202-4158 

2  DENVER  RESEARCH  INSTITUTE 
ATTN  TECH  LIB 

L  BROWN 
PO  BOX  10127 
DENVER  CO  80210 

2  DYNAMIC  SCIENCE  INC 
ATTN  S  ZARDAS 
P  NEUMAN 
POBOXN 

ABERDEEN  MD  21001 

2  ELECTROSPACE  SYSTEMS  INC 
ATTN  TECH  LIB 
S  PAREKH 
PO  BOX  831359 
RICHARDSON  TX  75083-1359 

1  EMTEC  SYSTEMS  INC 
ATTN  J  LATTERY 
4500  ANAHEIM  AVE  NE  B6 
ALBUQUERQUE  NM  87113 

1  ETHYL  CORPORATION 

ATTN  TECH  LIB 
HOUSTON  PLANT 
BOX  472 

PASADENA  TX  77501 

1  FMC  CORPORATION 

ATTN  TECH  LIB 
1105  COLEMAN  AVE 
SAN  JOSE  CA  95108 

1  FRANKLIN  RESEARCH  CENTER 
ATTN  TECH  LIB 
BENJAMIN  FRANKLIN  PKY 
PHILADELPHIA  PA  19103 


2  GENERAL  SCIENCES  INC 
ATTN  TECH  LIB 
M  RILEY 

655  S  GRAVERS  RD 
PLYMOUTH  MEETING  PA  19462 

1  ITT  RESEARCH  INSTITUTE 
ATTN  TECH  LB 

10  W  35TH  ST 
CHICAGO  IL  60616 

2  ITT  AEROSPACE  OPTICAL  DIVISION 
ATTN  TECH  LB 

KRUSTER 
3700  E  PONTIAC  ST 
PO  BOX  3700 
FT  WAYNE  IN  46803 

3  THE  JOHNS  HOPKINS  UNIVERSITY 
ATTN  TECH  LB 

T  COUGHLIN 
J  KOUROUPIS 
APPLIED  PHYSICS  LAB 
JOHNS  HOPKINS  RD 
LAUREL  MD  20707 

1  KAMAN  SCIENCES  CORPORATION 
ATTN  F  MCMULLAN 

6400  UPTOWN  BLVD  STE  300E 
ALBUQUERQUE  NM  87110 

2  MARTIN  MARIETTA  AEROSPACE 
ATTN  TECH  LB 

M BAUER 
PO  BOX  179 
DENVER  CO  80201 

2  MCDONNELL  DOUGLAS  CORPORATION 
ATTN  TECH  LB 
C COREY 

BALLISTIC  MISSILE  DEFENSE 
5301  BOLSA  AVE 
HUNINGTON  BEACH  CA  92647 

1  MEDTHERM  CORPORATION 
ATTN  LARRY  JONES 
PO  BOX  412 
HUNTSVILLE  AL  35804 


NO.  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


2  NEW  MEXICO  ENGRG  RSRCH  INST 
ATTN  TECH  LIB 
R  ROBEY 
UOFNM 

ALBUQUERQUE  NM  87131-1376 

1  DEPARTMENT  OF  STATISTICS 
ATTN  DR  JAMES  R  THOMPSON 
RICE  UNIVERSITY 
HOUSTON  TX  77251-1892 

2  OLIN  ORDNANCE 
ATTN  TECH  LB 

J  KBIGER 

PRODUCT  MATERIAL  CONTROL 
10101  9TH  ST  N 
ST  PETERSBURG  FL  33716 

1  PCI 

ATTN  V  SCHMIDT 
900  19TH  ST  NW 
STE600 

WASHINGTON  DC  20006 

2  REYNOLDS  METALS  COMPANY 
ATTN  TECH  LB 
NKOOPMAN 

PLANT  #3 

4101  CAMP  GROUND  RD 
LOUISVILLE  KY  40211 

3  S  CUBED  MAXWELL  LABS 
ATTN  TECH  LB 

C  NEEDHAM 
K  SCHNIEDER 
ALBUQUERQUE  NM  87131 

1  SI  DIVISION  OF  SPECTRUM  39 
ATTN  W  SCHUMAN 

8831  SATYR  HILL  RD  STE  312 
BALTIMORE  MD  21234 

2  SCIENCE  APPLICATIONS  INTNL  CORP 
ATTN  TECH  LB 

J  BRYARS 

11526  SORRENTO  VALLEY  RD 
STE  A 

SAN  DBGO  CA  92121 


4  SCIENCE  APPLICATIONS  INTNL  CORP 
ATTN  J  SIMMONS 
J GUEST 
J  DISHON 
P  VERSTEEGAN 
PO  BOX  1303 
1710  GOODRICH  DR 
MCLEAN  VA  22102 

1  SCIENCE  APPLICATIONS  INTNL  CORP  ^ 
ATTN  S  DOERR 
2109  AIR  PK  RD  SE 
ALBUQUERQUE  NM  87106 

1  SBERLINE  MANUFACTURING  CO  INC 
ATTN  TECH  LB 

POBOX  A 

LANSFORD  PA  18232 

2  SIMULA  GOVERNMENT  PRODUCTS  INC 
ATTN  TECH  LB 

W  PERCIBALLI 
10016  S  51ST  ST 
PHOENIX  AZ  85044-5299 

1  SOUTHWEST  RESEARCH  INSTITUTE 
ATTN  TECH  LB 
PO  DRAWER  28510 
SAN  ANTONIO  TX  78284 

1  SRI  INTERNATIONAL 
ATTN  TECH  LB 

333  RAVENSWOOD  AVE 
MENLO  PK  CA  94025 

2  TECH  REPS  INC 
ATTN  F  MCMULLEN 
B  COLLINS 

5000  MARBLE  NE  STE  222 
ALBUQUERQUE  NM  87110 

1  TELEDYNE  MCCORMICK  SEIPH 
ATTN  C  GARRISON 
POBOX  6 

HOLLISTER  CA  95023-0006 

1  TERA 

ATTN  TECH  LB 
NM  nsrSTITUTE  OF  TECHLGY 
SOCORRO  NM  87801 
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1  THERMOGAGE  INC 

ATTN  CHARLES  BROOBtLY 
330  ALLEGHANY  ST 
FROSTBURG  MD  21532 

1  TRANSMET  CORPORATION 

ATTN  TECH  LIB 
4290  PERIMETER  DR 
COLUMBUS  OH  43228 

1  US  BRONZ  POWDERS  INC 

ATTN  TECH  LIB 
PO  BOX  31 
RT202 

FLEMINGTON  NJ  08822 

1  UNYSIS  CORPORATION 
ATTN  H  ZADEH 

A  J  FREDERICK 
4000  S  MEMORIAL  PKY 
HUNTSVE.LE  AL  35802 

2  UNIVERSITY  OF  MD  AT  COLLEGE  PK 
DEPT  OF  MECHL  ENGRG 

ATTN  ENME  J  WALLACE 
U  PIOMELLI 
RM  2168 

COLLEGE  PK  MD  20742-5121 

1  UNIVERSITY  OF  MD  AT  COLLEGE  PK 
ATTN  ENGRG  AND  PHYS  SCI  LIB 
MATHEMATICS  BLDG 
COLLEGE  PK  MD  20742 

1  VALIMET  INC 
ATTN  TECH  LIB 
431  E  SPERRY  RD 
STOCKTON  CA  95206 

5  WALCOFF  AND  ASSOCIATES  INC 
ATTN  T  TOYA 
C  PAQUETTE 
C  WALCOFF 
W  DEAL 
R  WALDEN 

12015  LEE  JACKSON  MEM  HWY 
#500 

FAIRFAX  VA  22033-3300 


1  STRICOM 

ATTN  AMSTI  MC  L  LAWRENCE 
2350  RESEARCH  PKY 
ORLANDO  FL  32826-3276 

ABERDEEN  PROVING  GROUND 

1  CMD.  AMCCOM 

ATTN:  SMCAR-ACW  (E3516) 
WEAPONS  SYSTEMS  CONCEPT  TEAM 

9  DER,  USAMSAA 

ATTN:  AMSXB-GB,  ABEL 

AMXSY-S,  CARROLL 
AMXSY-GC, 

B  BRAMWELL 
C  EISSNER 
L  MEREDITH 
W  WIDERMAN 
A  WONG 
LTC  HASSELL 
JTCG-ME,  LAGRANGE 

5  CMD,  USACSTA 
ATTN:  STECS-AE-CL, 

R.  scum 

D.  GRIFFEN 
STECS-AE-TL,  BINDEL 
STECS-AE-TH,  WILEY 
STECS-LI,  R.  HENDRICKSEN 

4  DIR,  ERDEC 

ATTN:  SCDRB-RTT, 

L  BICKFORD 
SFUNK 
G  GOLDSMITH 
I  SWANN 

1  CDR,  USAOC&S 
ATTN:  TECH  LIB 

6  CDR,  USATECOM 
ATTN:  AMSTE-SI-F 

AMSTE-CL 

AMSTE-EV-0 

AMSTE-TA-F 

AMSTE-ML 

AMSTE-TE-V 


NO  OF 

COPIES  ORGANIZATION 


NO.  OF 

COPIES  ORGANIZATION 


86  DIR,  USARL 

ATTN;  AMSRL-CI-S,  A.  MARK 
AMSRL-CP-TI,  J.  POLK 
AMSRL-SL, 

J.  WADE  (433) 

J.  SMITH  (433) 

AMSRL-SL-I,  M.  STARKS  (433) 
AMSRL-SL-B,  P.  DEITZ  (328) 
AMSRL-SL-BA,  J.  WALBERT  (1068) 
AMSRL-SL-BG,  A.  YOUNG  (433) 
AMSRL-SL-BL,  M.  RITONDO  (328) 
AMSRL-SL-BS,  D.  BELY  (328) 
AMSRL-SL-BV, 

J.  MORRISSEY  (247) 

W.  BAKER  (247) 

J.  COLLINS  (247) 

L.  MOSS  (247)  (5  CYS) 

W.  WINNER  (247) 

AMSRL-a-C, 

M.  TAYLOR 
AMSRL-WT-NC, 

R.  LOTTERO 
R.  LOUCKS  (20  CYS) 

P.  MULLER 
R.  THANE 
J.  SULLIVAN 

R.  RALEY 

G.  FERGUSON 

S.  SCHRAMEL 

C.  MERMAGEN 

A.  MIHALCIN 

AMSRL-WT-TB,  J.  CONDON 
AMSRL-WT-TD,  P.  KINGMAN 
AMSRL-WT-W,  C.  MURPHY 
AMSRL-WT-WA, 

B.  MOORE 

H.  ROGERS 
F.  BRANDON 

T.  GORDON  BROWN 
W.  D’AMICO 

B.  DAVIS  (20  CYS) 

D.  HEPNER 
M.  HOLLIS 
R.  MCGEE 

AMSRL-WT-WB,  A.  BARAN 
AMSRL-WT-WD,  G.  THOMSON 


AMSRL-WT-WE, 

R.  HAUSER 
W.  JOHNSON 
J.  THOMAS 
J.  WITAS 

AMSRL-WT-WG,  L.  PUCKETT 
AMSRL-W-WC,  T.  HAUG 
AMSRL-HR-SB,  W.  HANLON 
AMSRL-HR-SD,  J.  KALB 
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2  BUNDESAMT  FUR  WEHRTECHNIK 
UND  BESCHAFFUNG 
ATTN:  K  KOHLER 
TECH  LIB 
POSTFACH  7360 
5400  KOLENZ,  GERMANY 

2  DEFENSE  RESEARCH  ESTABLISHMENT 
SUFFIELD 
ATTN:  D.  RTTZEL 
TECH  LIB 
PO  BOX  4000 
MEDICINE  HAT 
ALBERTA,  TIA  8K6  CANADA 

4  DSTO,  MATERIALS  RESEARCH 
LABORATORY 
ATTN:  N.  BURMAN 
D.  SAUNDERS 
M.  BUCKLAND 
TECH  LIB 
PO  BOX  50 
ASCOT  VALE 

VICTORIA,  AUSTRALIA  3032 

2  INDUSTRIEANLAGEN- 

BETRIEBSGESGELLSCHAFT  MBH 
ATTN:  TECH  LIB 

H.  DIEKHOFF 
ABTEILUNG  FINITE 
BERECHNUNGSVERFAHREN 
EINSTEINSTRApE  20,  D-8012  OTTOBRUNN 

8  MINISTRY  OF  DEFENSE 

ATOMIC  WEAPONS  ESTABLISHMENT 
ATTN:  M.  GERMAN 
J.  THREAGOLD 
W.  BABBAGE 
B.  BOGARTZ 

I.  SMITH 

J.  TATE 
M.  KING 
TECH  LIB 

FOU124ESS,  ESSEX,  SS3  9XE,  UK 


4  MINISTRE  DE  LA  DEFENSE 

CENTRE  D’ETUDE  DE  GRAMAT 
ATTN:  S.  GRATIAS 
E.  CANTON 
D.  MERGNAT 
TECH  LIB 

46500  GRAMAT,  FRANCE 

2  MINISTERE  DE  L’EQUIPMENT 

LABORATOIRE  D’ESSAIS  D’EQUIPEMENTS 
D’ABRIS 

ATTN:  TECH  LIB 
D.  FAU 

BASE  DE  VmOULOU-ALVIGNAC 
46500  GRAMAT,  FRANCE 

2  NATIONAL  DEFENSE  RESEARCH  INSTITUTE 
WEAPONS  TECHNOLOGY  DEPARTMENT 
ATTN:  H.  AXELSSON 

TECH  LIB 
PO  BOX  98 

S-147  00  TUMBA,  SWEDEN 

3  NORWEGIAN  IHENSE  Oa^STRUCnOJ  SERVICE 
TEST  AND  DEVELOPMENT  SECTION 
ATTN:  TECH  LIB 

A.  JENSSEN 

OSLO  MIL/AKERSHUS 

OSLO,  NORWAY 

3  STORES  AND  CLOTHING  RESEARCH  AND 
DEVELOPMENT  ESTABLISHMENT 
ATTN:  S.  ELTON 
S.  CROSS 
TECH  LIB 
FLAGSTAFF  RD 
ESSEX,  C02  7SS,  UK 

3  wehrwissenschafiuchediensteiiefOr 

SPRENGMiriEL  UND  SONDERTECHNIK 

ATTN:  M.  KLAUS 

L.  KLUBERT 

TECH  LIB 

OBERJETIENBERG 

8230  SCHNIZLREUTH,  GERMANY 
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NO.  OF 

COPffiS  ORGANIZATION 


2  WEHRWISSENSCHAFTLICHEDIENSTELLE 
DER  BUNDESWEHR  FUR  ABC-SCHUTZ 
ATTN;  W.  REHMANN 
TECH  LIB 

HUMBOLTSTRApE/POSTFACH  1142 
3042  MONSTER,  GERMANY 

1  DEWEY  MCMILLAN  &  ASSOCIATES,  INC. 
ATTN:  J.  DEWEY 
1741  FELTHAM  RD 
VICTORIA,  BC  V8N  2A4,  CANADA 


206 


USER  EVALUATION  SHEET/CHANGE  OF  ADDRESS 


This  Laboratory  undertakes  a  continuing  effort  to  improve  the  quality  of  the  reports  it  publishes.  Your  comments/answers 
to  the  items/questions  below  will  aid  us  in  our  efforts. 

1.  ARL  Report  Number  ARL-TR-859 _ Date  of  Report  September  1995 _ 

2.  Date  Report  Received _ _ _ 

3.  Does  this  report  satisfy  a  need?  (Conunent  on  purpose,  related  project,  or  other  area  of  interest  for  which  the  report 

will  be  used.)  _ _ 


4.  Specifically,  how  is  the  report  being  used?  (Information  source,  design  data,  procedure,  source  of  ideas,  etc.) 


5.  Has  the  information  in  this  report  led  to  any  quantitative  savings  as  far  as  man-hours  or  dollars  saved,  operating  costs 
avoided,  or  efficiencies  achieved,  etc?  If  so,  please  elaborate. _ 


6.  General  Comments.  What  do  you  think  should  be  changed  to  improve  future  reports?  (Indicate  changes  to 
organization,  technical  content,  format,  etc.) _ _ 


Organization 


CURRENT  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 

7.  If  indicating  a  Change  of  Address  or  Address  Correction,  please  provide  the  Current  or  Correct  address  above  and  the 
Old  or  Incorrect  address  below. 


Organization 


OLD  Name 

ADDRESS  _ 

Street  or  P.O.  Box  No. 


City,  State,  Zip  Code 


(Remove  this  sheet,  fold  as  indicated,  tape  closed,  and  mail.) 
(DO  NOT  STAPLE) 


